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Abstract 

New techniques, both for generating and detecting turbulence in the helium su- 
perfluids 3 He-B and 4 He, have recently given insight in how turbulence is started, 
what the dissipation mechanisms are, and how turbulence decays when it appears 
as a transient state or when externally applied turbulent pumping is switched off. 
Important simplifications are obtained by using 3 He-B as working fluid, where the 
highly viscous normal component is practically always in a state of laminar flow, or 
by cooling 4 He to low temperatures where the normal fraction becomes vanishingly 
small. We describe recent studies from the low temperature regime, where mutual 
friction becomes small or practically vanishes. This allows us to elucidate the mech- 
anisms at work in quantum turbulence on approaching the zero temperature limit. 
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1. Introduction 



The transition to turbulence is the most well known example of all hydrody- 
namic transitions. It has been marveled for centuries, since dramatic demon- 
strations can be seen everywhere where a sudden change in the flow occurs, 
owing to a constriction in the flow geometry, for instance. For fifty years it has 



been known that turbulence also exists in superfluids (Vinen and Donnelly 



2007), although by its very nature a superfluid should be a dissipation-free 



system. In many situations it is found that on the macroscopic level super- 
fluid vortex dynamics mimics the responses of viscous hydrodynamics. This is 
one of the reasons why it has been thought that superfluid turbulence might 
provide a short cut to better understanding of turbulence in general. From 
the developments over the past fifty years we see that this has not become the 
case, superfluid turbulence is a complex phenomenon where experiments have 
often been clouded by other issues, especially by vortex formation and vortex 
pinning. Nevertheless, the topic is fascinating in its own right: When the flow 
velocity is increased, the inherently dissipation-free superfluid is observed to 
become dissipative and eventually turbulent. This is particularly intriguing 
in the zero temperature limit where the density of thermal excitations ap- 
proaches zero and vortex motion becomes undamped down to very short wave 
lengths (of the order of the vortex core diameter). 

There are two isotropic helium superfluids in which turbulence has been 
studied, namely the B phase of superfluid 3 He ( 3 He-B) and superfluid 4 He 
( 4 He II). In the anisotropic A phase of superfluid 3 He ( 3 He-A) dissipation is 
so large that conventional superfluid turbulence is not expected at the now 
accessible temperatures above 0.1 T c (Finne et al. 2003). Instead rapid dy- 
namics and large flow velocities promote in 3 He-A a transition in the topology 
and structure of the axially anisotropic superfluid order parameter field, a 



transition from linear line-like vortices to planar sheet-like vortices (Eltsov et 



aT||2002 ). Turbulence has also been studied in laser-cooled Bose-Einstein con- 



densed cold atom clouds, although so far only theoretically (Parker and Adams 



2005, Kobayashi and Tsubota 2008), but it is expected that experiments will 



soon follow. Here we are reviewing recent work on turbulence in rotating flow 
in both 3 He-B and 4 He II, emphasizing similarities in their macroscopic dy- 
namics. 

A number of developments have shed new light on superfluid turbulence. 
Much of this progress has been techniques driven in the sense that novel 
methods have been required, to make further advances in a field as complex 
as turbulence, where the available techniques both for generating and detecting 
the phenomenon are not ideal. Three developments will be discussed in this 
review, namely (i) the use of superfluid 3 He for studies in turbulence, which 
has made it possible to examine the influence of a different set of superfluid 
properties in addition to those of superfluid 4 He, (ii) the study of superfluid 
4 He in the zero temperature limit where the often present turbulence of the 
normal component does not complicate the analysis, and (iii) the use of better 
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numerical calculations for illustration and analysis. 

From the physics point of view, three major advances can be listed to 
emerge: In superfluid 3 He one can study the transition to turbulence as a 
function of the dissipation in vortex motion (Eltsov et al. 2006a), known as 
mutual friction. The dissipation arises from the interaction of thermal excita- 
tions with the superfluid vortex, when the vortex moves with respect to the 
normal component. In classical viscous flow such a transition to turbulence 
would conceptually correspond to one as a function of viscosity. This is a new 
aspect, for which we have to thank the 3 He-B Fermi superfluid where the easily 
accessible range of variation in mutual friction dissipation is much wider than 
in the more conventional 4 He II Bose superfluid. We are going to make use of 
this feature in Sec. 2, where we examine the onset of superfluid turbulence as 
a function of mutual friction dissipation (Finne et al. 2006a). 

Secondly, in Sec. 3 we characterize the total turbulent dissipation in super- 
fluid 3 He as a function of temperature, extracted from measurements of the 
propagation velocity of a turbulent vortex front (Eltsov et al. 2007). A par- 
ticular simplification in this context is the high value of viscosity of the 3 He 
normal component, which means that in practice the normal fraction always 
remains in a state of laminar flow. 

Finally, our third main topic in Sec. 4 are the results from recent ion trans- 
mission measurements in superfluid 4 He ( Walmsley et a/.|2007a Walmsley and 



Golov 2008a), where the decay of turbulence is recorded from 1.6 to 0.08 K. 



Here turbulent dissipation can be examined in the true zero temperature limit 
with no normal component. As a result we now know that turbulence and 
dissipation continue to exist at the very lowest temperatures. Although the 
dissipation mechanisms of 4 He or 3 He-B in the T — > limit are not yet firmly 
established (Vinen 2000, 2001, Kozik and Svistunov 2005b), it is anticipated 



that these questions will be resolved in the near future (Vinen 2006, Kozik 



and Svistunov 2008b). 



Phrased differently, our three studies address the questions (i) how turbu- 
lence starts from a seed vortex which is placed in applied vortex-free flow in 
the turbulent temperature regime (Sec. 2), (ii) how vortices expand into a 
region of vortex-free flow (Sec. 3), and (iii) how the vorticity decays when the 
external pumping is switched off (Sec. 4). The common feature of these three 
studies is the use of uniformly rotating flow for creating turbulence and for 
calibrating the detection of vorticity. Turbulence can be created in a superfluid 
in many different ways, but a steady state of constant rotation does generally 
not support turbulence. Nevertheless, at present rotation is the most practical 
means of applying flow in a controlled manner in the millikelvin temperature 
range. In this review we describe a few ways to study turbulence in a rotating 
refrigerator. Superfluid hydrodynamics supports different kinds of flow even 
in the zero temperature limit, so that turbulent losses can vary greatly both 
in form and in magnitude, but generally speaking, the relative importance of 
turbulent losses tends to increase with decreasing temperature. Two opposite 
extremes will be examined: highly polarized flow of superfluid 3 He-B, when 
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Fig. 1. Mutual friction parameter £ = (1 — a') /a as a function of temperature. 
In superfiuid dynamics this parameter, composed of the dissipative mutual friction 
a(T) and the reactive mutual friction a'(T) corresponds to the Reynolds number Re 
of viscous hydrodynamics. Typically, when Re > 1, turbulence becomes possible in 
the bulk volume between interacting evolving vortices. This transition to turbulence 
as a function of temperature can readily be observed in 3 He-B (at 0.59 T c ), while in 
4 He II it would be within ~ 0.01 K from the lambda temperature and has not been 
identified yet. 

a vortex front propagates along a rotating cylinder of circular cross section 
(Sec. 3), and the decay of a nearly homogeneous isotropic vortex tangle in su- 
perfiuid 4 He (Sec. 4), created by suddenly stopping the rotation of a container 
with square cross section. 

Turbulent flow in superfiuid 3 He-B and 4 He is generally described by the 
same two- fluid hydrodynamics of an inviscid superfiuid component with singly- 
quantized vortex lines and a viscous normal component. The two components 
interact via mutual friction. There are generic properties of turbulence that 
are expected to be common for both superfluids. However, there are also inter- 
esting differences which extend the range of the different dynamic phenomena 
which can be studied in the He superfluids: 

• In typical experiments with 3 He-B, unlike with 4 He, the mutual friction 
parameter a can be both greater and smaller than unity (Fig. 1) - this 
allows the study of the critical limit for the onset of turbulence at a ~ 1 
(Sec. 2); 

• The viscosity of the normal component in 3 He-B is four orders of magni- 
tude higher than in 4 He, hence the normal component in 3 He-B is rarely 
turbulent, which amounts to a major simplification at finite temperatures 
(but not in the T = limit with a vanishing normal component); 

• While the critical velocity v c for vortex nucleation is much smaller in 3 He-B, 
pinning on wall roughness is also weaker; this makes it possible to create 
vortex-free samples, which are instrumental in the transitional processes 
studied in Sees. 2 and 3; on the other hand, the ever-present remanent vor- 
tices in superfiuid 4 He are expected to ease the production of new vortices, 
which becomes important in such experiments as spin- up from rest; 
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• The vortex core diameter in both liquids is small (which allows us to use 
the model of one- dimensional line filaments), but in 3 He-B it is up to three 
orders of magnitude larger than in 4 He; hence the dissipation mechanisms 
in the T — > regime, which ultimately rely on the emission of excitations, 
are expected to work in 3 He-B at larger length scales and to lead to more 
significant energy loss in vortex reconnections. 

A comparison of the turbulent dynamics in these two superfluids allows 
one to identify generic properties that are common for both superfluids, and 
also to pinpoint specific reasons when there are differences. The main quantity 
controlling dissipation is the mutual friction dissipation a(T), which dominates 
the temperature dependence of the dynamic mutual friction parameter ( = 
(1 — a') /a, shown in Fig. 1. Experimental values are plotted with filled symbols 
for 3 He-B (Bevan et al. 1997) and with a solid line for 4 He (Donnelly and 



Barenghi 1998). At low temperatures, the following extrapolations are used 
(shown as dashed lines): For 3 He-B at a pressure P = 29 bar, we use a = 
37.21 exp(— 1.968T C /T), where the value for the superfluid gap A = 1.968T c /cb 
is a linear interpolation as a function of density p between the weak coupling 
value at zero pressure and that measured by Todoshchenko et al. (2002) at 
melting pressure. For 4 He at saturated vapour pressure (svp), we follow Kozik 
and Svistunov] (|2008b[) and use a = 25.3 exp(-8.5/T)T- 1 /2 + 5 78 . io^T 5 , 



where T is in K. 



2. Dynamic instability — precursor to turbulence 



2. 1 . Introduction 



In practice, superfluid flow remains dissipation-free only as long as there 
are no quantized vortices (or the existing vortices do not move, which is more 
difficult to arrange). The classic question in superfluid hydrodynamics is there- 
fore: How is the quantized vortex formed ( Feynman|1955 Vinen|196 3)? In flow 
measurements with bulk liquid the understanding about the origin of the first 
vortex has been improving in recent times. Whether it is created in an intrinsic 



nucleation process (Ruutu et al. 1997) or from remanent vortices (Solntsev et 



q/.|2007fl , which were created earlier in the presence of flow or while cooling the 
sample to the superfluid state (Hashimoto et al. 2007), these questions we are 



not going to address here. Instead we assume that the first vortex is already 
there, for instance as a remanent vortex. We then ask the question: How is 
turbulence started when the flow velocity is suddenly increased by external 
means? After all a turbulent vortex tangle is created through the interaction 
of many vortices: So how can turbulence start from a single seed vortex? 

In rotating 3 He-B one can create reliably a meta-stable state of vortex-free 
flow. It is then possible to inject a single vortex ring into the flow with neutron 
irradiation. When a slow thermal neutron undergoes a capture reaction in 
liquid 3 He with a 3 He nucleus, a vortex ring may escape from the overheated 
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Fig. 2. Principle of measurements on seed vortex injection at constant rotation Q 
and temperature T. Well separated isolated seed vortices are introduced in rotat- 
ing vortex-free counterflow. The initial high-energy state may then relax to the 
equilibrium vortex state via vortex generation processes which become possible at 
temperatures below the hydrodynamic transition at 0.59 T c . The Kelvin-wave in- 
stability of a single evolving seed vortex is the first step in this process. It is then 
followed by a turbulent burst which is started if the density of newly created vortices 
grows sufficiently. The combined process depends on the dynamic mutual friction 
parameter £ = (1 — a') /a which is shown on the top. On the very top the range of 
variation for this parameter in 4 He-II is indicated, i.e. the temperature regime of 
conventional 4 He-II measurements. 



reaction bubble into the flow if the flow velocity is above a critical value 
(Eltsov et al. 2005). Making use of this phenomenon one can inject a single 
vortex ring in vortex-free flow at different temperatures. At low temperatures 
it is observed that a turbulent vortex tangle is spontaneously formed from the 
injected ring (Finne et al. 2004a), while at high temperatures only a single 
vortex line results (Ruutu et al. 1998a). What is the explanation? 

This demonstration in 3 He-B shows that in addition to the applied flow 
velocity also mutual friction matters importantly in the formation of new 
vortices, in their expansion, and in the onset of turbulence. In 3 He-B, mutual 
friction dissipation a(T) is strongly temperature dependent (Fig. 1) and it so 
happens that ot(T) drops to sufficiently low value for the onset of turbulence in 
the middle of the accessible temperature range. The principle of seed vortex 
injection experiments is summarized in Fig. 2. In these measurements the 
number and configuration of the injected vortex loops can be varied. It turns 
out that the highest transition temperature is observed when turbulence starts 
in bulk volume from many small vortex loops in close proximity of each other. 
This transition has been found to be at T^ lk ~ 0.6 T c (Finne et al. 2004b) 
and to be independent of flow velocity over a range of velocities (3-6 mm/s). 
In viscous hydrodynamics the Reynolds number is denned as Re = UD/u c \, 
where U is the characteristic flow velocity, D the relevant length scale of the 
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Fig. 3. Vortex instability and turbulence in a rotating column of 3 He-B in the 
turbulent temperature regime, T < T^ lk . A seed vortex loop is injected in applied 
vortex-free flow and the subsequent evolution is depicted. Different transient states 
are traversed, until the stable rotating equilibrium vortex state is reached. 

flow geometry, and v c \ = rj/ ' p the kinematic viscosity. In an isotropic superfluid 
the equivalent of the Reynolds number proves to be £ = (1 — ac') /a. It defines 
the boundary between laminar and turbulent flow as a function of dissipation 



and is independent of flow velocity or geometry ilFinne et al. 2006b). 



However, if we inject instead of several closely packed vortex loops only one 
single seed loop in vortex-free flow (or several loops but so far apart that they 
do not immediately interact) then the transition to turbulence is found to be 
at a lower temperature and to depend on the flow velocity. Thus the onset of 
turbulence must also have a velocity dependent boundary. 

To explain all these observations, one has to assume that an independent 
precursor mechanism exists which creates more vortices from a seed vortex 
evolving in externally applied flow. The characterization of this instability 
is the topic of this section. It turns out that this can be done in 3 He-B in 
the temperature regime close to the onset T^ lk of turbulence in the bulk 
volume. Here the precursor often progresses sufficiently slowly so that it can 
be captured with present measuring techniques, while at lower temperatures 
turbulence starts too rapidly. This latter case is exactly what happens in 
4 He II: mutual friction dissipation is always so low in the usual experimental 
temperature range that the instability has not been explicitly identified. 

The central question is the reduced stability of the evolving seed vortex loop 
when mutual friction dissipation is decreasing on cooling to lower temperature: 
At sufficiently low a(T) an evolving vortex becomes unstable with respect to 
loop formation, so that one or more new vortex loops are split off, before the 
seed vortex has managed to reach its stable state as a rectilinear line parallel 
to the rotation axis. The evolution during this entire process, from injection to 
the final state, is depicted in Fig. 3 in a rotating cylindrical sample. The final 



2 



9 



4He-II 



0.14 0. 


65 1.8 


a(T) 


bulk turbulence exists 


vortex number con- 
served at Q = const 


► 


"rectilinear" vortices ^ n rcctilinc 
destabilized by £2 


ar vortices are stable 




evolving vortex onset regime of 

always unstable dynamic vortex 

generation 

i i i i i 


ideal regular 
vortex formation at v c 

i i i i 


T!T Q 
► 



0.2 04 6.6 o'.8 7c 



Fig. 4. Summary of vortex stability in rotating counterflow of 3 He-B, as a function 
of temperature. (Top row) The hydrodynamic transition at ~ 0.59 T c (at 

P = 29 bar pressure) separates regular and turbulent vortex dynamics. Above the 
transition vortices are stable in all situations which have been studied, while below 
turbulence becomes possible. (Middle row) In rotation at constant $7 rectilinear 
vortices are stable. In time-dependent rotation ^ 0) the "rectilinear" vortex 
turns out to be an idealization, presumably because of the experimentally inevitable 
slight misalignment between the rotation and the sample cylinder axes and because 
of surface interactions. In practice the "rectilinear" vortices are found to be stable 
above ~ 0.3 T c in time-dependent rotation, while at lower temperatures they tend 
to transform to increasingly turbulent configurations with increasing |0|. (Bottom 
row) Dynamically evolving vortices are stable above the transition, but at lower 
temperatures an evolving vortex may become unstable, generate a new vortex, and 
eventually bulk turbulence. The conditions at seed vortex injection determine the 
onset temperature T on below which turbulence follows. The onset temperatures 
have been found to concentrate in the regime 0.35 T c < T on < 0.59 T c . The very low 
temperatures below 0.3 T c display consistently turbulent response. 



state is the equilibrium vortex state, with an array of rectilinear vortex lines, 
where their areal density n v in the transverse plane is given by the rotation 
velocity Q: n v = 20/k. Here k = 2T[h/2m^ = 0.066 mm 2 /s is the superfluid 
circulation quantum of the condensate with Cooper-pairs of mass 2m^. In this 
equilibrium state at constant rotation the superfluid component is locked to 
solid-body rotation with the normal component, when averaged over lengths 
exceeding the inter- vortex distance t ~ 1/ ■Jtu,. In the ideal case all vortices 
are here rectilinear, while in the transient states in Fig. 3 vortices can exist in 
many different configurations. 

In Fig. 4 a rough classification is provided of the stability of vortices as a 
function of temperature (or more exactly mutual friction) in different config- 
urations and rotating situations. The lowest temperatures below 0.3 T c are in 
the focus of current research and have by now been probed with a few different 
types of measurements. The most extensive work has been performed by the 
Lancaster group. They create with various vibrating oscillators in a quiescent 
3 He-B bath a vortex tangle and then monitor the decay of the tangle with a 

. The total turbulent dissipation 



vibrating wire resonator (Bradley et a/.||2006 
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Fig. 5. Numerical calculation of the evolution of remanent vortices in rotating flow 
( de Graaf et al.\2007 ). (t < 0) Initial state with 22 vortices at 0.1 rad/s rotation. The 



vortices have been artificially tilted by 1 .4° , by displacing their end points uniformly 
by 1 mm at both end plates of the cylinder, to break cylindrical symmetry, (t = 0) 
Rotation is abruptly reduced to zero, to allow vortices to annihilate, (t = 600 s) 
After a waiting period At = 600 s at zero rotation, 12 remanent vortices remain in 
dynamic state. Rotation is then suddenly increased to Of = 0.5 rad/s. (t > 600 s) 
The 12 remnants start evolving towards rectilinear lines. This requires that the 
vortex ends on the cylindrical wall travel in spiral motion to the respective end 
plates. The parameters are: radius R = 3 mm and length h = 80 mm of cylinder, 
T = 0.4T C , P = 29.0bar, a = 0.18 and a 1 = 0.16 ( |Bevan et a^|1997| ). In the figure, 



the radial lengths have been expanded by two compared to axial distances. 

in a vortex front propagating along a rotating column (see Fig. 3) has re- 
cently been measured ( Eltsov et a/.||2007 ) and will be discussed in Sec. 3. Also 



the response of the superfluid component has been studied to rapid step-like 
changes in rotation, when Q is changed from one constant value to another. 
This type of measurement is commonly known as spin up or spin down of the 
superfluid fraction and will be extensively described for spin down in the case 
of superfluid 4 He later in this review. 



2.2. Seed vortex evolution in rotating column 

The motion of a seed vortex follows a distinctive pattern, while it expands in 
a rotating cylinder. Numerically calculated illustrations are shown in Figs. 5 
and 6 which depict the evolution of the seeds to stable rectilinear lines. In 
Fig. 5 an example with remanent vortices is examined, while in Fig. 6 the 




Fig. 6. Numerical calculation of the evolution of two curved peripheral vortices in 
an inclined rotating cylinder, when rotation is suddenly increased at t = from 



the equilibrium vortex state at ~ 0.03rad/s to f2f = 0.2rad/s (Hanninen et al. 



2007a). There are 22 vortices in this sample, of which two in the outermost ring 



(lying opposite to each other) have been initially bent to the cylindrical wall, to 
mimic an inclined cylinder. In the later snapshots at Of, the two short vortices 
expand towards the top and bottom end plates of the cylinder, to reach their final 
stable state as rectilinear lines. Parameters: R = 3mm, h = 30 mm, P = 29.0 bar, 
and T = 0.4 T c (which corresponds to a = 0.18 and a' = 0.16). 



initial configuration is an equilibrium vortex state in the usual situation that 
the rotation and cylinder axes are inclined by some small residual angle ~ 1°. 
These calculations (in a rotating circular cylinder with radius R and length h) 
describe the situation at intermediate temperatures when the vortex instabil- 
ity does not yet occur. The purpose is to focus on the motion of the expanding 
vortices and the transient configurations which thereby evolve. The character- 
istic property is the spiral trajectory of a vortex end along the cylindrical 
wall and the strong polarization on an average along the rotation axis. The 
calculations have been performed using the numerical techniques described in 
Sec. 2.5 Qde Graaf et aI1[2007| ). 

In Fig. 5 the remnants are obtained from an equilibrium vortex state ro- 
tating at O.lrad/s, by reducing rotation to zero in a sudden step-like decel- 
eration. The vortices are then allowed to annihilate at zero rotation for a 
period At. Some left-over remnants, which have not yet managed to annihi- 
late, still remain after this annihilation time. Ideally smooth walls are assumed 
without pinning. By suddenly increasing rotation from zero to a steady value 
Q{ = 0.5rad/s, the remaining remnants are forced to expand. The configu- 
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ration 100 s later shows how the spiral vortex motion has created a twisted 
vortex cluster in the center, with a few vortex ends still traveling in circular 
motion around the cluster. This motion thus winds the evolving vortex around 
the straighter vortices in the center. On the far right 200 s after the start of 
the expansion the cluster is almost completed. Nevertheless, this state is still 
evolving since ultimately also the helical twist relaxes to rectilinear lines, while 
the vortex ends slide along the end plates of the container. 

In Fig. 6 a calculation is presented with 20 rectilinear vortex lines and 
two short vortices which connect at one end to the cylindrical wall. This con- 
figuration mimics the equilibrium vortex state in a real rotating experiment 
where there exists some residual misalignment between the rotation and sam- 
ple cylinder axes. Depending on the angle of misalignment and the angular 
velocity of rotation some of the peripheral vortices may then end on the 
cylindrical side wall in the equilibrium vortex state, as shown on the far left. 
At t — rotation is increased in step-like manner from f2; to a higher value 
f2f. Two types of vortex motion are started by the rotation increase. First, 
the N = 20 rectilinear vortices are compressed to a central cluster with an 
areal density n v = 20f/« by the surrounding azimuthally flowing counterflow. 
Outside the vortex cluster the counterflow has the velocity 

v(Qf, r, N) = v n - v s = Q f r - KN/(2irr) , (1) 

The normal excitations are in solid-body rotation and thus v n = Qfr, while 
the superflow velocity around a cluster of kN circulation quanta is decaying 



as v s = KN/(2irr), where r > R v and the cluster radius R v ~ R\JVti/Vti. 
Experimentally it is convenient to define the number of vortices N via the 
initial equilibrium vortex state: N = N eq (fli), where the externally adjusted 
rotation velocity f2; ~ kN cci /(2ttR 2 ) uniquely defines A eq in a given exper- 
imental setup. It is customary to denote this specially prepared calibration 
value with Q v = The maximum counterflow velocity at fif at the cylinder 
wall is then given by v ~ (f2f — Q V )R. This definition has been used in Fig. 7 
to characterize the number of vortices in the central cluster. 

The second type of vortex motion, which is enforced by the increased rota- 
tion in Fig. 6, is the spiral motion of the two short vortices, as they become 
mobile and start expanding towards the top and bottom end plates, respec- 
tively. Let us now examine this motion in more detail. 

A vortex moving with respect to the superfluid component is subject to the 



influence from the Magnus lift force. This force can be written as (Donnelly 



1991) (p s is the density of the superfluid fraction) 

fM = p s KSX (v L - V s ) , (2) 

which acts on a vortex element s(£, t) with a tangent unit vector s = ds/d^ 
moving with velocity vl = ds/dt with respect to the superfluid component, 
which locally has the velocity v s . The motion from the Magnus force is damped 
by mutual friction which arises when the vortex moves with respect to the 
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surrounding cloud of normal excitations: 

fmf = -10PsK S X [§ X (v n ~ V L )] + 7oPs K S X (v n - V L ) . (3) 

The mutual friction force has dissipative and reactive components, which here 
are expressed in terms of the two parameters 70 and 7^. Balancing the two 
hydrodynamic forces, fvi+fmf = 0, one gets the equation of motion for the vor- 
tex line element, which when expressed in terms of the superfluid counterflow 
velocity v = v n — v s has the form 

v L = v s + as x (v n - v s ) - as x [s x (v n - v s )] . (4) 

Here the dissipative and reactive mutual friction coefficients a and a' appear. 
Conversion formulae between different sets of friction parameters are listed by 



Donnelly ( 1991 ). Evidently solutions of the equation of motion can be classified 
according to the ratio of the two components. The important parameter proves 
to be ( = (1 — a') /a (Fig. 1), which is the equivalent of the Reynolds number 



of viscous fluid flow (Finne et al. 2003). 



Two elementary examples are useful to consider. A single rectilinear vortex 
in rotating counterflow at w = Or moves such that its velocity components 
in the transverse plane consist of the radially oriented dissipative part aflr 
and the azimuthally oriented reactive part —(1 — a')Qr, when expressed in 
the rotating coordinate frame. The former is responsible for the contraction of 
the rectilinear vortices to a central cluster in Figs. 5 and 6. The latter causes 
the rectilinear vortex to rotate with the azimuthal flow with respect to the 
cylindrical wall. 

The second simple consideration concerns the end point motion of an evolv- 
ing vortex along the cylindrical wall. Since the vortex end is perpendicular 
to the cylindrical wall, it has from Eq. (4) a longitudinal velocity v^ z = 
av(Qf, R, N) and an azimuthal component = —(1 — a')v(Q{, R, N). Ev- 
idently other parts of the vortex also contribute to its motion, in particular 
its curvature where it connects to the cylindrical wall. However, it turns out 
that the end point velocity is an approximate guide for the expansion of a 
single vortex in vortex-free rotation. For comparison, the calculated veloci- 
ties of the two vortex ends in Fig. 6 are v^ z ~ 0.84 aQR ~ 0.96 av(R) and 
v L(j) w 0.73(1 - a')£lR w 0.83(1 - a')v(R). The wave length of the spiral tra- 
jectory is A = 2irRvi jZ /v\ J( j ) ~ 5 mm and the period p = 2nR/v^ ~ 50s. 
Thus the end point motion can be used to construct a simplified model of the 
motion of the two short vortices in Fig. 6. 

As seen in Fig. 6, the spiral motion of the vortex end point along the 
cylindrical wall winds the rest of the evolving vortex around the central vortex 
cluster with a wave vector Q, such that QR = v^/v^z = (. The other almost 
straight end of the evolving vortex is fixed to a flat end plate of the cylinder and 
resides there at the edge of the vortex cluster, where the counterflow velocity 
is close to zero. Therefore the helical twist is removed only by a slow sliding of 
the vortex end along the end plate. As seen in Fig. 6, occasional reconnections 
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between the twisted evolving vortex and a straight outer vortex in the cluster 
or with a second oppositely twisted vortex can help to reduce the twist. Finally, 
we see in Fig. 6 that while the evolving vortex is wound tightly around the 
cluster this induces Kelvin wave oscillations which propagate vertically along 



the vortices in the cluster (Hanninen et al. 2007a) 



2.3. Onset temperature of turbulence 

When the dynamics calculated in Figs. 5 and 6 is probed with measure- 
ments, the final state proves to depend crucially on temperature: At tem- 
peratures above the transition to turbulence, T > T^ 114 , the calculations are 
confirmed and the number of vortices remains constant. At low temperatures, 
in contrast, evolving vortices may become unstable in applied counterflow and 
trigger a turbulent burst. The evolution after the burst continues as illustrated 
in Fig. 3. The final state is then consistently close to the equilibrium vortex 
state. Interestingly, it turns out that for each initial configuration, such as 
those in Figs. 5 or 6, there exists a specific temperature T on , which character- 
izes the onset of turbulence: well above T on no turbulent burst occurs while 
well below T on a burst always occurs. 

In these onset measurements only two different types of final states are ob- 
served: sufficiently far above T on the number of vortices is conserved, while well 
below T on close to the equilibrium number of vortices is formed. The change 
over takes place within a narrow temperature interval around T on , typically 
within ±0.02 T c . Within this interval either of the two final states can emerge. 
The reason for the narrow width is the strong, nearly exponential temperature 
dependence of the mutual-friction parameter £ = (1 — a') /a, which controls 
the dynamic instability of seed vortices evolving in the applied counterflow 
(Fig. 1). As sketched in Fig. 3, two sequential processes are needed to start 



the low-temperature evolution: First the single- vortex instability (Finne et al 



2006a), the precursory process which becomes possible only at temperatures 
below T££ and which is responsible for generating a bunch of new evolv- 
ing vortices. Secondly, the turbulent burst has to be triggered as a collective 
process in which several evolving vortices interact and generate a sudden lo- 
calized event of turbulence, which expands across the entire cross section of 
the rotating column, but only over a short section of its length (of order ~ R). 

The experimental confirmation of this scenario is obtained by examining 
the final state as a function of temperature. By recording the line shape of 
the NMR absorption profile, when the magnitude of the magnetic polarization 
field is swept across the resonance region, the number of vortices in the final 
state can be determined. The line shapes of the two types of final states differ 
drastically, as seen in Fig. 7, where all the spectra have been recorded at the 
same temperature of 0.40 T c and where thus the integrated area under each 
line shape is the same. The characteristic feature are the large NMR shifts. 
These are controlled by the temperature and pressure dependent spin-orbit 
coupling. If the central vortex cluster is surrounded by applied counterflow at 
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Fig. 7. On the left the NMR absorption spectra constitute an image of the "flare-out r 



order parameter texture in the long rotating cylinder (Kopu et al. 2000 1 . The Larmor 
field, around which the NMR absorption is centered in the normal phase, is here 
at 48.1 mT. In the B phase the Larmor value becomes a sharp edge beyond which 
at higher fields the absorption vanishes, while most of the absorption is shifted 
to lower fields. The dominant absorption maximum on the left is the counterflow 
peak. Its height depends on the number of vortices N in the central cluster. N 
can be defined in terms of the rotation velocity Q v at which the vortices form the 
equilibrium state: N = N cq (£l v ). Here all counterflow peaks have been recorded at 
the counterflow velocity O, — £l v = 0.8 rad/s. The conversion from peak height to 7Y 
in given conditions (T, O, P) can be obtained experimentally or from a numerical 



calculation of the order parameter texture and the corresponding line shape ( Kopu 



2006|. The two line shapes without a counterflow peak, but with a large truncated 
maximum bordering to the Larmor edge, are for the non-rotating sample (O = 0) 
and for the equilibrium vortex state (N = N eq ) at £l v = 0.9 rad/s. For both of them 
the absorption at the site of the counterflow peak is close to zero. On the right a 
NMR setup is shown which was used to study the onset temperature of turbulence 
(cf. also Fig. 13). Two different contours of the AB interface are shown when the 
A-phase barrier layer is present at 0.6 T c . The contours correspond to a current of 
4 A (narrow A-phase sliver with curved concave interfaces) and 8 A (wider A-phase 
layer with flat interfaces) in the superconducting A-phase stabilization solenoid 



(Finne et al. 2004b). 



some sizeable velocity, a large sharp peak is formed which is shifted down field 
from the Larmor value. The number of vortices in the central cluster can be 
determined from the height of this so-called counterflow peak. At small vortex 
numbers (N N eq ), the reduction in the counterflow peak height is directly 
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Fig. 8. Number of remanent vortices JV( At) for A-phase separated top and bottom 
sample sections, measured as a function of the annihilation time At in the temper- 
ature regime of laminar vortex motion. The results of these two independent mea- 
surements can be fitted in both cases with the solid curve jV(At) ~ 2- 10 3 / (a At + 6) 
[At in sec, b ~ 7 s]. (Inset) In the upper right corner the sequence of rotations f2(t) 
is shown which was used to perform the measurement. The data are for 0.57 T c with 
a = 0.60 and 0.7T C with a = 1.1. Parameters: Of = 0.9rad/s, R = 3mm, length of 
top sample section ht = 44 mm (41mm) and = 54 mm (51mm) for the bottom 
section at 0.57 T c (0.70 T c ). 



proportional to the number of rectilinear vortices N in the cluster. At larger 
N the dependence becomes nonlinear and ultimately the peak height drops 



to zero well before N eq is reached, in practice around v(R) < lmm/s (Kopu 



2006). Accordingly, in the equilibrium vortex state the line shape is radically 



different and easily distinguished from a state with sizeable counterflow. 

A measurement of T on for any particular initial setup, such as in Fig. 5 for 
remanent vortices, requires repeating the measurement at different tempera- 
tures and recording the line shape in the final state. Surprisingly, it turns out 
that in the final state the vortex number is either preserved or it has increased 
close to that in the equilibrium vortex state. Practically no intermediate cases 
are observed. As a result, in practice a measurement of T on requires simply a 
visual check of the measured line shape in the final state. This feature about 
the turbulent burst is similar to recent observations from measurements with 
a closely spaced pair of vibrating wires in superfluid 4 He (iGoto et a/.||2007ft. 
One of the wires is driven at high oscillation amplitude as generator while the 
second is operated at low amplitude as detector. Once turbulence has been 
switched on by running the generator at high drive, the generator can be 
switched off and turbulent flow will still be maintained around the detector. 
Only if the detector drive is reduced to sufficiently low level, turbulence ceases 
and the flow around the detector returns to the laminar state. This shows that 
once turbulence has been switched on it can be sustained at much lower flow 
velocities. Similarly, once the turbulent burst is started in the rotating col- 
umn, turbulent vortex formation will continue until the counterflow velocity 
has dropped close to zero in a section of the column of height ~ R. 
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The situation at temperatures above T on is illustrated by the measurements 
on vortex remanence in Fig. 8 ( Solntsev et a/.|2007 ). These measurements have 
been performed at two different temperatures above onset, where no increase 
in the number of vortices is expected. This is confirmed by extracting from 
the counterflow peak height the number of vortex lines in the final state. The 
measurement proceeds as sketched in the inset on the top right of Fig. 8: An 
equilibrium vortex state is decelerated to zero and the vortices are allowed 
to annihilate for a time interval At before rotation is turned back on. The 
measurement is repeated many times by varying the annihilation time At 
at zero rotation. The annihilation time is found to govern the number of 
rectilinear vortex lines in the final state and thus the number of remnants at 
the end of the annihilation period: N(At) oc (l + A£/r m f) _1 , where the mutual- 
friction-controlled time constant is r m f = [2afi ] _1 an d — < 0). This 
is exactly as expected for the mutual-friction damped motion of vortices in 
the radial direction, when straight vortices move outward to annihilate on the 
cylindrical wall at zero rotation. For this to apply, the vortices have to be 
polarized along the cylindrical symmetry axis dKrusius et oT]l993^ As seen m 
Fig. 5 (second from left, at t < 600 s) this is the case: The polarization remains 
at high level even in the remanent state at zero rotation. Consequently, the 
measurements in Fig. 8 confirm that at constant temperature above T on the 
number of vortices in this experiment is controlled by the annihilation period 
At and no uncontrolled increase occurs. 

It is useful to note some additional features about vortex remanence in 
the measurements of Fig. 8. Let us denote the number of remnants after the 
annihilation period with Af(At) and the number of rectilinear lines in the 
final state with N. Although the annihilation time At controls the number 
of remnants Af(At), the result N = Af is independent of At. It is also to a 
large extent independent of how the measurement is performed, i.e. what the 
rotation velocity Qf in the final state is or what acceleration Cl is employed 
to reach Qf (as long as Af(At) < N cq (Qf) or the critical velocity for vortex 
formation, Qf — Q v < v c /R, is not exceeded). Furthermore, the result N = Af 
is established separately both for the top and bottom sections of the cylinder, 
when these are separated by a magnetic-field stabilized A-phase barrier layer 
(cf. Fig. 7), and for the entire cylinder without A-phase barrier. The A-phase 
barrier layer prevents vortices from traversing across the AB interfaces at low 
counterflow velocity ( |Blaauwgeers et al. 2002). In this way the number of 
remanent vortices Af(At) has been found to be proportional to the length h of 
the cylinder (as long as h ^> R). All these properties are consistent with the 
conclusions that when T > T on , the vortex number is conserved in dynamical 
processes, the annihilation decay of remnants is a laminar process regulated by 
mutual friction damping, and that pinning is weak. For simplicity, we neglect 
vortex pinning altogether and assume ideal wall properties throughout this 
review. 

The situation at temperatures around T on is illustrated by the measure- 
ments in Fig. 9 which determine T on for this particular choice of initial state 
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Fig. 9. Measurements on the onset temperature T on of the transition to turbulence. 
The measurements are performed similar to those in Fig. 8 and start from an initial 
state which is obtained by decelerating an equilibrium vortex state at 1.7rad/s 
to zero at a rate 0.01rad/s 2 . The remaining vortices are left to annihilate for a 
period At at f2 = 0. Rotation is then increased to Of at a rate 0.02rad/s 2 . When 
all transients have decayed the number of vortices is measured in the final steady 
state at Of. The result is plotted as a function of temperature with 30 - 40 data 
points per panel. The solid curve is a gaussian fit which represents the probability 
for turbulence with a half width <tt = 0.02 T c , centered around T on . Comparing 
results in the two panels for At = 20 min and 2 min, we see that T on decreases with 
increasing At, since the number, average size, and density of remnants is reduced 
as At increases. Both panels have been measured for the upper sample section 
which is separated from the bottom half with an A-phase barrier layer. Parameters: 
R = 3mm, h = 45 mm, and P = 29.0 bar. 



(Solntsev et al. 2007, de Graaf et al. 2007). The probability of the turbu- 
lent burst is plotted as a function of temperature, when the annihilation time 
At = 20 min (on the left) and At = 2 min (on the right). The striking feature 
is the abrupt change over from the laminar behaviour, where the vortex num- 
ber is conserved, to turbulence, where the vortex number surges close to iV eq 
and the system relaxes to its minimum energy state. The center of the narrow 
transition defines the onset temperature of turbulence T on , which proves to be 
different for the two cases studied in Fig. 9. 

As seen in Fig. 9, T on depends on the annihilation time At and thus on the 
initial number and configuration of evolving remnants A/"(At) at the moment 
when the step increase in rotation from zero to Of is applied. Calculating 
from the results in Fig. 8 one finds that the number of remnants at the start 
of acceleration is AA ~ 40, when At = 2 min and T m T on = 0.44 T c , and 
A/j ~ 10, when At = 20 min and T = 0.39 T c . Thus at a higher temperature 
a larger number of remnants is needed to achieve the turbulent burst. Both 
panels in Fig. 9 refer to the top sample half where, with no orifice, there 
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is no preferred site for the remnants to accumulate and the turbulent burst 
occurs randomly at any height z in the column ( de Graaf et a/.|[2007 ). Similar 
measurements at different values of At and final rotation velocity Qf show that 
the onset temperature depends weakly on both the initial number of remnants 
A/j and the applied flow velocity. These dependences can be summarized in 
the form 

CCTon)- 1 ^^ 3 ^ 1 - 3 . (5) 

Thus the onset temperature T on depends primarily on the mutual friction 
parameter ((T), but also weakly on other factors which influence the likelihood 
of achieving locally somewhere in the maximum available counterflow velocity 
a density of evolving vortices which allows to trigger the turbulent burst. 
Among these additional factors most important are (i) the applied counterflow 
velocity v — v n — v s , (ii) the number and configuration of the injected seed 
vortices, and (iii) the sample geometry. In Figs. 8 and 9 we examined the 
response of remanent vortices to a step-like increase in rotation. The same 
measurements can also be performed by starting from the equilibrium vortex 
state at finite rotation, as discussed in the context of Fig. 6. In fact, the 
most extensive study of the scaling law in Eq. (5) was performed using this 
approach. 

Finally we note that in the rotating container all measured onset temper- 
atures, which depend on the presence of the precursor, are found to be below 
the transition to turbulence in the bulk: T on < T^ lk . Furthermore, since the 
onset also depends on the applied counterflow velocity in Eq. (5), the insta- 
bility is expected to occur first close to the cylinder wall, where the applied 
velocity v = Qr reaches its maximum value at r = R. Thus the reconnection 
of the expanding loop will most likely occur with the wall. Surprisingly it is 
also found that once the instability is triggered, the turbulent burst essen- 
tially always follows next, since little if no increase in the vortex number is 
detected at T ~ T on in such cases where the turbulent burst does not switch 
on (cf. Fig. 9). To provide more understanding on the role of the single- vortex 
instability as the precursor mechanism to turbulence, we next examine it in 
the onset temperature regime, T ~ T on , where the instability proceeds suffi- 
ciently slowly in time so that it can be monitored with continuous- wave NMR 
measurement. 



2.4- Single-vortex instability in applied flow 



Since the time when it was first understood that superfluid turbulence is 



made up of tangled quantized vortices (Hall and Vinen 1956, Vinen 1961) 



the most basic question has been its onset as a function of applied counter- 
flow velocity: How is turbulence started and what defines its critical velocity? 



An important clue was provided by the rotating experiments of Cheng et al. 



(1973) and Swanson et al. (1983), who found that rectilinear vortex lines in 



rotation are broken up in a turbulent tangle if a heat current is applied par- 
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allel to the rotation axis. The thermal current is transported as a counterflow 
of the normal and superfluid components along the rotation axis. Rectilinear 
vortices become unstable in this parallel flow and above a low critical velocity 
transform to a tangle which tends to be aligned in the plane transverse to the 
heat current. 



This phenomenon was explained by Glaberson and his coworkers (Glaber 



son et al. 1974 ) who showed that an array of rectilinear vortices becomes un- 
stable in longitudinal counterflow above the critical velocity v = 2(2VLk c ) 1 I 2 , 
where K e « k is an effective circulation quantum (supplemented with the log- 
arithmic cutoff term: K e = (k/Att) ln(£/ao), where the average inter- vortex 
distance is I ~ (k/2VL) 1 / 2 and the vortex core radius do). The instability ap- 
pears when a Kelvin-wave mode with wave vector k = (2f2/K e ) 1//2 starts to 
build up, whose amplitude then grows exponentially in time. The Glaberson 
instability has also been examined in numerical calculations which qualita- 
tively confirm the instability and the vortex tangle in the transverse plane 



which starts to form above a first critical velocity (Tsubota et al. 2003) 



In general, the dispersion relation of a helical Kelvin wave disturbance cx 



exp[-i(o;ki — kz)] can be written as (Donnelly 1991, Finne et al. 2006b) 



^k(k) = n e k 2 — a'(n e k 2 — kv) — ia(n e k 2 — kv) . 



(6) 



In the absence of flow (v = 0), these Kelvin modes are always damped, at 
high temperatures they are actually overdamped, but at low temperatures 
(a < 1 — a') this is not the case. In applied flow (v > 0), the long wave 
length modes with < k < v/n e become exponentially unstable. If an evolv- 
ing vortex accumulates enough length L» parallel to the applied flow, then a 
disturbance with wave length A m i n ~ L» ~ K c jv ~ 1/A; max may start to grow. 
The expanding loop may reconnect, either with the wall of the container, with 
itself, or with another vortex. This leads to a growing number and density of 
evolving vortices, which ultimately start interacting and trigger the onset of 
turbulence in the bulk. 

No rigorous analytical calculation has been presented of the single-vortex 
instability in the rotating container, but a simple scaling model illustrates the 
problem. Consider a vortex ring in vortex-free counterflow, which is initially 
perpendicular to the plane of the ring in a rotating cylinder of radius R. If the 
ring is large enough, then it expands until it reaches the container size R. The 
time needed for this expansion is of order 5t ~ R/av, where v is the average 
normal velocity through the ring. The ring also has a self-induced velocity 
component v T ~ K e /R, which arises from its own curvature and is directed 
along the normal of the plane of the ring. Because of this velocity component, 
the plane of the ring is rotated away from being perpendicular to the azimuthal 
flow in the cylinder, while it drifts in the flow. During the time 5t, the vortex 
length parallel to the flow becomes of order (1 — a')v T 8t. Equating this to 
it is seen that the instability condition Ly > A min leads to the requirement 
( = (1— a')/a > 1. This condition is virtually independent of velocity; the only 
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restriction is imposed by the finite container radius, Lii < R, which defines 
a critical velocity v c ~ n e /R. Typically the time spent by a vortex in radial 
motion before reaching the sample boundary is of order R/(aQR) = (aQ)^ 1 . 
At 0.45 T c and Q = 0.6rad/s, the inverse of this quantity equals 0.20 s -1 , which 
fits with the measured vortex generation rate of dN/dt = N = 0.23 s" 1 in 
Fig. 10. However, numerical calculations confirm that the presence of surfaces 
is required to demonstrate the single- vortex instability in usual experimentally 
relevant flow conditions (Finne et al. 2006a). Also the calculations demonstrate 
that the instability is not characterized by a unique critical velocity, since it 
depends on the relative orientation of the flow with respect to the vortex, 
while the vortex expands in helical motion in the rotating cylinder. Thus the 
above model is incomplete. 

Let us now examine direct observations of the single-vortex instability in 
rotating flow. In the onset temperature regime, T ~ T on , in about half of the 
measured cases, which lead to a turbulent burst, a slow increase in the number 
of vortex lines N(t) can be observed to precede the turbulent burst. If present, 
the increase is invariably followed by a turbulent burst. Thus it appears rea- 
sonable to associate the slow increase with the single-vortex instability. In 
Fig. 10 the number of vortex lines N(t) is plotted as a function of time while 
the precursor generates new vortices at slow rate. In this example the increase 
in N(t) is almost linear (solid line), until the turbulent burst sets in and starts 
the vortex front motion along the rotating column both upwards and down- 
wards from the site of turbulence. At flf = 0.6rad/s the slow increase lasts in 
this example for about 140 s, generating approximately one vortex every five 
seconds, until some 30 new vortices have been created and the turbulent burst 
manages to switch on. The time interval from t — to the turbulent burst 
is called the burst time which here is % = 140 s. At larger Of the burst time 
is shorter in duration, eg. in a repetition of the measurements in Fig. 10 at 
fit = lrad/s the turbulent burst was found to start in less than 30 s. 

Two further observations about the precursor can be made from Fig. 10. 
First, vortex formation proceeds independently in different parts of the sample. 
At flf = 0.6 rad/s it takes more than 300 s for a vortex created at one end of the 
sample to reach the other end. Still, vortex formation at the top and bottom 
ends is observed to proceed at roughly the same rate. Thus vortex generation 
by the single-vortex instability is not localized, in contrast to the turbulent 
burst. The random occurrence of the single- vortex instability agrees with the 
notion of ideal walls (or at least weak pinning), as opposed to a vortex mill 
localized at a surface defect on the cylinder wall. 

Secondly, in Fig. 10 the equilibrium vortex state at low initial rotation 
Qi = 0.05 rad/s has been used to introduce evolving vortices in the applied 
flow. This approach provides a more reproducible initial vortex configuration 
than remanent vortices, since the number of those vortices, which connect to 
the cylindrical side wall, is primarily determined by the misalignment between 
the cylinder and rotation axes (Fig. 6). In a given experiment the residual angle 
between the two axes is generally a constant. 
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Fig. 10. Experimental illustration of single-vortex instability as precursor of bulk 
turbulence. The number of vortices N(t) is recorded with NMR coils at the top 
and bottom ends of the sample. As seen in the inset, initially the sample is in the 
equilibrium vortex state at Oi = 0.05 rad/s with N ~ 30 vortices, of which close to 
one half connect to the cylindrical side wall. Rotation is then increased to a new 
stable value Of = 0.6 rad/s, which is reached at t = 0. During the ramp to Of the 
counterfiow builds up, compresses the rectilinear sections of all vortices to a central 
cluster, and starts the spiral motion of the vortex ends connecting to the side wall. 
Eventually in the increased applied flow at Of the instability starts to generate new 
vortices which contribute to the average slow rate N of vortex formation, shown by 
the solid straight line. After about 140 s the turbulent burst occurs 63 mm above 
the bottom end plate. Vortex fronts traveling up and down along the column then 
approach the two detector coils and reach their closer ends as indicated by vertical 
arrows (at 230 s and 275 s). The filled data symbols are derived from experimentally 
calibrated counterfiow peak heights and the open symbols from order parameter 
texture calculations fitted to the NMR signal in the non-rotating state. 



To appreciate the influence of the vortices curving to the side wall, the 
experiment was repeated differently. A cluster with only rectilinear vortices 
(N < N cq ) was prepared at higher temperatures and was then cooled below 
0.5 T c . As long as this cluster is separated by a sufficiently wide vortex-free 
counterfiow annulus from the cylindrical boundary, Q can be increased or 
decreased without change in N at any temperature down to 0.35 T c (which 
is the lower limit of the so far measured onset temperatures T on ). If Q is 
reduced too much, the cluster makes contact with the cylindrical side wall, 
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some outermost vortices become curved, and during a subsequent increase of 
Q, while T < 0.5 T c , the behavior in Fig. 10 is reproduced. Therefore we are 
led to assume that, to observe the vortex instability, at least one curved vortex 
connecting to the cylindrical side wall needs to be present. At temperatures 
below 0.35 T c this may not be the case, since in rapid changes of rotation even 
rectilinear vortices seem to be destabilized (Fig. 4). 

More statistics on the properties of the precursor has been collected from 



measurements similar to that in Fig. 10 by de Graaf et al. (2007). Important 
characteristics are the initial rate of vortex generation N(t = 0) and the burst 
time tb- These can be examined for events with sufficiently long burst times 
tb > 20 s, so that the rate of the counterflow peak height decrease with time 
can be adequately resolved. In general it is found that N increases and % 
decreases rapidly with decreasing temperature below T on . To find events with 
well-resolved pre-turbulent vortex generation and long burst time one has to 
scan for data (i) in the onset temperature regime, T « T on , (ii) with a low 
initial formation rate N < lvortex/s, and (iii) by starting from a state with 
a small number of seed vortices (iv) at low applied flow velocity. 

These measurements demonstrate that the precursor generates new inde- 
pendent vortex loops which start to evolve along spiral trajectories towards 
the final state of a rectilinear vortex line. When the density of evolving vortices 
rises sufficiently, so that interactions between them in the bulk volume become 
possible, then the process is terminated in a turbulent burst. The burst is a 
localized event which from one measurement to the next happens randomly 
at different heights z of the sample (Fig. 7). The measured properties of the 
precursor are consistent with those expected for a single-vortex instability 
based on the excitation of Kelvin-wave modes of sufficiently long wave length. 
Overall, measurements in the onset regime reveal the precursor mechanism, 
owing to the strongly temperature dependent mutual friction of 3 He-B, which 
makes the precursor observable within a narrow temperature interval around 
the onset temperature. At lower temperatures the turbulent burst develops so 
rapidly that the measuring techniques, which have been employed so far, are 
not fast enough to capture the details. The latter case is the typical situation 
in superfluid 4 He experiments. 



2. 5. Numerical calculation of dynamic vortex generation 

Numerical calculations on vortex dynamics are carried out with the vortex 
filament model introduced by Schwarz (1988). With today's computing power, 
one uses Biot-Savart integration along all vortex lines, so that the superfluid 
velocity field from vortices is obtained from (Hanninen et al. 2005) 



(7) 
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t= 50.5 s f = 77.5s 



Fig. 11. Two snapshots from a calculation of vortex generation in a rotating cylinder. 
The summary of these calculations with results accumulated over more than 100 s 
is shown in Fig. 12. 



The line integral is taken along all vortices in the system, s(£, t) denotes the 
location of the vortex core at time t, and £ is measured along the arc length 
of the vortex core. In the presence of solid boundaries the total superfluid 
velocity field, v s = v Sj£J + v b , is modified by the boundary induced velocity 
Vb- At a plane boundary one can use image vortices to satisfy the requirement 
of zero flow through the boundary, n ■ v s = 0, where fi is the unit vector 
along the surface normal. More generally we obtain v b = V$ by solving 
the Laplace equation V 2 $ = combined with the requirement that at the 
boundary h ■ V$ = — h ■ v S)£J . 

No surface pinning or even surface friction is generally included, the bound- 
aries are assumed ideal, as indicated so far by measurements on 3 He-B in 
smooth-walled simple cylindrical containers. Mutual friction in the bulk su- 
perfluid is included using the equation of motion (4) for the vortex element 
at s(£, i), which moves with the velocity vl = ds/dt. For the mutual friction 



parameters a(T, P) and az'(T, P) one uses the 3 He-B data measured by Bevan 



et al. (1997) at 10 and 29 bar pressures. A reconnection between two vortex 



segments is enforced if they have drifted within a distance from each other 
which is less than the minimum spatial resolution of the calculation (usually 
~ 0.05 mm). The configuration after reconnection should correspond to shorter 
overall vortex length than the initial state. In practice, the computing time 
limits severely what can be calculated and what becomes too time consuming. 
Therefore the practical implementation becomes of great importance, how the 
Biot-Savart integration and the proper solution for the boundary conditions 
are worked out. For details we refer to de Graaf et al. (2007). 

In Fig. 11 two snapshots are shown from calculations on vortex formation 



and the configurations which evolve in a rotating cylinder (de Graaf et al. 
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2007). Recently formed younger vortices are here in helical configurations on 



the outer circumference closer to the cylindrical wall. There in the outer re- 
gions one can see loops of Kelvin waves, small separated loops with both ends 
of the vortex on the cylindrical wall, and even closed vortex rings (lower right 
corner at t — 50 s). Since it is primarily surface reconnections at the cylindrical 
wall which contribute to the formation of new vortices in the early stages of 
the calculation (at t < 100 s), the many newly formed short loops are still close 
to the side wall. Further inside the cluster one can see older and straighter 
vortices which congregate within the central parts. 

The general observation from these calculations is that evolving vortices 
in a rotating sample are more stable in the numerical experiment than in 
measurements. For instance, in Fig. 11 vortex formation has to be started 
from an artificial initial configuration (Finne et al. 2006a). This consists from 
an initial single vortex ring which is placed in the plane perpendicular to the 
rotation axis at height 0.2 h slightly off center, to break cylindrical symmetry 
Qde Graaf etaT^mf) . This is an unstable configuration where Kelvin waves of 
large amplitude immediately form and then reconnect at the cylindrical wall. 
The end result is the sudden formation of roughly 30 vortices which have one 
end on the bottom end plate and the other moving in spiral trajectory along 
the cylindrical wall. After the initial burst the later evolution is followed as a 
function of time t, the number of vortices N(t) is listed, and the reconnections 
of different type are classified. The results are shown in Fig. 12. 

In Fig. 12 one keeps account of all reconnection processes which occur in 
the rotating sample as a function of time while it is evolving towards its fi- 
nal stable state with an array of rectilinear vortices and iV — * N ccl . After the 
initial burst of the first ~ 30 vortices N increases first gradually, but after 
about 50 s the rate N picks up. During the first 50 s reconnections in the bulk 
do not contribute to the generation of new vortices, but later such processes 
also start to appear. However, even during the later phase a reconnection of 
a single vortex at the cylindrical wall, while Kelvin waves expand along this 
vortex, remains the dominant mechanism of vortex generation. This is seen 
from the fact that the curve for N follows closely that of the successful surface 
reconnections (dashed curve marked as "AiV = +1"). The most frequent re- 
connections after the first 40 s are denoted by the solid "AiV = 0" curve and 
occur in the bulk between two different vortices. These inter-vortex reconnec- 
tions do not lead to changes in N and are primarily associated with processes 
occurring between the twisted vortices in the bundle further away from the 
wall. 

The inset in Fig. 12 compares the rates of vortex generation from reconnec- 
tions at the wall and in the bulk: The reconnection of a single vortex at the 
cylindrical wall is clearly the most important mechanism for the generation 
of new independent vortex loops in the early stages of the calculation. The 
dominant role of such wall reconnections is compelling. A second important 
consideration is correspondence with measurement. The obvious difference is 
the higher stability of evolving vortices in the calculation as compared to ex- 
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Fig. 12. Calculation of the cumulative number of reconnections and vortices in a 
rotating cylinder. The different curves denote: (AN = 0, solid curve) reconnections 
in the bulk which do not change N; (+1, dashed) reconnections with the cylindrical 
wall which add one new vortex loop; (N, solid) total number of vortices; (removed, 
dash-dotted) small loops which form in reconnections mainly close to the cylindrical 
wall, but which are contracting and are therefore removed; (+1, solid) reconnections 
in the bulk which add one vortex and (-1, solid) which remove one vortex; (0, dashed) 
reconnections at the cylindrical wall which do not change N. (Inset) Averaged rate 
of increase in N owing to reconnections on the cylindrical wall and in the bulk. The 
large initial peak in the boundary rate represents the starting burst, which is used 
to start vortex formation. Parameters: R = 3 mm, h = 10 mm, $7 = 0.9rad/s, and 
T = 0.35 T c (where a = 0.095 and a' = 0.082). 



periment. In Fig. 12 the rate of vortex generation remains modest, no clearly 
identifiable turbulent burst can be distinguished, and the vortex number ap- 
proaches the equilibrium value from below. After 115 s of evolution the vortex 
number has progressed to A" w 400 , where the increase is almost stopped, well 
below the saturation value of N eq « 780 ( de Graaf et a/.|2007 ). This is the gen- 
eral experience from calculations on an ideal rotating cylinder, with smooth 
surfaces and no surface friction or pinning. The calculations become more and 
more time consuming with decreasing temperature, which limits the possibil- 
ities to obtain a more comprehensive understanding of their predictions and 
of the origin of the differences with measurement. The low probability of the 
single-vortex instability in the calculations appears to be a particular prop- 
erty of rotating flow in a circular cylinder, since linear pipe flow, for instance, 
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displays a steady rate of vortex generation (de Graaf et al. 2007). 

Clearly numerical calculations provide important illustrations and guidance 
in situations where measurements answer only specific limited questions. The 
calculations take full account of interactions between vortices and between a 
vortex and the ideal container wall. Nevertheless, the correspondence between 
calculation and measurement is not satisfactory at present, when we speak 
about the single-vortex instability and the onset of turbulence in a rotating 
cylinder. It appears that some mechanism, which makes vortices more unstable 
and adds to the vortex generation rate, is missing from the calculations. The 
difficulty is likely to reside on the cylindrical wall, where the assumption of 
ideal conditions should be examined closer. Attempts in this direction have 
so far not produced more clarification. However, these uncertainties about the 
mechanisms behind the single- vortex instability in rotating flow do not change 
the fact that at low vortex density Kelvin-wave formation on a single vortex, 
followed by a reconnection at the surface, is the only efficient mechanism for 
generating new vortices. 



2.6. Summary: onset of turbulence 



Since the advent of 3 He-B new possibilities have appeared to study tur- 
bulence. Firstly, it has become possible to distinguish and characterize, in 
measurements with large samples, vortex formation at a stable reproducible 
critical velocity, vortex remanence and turbulent proliferation of vortices. Sec- 
ondly, the mutual friction dissipation a(T) with strong temperature depen- 
dence around a ~ 1 has made it possible to evaluate the role of mutual friction 
in the onset of turbulence. The important dynamic parameter proves to be 
£ = (1 — a') I a. It controls the onset of the single-vortex instability, where 
an evolving vortex becomes unstable and generates, during a reconnection at 
the wall, a new vortex loop. After several such events the density of evolving 
vortices is sufficient to produce a turbulent burst. The necessary condition 
is ( > 1, to start the cascade of the single- vortex instability followed by the 
turbulent burst. 

The single-vortex instability becomes possible only at temperatures below 
the turbulent transition in the bulk volume and thus T on < . The onset 
temperature T on of these two series-coupled processes has been found to obey 
a power-law dependence which relates the mutual friction parameter ( to the 
magnitude of the "flow perturbations" in Eq. (5). Well above T on no new 
vortices are detected (with a resolution < 10 new vortices), while well below 
T on all final states are found to be equilibrium vortex states with close to the 
equilibrium number of vortices, N < iV eq . In the onset regime itself, T ~ T on , 
one finds events with and without turbulent burst, but surprisingly practically 
no incomplete transitions with Ni < N <C N eq . 

In the intermediate temperature regime 0.3 T c < T < 0.6 T c , the equilib- 
rium vortex state is reached after a single turbulent burst. In fact, in the 
measurements with the sample setup of Fig. 7 no case of two or more almost 
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simultaneous bursts was identified above 0.35 T c . Apparently the probability 
of the single vortex instability to start a turbulent burst is still low at these 
temperatures. Secondly, after the burst the vortex front moves rapidly and 
removes the vortex-free flow. At these intermediate temperatures the burst is 
both spatially and temporally a localized event in a short section (of length 
~ R) of the column. From one measurement to the next, it occurs randomly 
at different heights of the column. Below 0.3 T c the longitudinal propagation 
velocity of vortices becomes slow and evolving vortices go rapidly unstable 
everywhere. As a result turbulence tends to be both spatially and temporally 
more extended, filling larger sections of the column. The later events, the 
evolution after the turbulent burst, are the subject of the next section. 



3. Propagating vortex front in rotating flow 

3. 1 . Introduction 

In rotation at constant angular velocity the steady state superfluid response 
is generally not turbulent. Nevertheless, transient states of turbulence can be 
formed by rapidly changing the rotation velocity, especially if the sample con- 
tainer does not have circular cross section or its symmetry axis is inclined by a 
larger angle from the rotation axis. The decay of turbulence and the approach 
to equilibrium can then be monitored at constant Q. The normal component 
relaxes back to solid body rotation by means of viscous interactions, while the 
superfluid component adjusts much slower, coupled only by mutual friction 
dissipation from vortex motion with respect to the normal component and (if 
any) by the deviations of the container walls from being axially-symmetric 
around the rotation axis. Such measurements on transient turbulence are gen- 
erally known as spin up or spin down of the superfluid component. This used 
to be an important topic in superfluid 4 He work in the fifties and sixties 



(Andronikashvili and Mamaladze 1967), but was then replaced (with few ex 



ceptions (Adams et al. 1985)) by other methods which one expected to lead 
to results with more straightforward interpretation. 

The turbulent burst, which suddenly starts the motion of N pa N m vortices 
along the rotating column at temperatures T < T on , as discussed in Sec. 2 (c/. 
Fig. 3), provides a novel technique to investigate transient turbulence in rota- 
tion. Originally it was assumed that this motion would take place as a tangle 
of vortices, which spreads longitudinally along the rotating column. It was 



soon realized from NMR measurements (Eltsov et al. 2006b) that this could 
not be the case; rather the propagating vortices were highly polarized and had 
to be coiled in a helical configuration owing to their spirally winding motion. 
This recognition presented a new problem: Is there any room at all for turbu- 
lence in this kind of motion and if there is, how is it expressed? Or perhaps 
the nature of the motion changes on approaching the zero temperature limit, 
when mutual friction dissipation vanishes a oc exp(— A/T)? These questions 
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provided the incentive to examine the propagation more closely and to mea- 
sure its velocity as a function of temperature. The results demonstrate that 
turbulent losses depend crucially on the type of flow, flow geometry, external 
conditions, the physical properties of the superfluid itself, etc. 

A measurement of the front propagation in the laminar and turbulent tem- 
perature regimes allows one to determine the rate of kinetic energy dissipation. 
The measurement proceeds as follows: The initial starting state is the rotating 
vortex-free state, the so-called Landau state, which is metastable with much 
larger free energy than the stable equilibrium vortex state. The latter consists 
of rigidly co-rotating normal and superfluid components, owing to the pres- 
ence of a regular array of rectilinear vortices, while in the vortex-free state the 
superfluid component is not rotating at all: it is at rest in the laboratory frame 
of reference. When the turbulent burst is triggered in the Landau state, a rapid 
evolution towards the equilibrium vortex state is started, where a boundary 
between the vortex-free and the vortex states propagates along the rotating 
column and displaces the metastable vortex-free counterflow. Particularly at 
temperatures below 0.4 T c the boundary has the form of a sharp thin vortex 
front which travels at a steady velocity Vf. The dissipation rate of the total 
kinetic energy, £(t), is related to Vf as 

dS jdt = -n Ps V{ n 2 R 4 /4 . (8) 

By measuring Vf one determines directly the energy dissipation dS/dt as a 
function of temperature. 

At high temperatures the motion is laminar and the front velocity is de- 
termined by mutual friction dissipation between the normal and superfluid 
components, Vf(T) ~ a(T)QR. Below 0.4 T c , Vf(T) deviates more and more 
above the laminar extrapolation (Eltsov et al. 2007[ ), in other words the dissi- 



pation becomes larger than expected from mutual friction in a laminar flow. At 
the very lowest temperatures a striking anomaly becomes apparent: d£(t) j dt 
does not go to zero, but the measured velocity Vf(T) appears to level off at 
a constant value which corresponds to an effective friction « c ff ~ 0.1, even 
though a{T) — > 0, when T — > 0. Evidence for a similar conclusion has been of- 
fered by the Lancaster group ( Bradley et a/.|[20 06), who measured the density 



of the vortex tangle created by an oscillating grid and found that this kind of 
turbulence decays at a temperature-independent finite rate below 0.2 T c . 

When mutual friction decreases and turbulent motions in the vortex front 
cascade downward to progressively smaller length scales, eventually individual 
quantized vortex lines must become important. This is the quantum regime of 
superfluid hydrodynamics. The energy cascade on length scales smaller than 
the inter-vortex distance and the nature of dissipation on these scales are 
currently central questions in superfluid turbulence ( Vinen and Niemela|2002 ). 



Theoretical predictions exist on the role of non-linear interactions of Kelvin 
waves and the resulting Kelvin-wave cascade, which is ultimately terminated 



in quasiparticle emission (Kozik and Svistunov 2004, 2005a, 2008a, Vinen 
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2000, Vinen et al. 2003), or on the importance of reconnections which could 



rapidly redistribute energy over a range of scales and also lead to dissipation 



(Svistunov 1995). From their front propagation measurements Eltsov et al. 



(2007) conclude that the Kelvin-wave cascade accounts for an important part 



in the increased dissipation below 0.3 T c . The different sources of dissipation 
in this analysis are discussed in Sees. 3.4 and 3.5. 

It is worth noting that a propagating turbulent vortex front has many 
interesting analogues in physics ( van Saarloos||2003 ) . For instance, it is similar 
to the propagation of a flame front in premixed fuel. Flame front propagation 
can also proceed in laminar or in turbulent regimes. In the latter case the 
effective area of the front increases and its propagation speed becomes higher 
than in the laminar regime. This property finds its practical use in combustion 



engines, but has also been used by Blinnikov et al. (2006) to describe intensity 



curves of type la supernovae. In all such metastable state of matter 

is converted to stable state in the front and Vf is determined by the rate of 
dissipation of the released energy. 

3.2. Measurement of vortex front propagation 



The velocity of the vortex front was measured by Eltsov et al. (2007) with 



the setup in Fig. 13. The initial vortex-free state was prepared by warming 
the sample above 0.7 T c , where remanent vortices annihilate rapidly, and by 
then cooling it in the vortex-free state at constant rotation to the target tem- 
perature. Two different procedures were used to trigger the turbulent burst 
at the target temperature. These are sketched in Fig. 13. In both cases the 
front velocity is determined by dividing the flight distance by the flight time, 
assuming that the front propagates in steady-state configuration. Although 
this is not exactly true, for instance owing to initial equilibration processes 
which follow injection, it is assumed for now that this simplification is jus- 
tified. Especially, since the two injection techniques for different propagation 
lengths give the same result. 

The first injection technique (depicted on the left, Fig. 13) makes use of 



remanent vortices (Solntsev et al. 2007). By trial and error it was found that 



one or more remnants can be freed with a small step increase in rotation from 
the region around the orifice on the bottom of the sample cylinder. This was 
done by increasing Q in small steps, until at some point usually above lrad/s 
a remnant starts expanding which below 0.35 T c immediately gives rise to a 
turbulent burst. The ensuing vortex front then propagates upwards along the 
entire column through both pick-up coils in succession. 

The second injection method (depicted on the right, Fig. 13) relies on the 
superfluid Kelvin-Helmholtz (KH) instability of the interface between the A 
and B phases of superfluid 3 He (Blaauwgeers et al. 2002). Two stable AB 



interfaces are formed by applying a specially configured magnetic field which 
stabilizes a narrow A-phase barrier layer over the mid section of the sample 
cylinder. The shear flow instability of these two AB interfaces is controlled 
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Fig. 13. Experimental setup for measuring the propagation velocity of the vortex 
front in the rotating column. Two methods are shown for measuring the front motion 
across different flight lengths in one single experimental setup. On the left the seed 
vortices are tiny remnants at the orifice. In increasing rotation at f2 > 1 rad/s 
they produce a turbulent burst in the volume around the orifice below the bottom 
detector coil. A single vortex front is then observed to pass first through the bottom 
coil and later through the top coil. The time difference separating the signals from 
the passing front over the flight path of 90 mm defines the front velocity V{. On the 
right the seed vortices are injected via the Kelvin-Helmholtz shear flow instability of 
the two AB interfaces. The injection event is followed instantaneously by a turbulent 
burst close to the AB interface on the B-phase side. A vortex front is then observed 
to propagate independently both up and down along the cylinder. The lengths of 
the flight paths are equal for the upper and lower halves. In Fig. 14 it is explained 
how the flight time is determined in this case. 



by rotation velocity, temperature, and the stabilization field. At the target 
temperature the instability can be triggered with a step increase of the rotation 

. The instability causes 
vortices from the A phase to escape across the AB interface into the vortex- 
free B-phase flow in the form of a bunch of small closely packed loops. Once in 
the B phase, at temperatures below 0.59 T c the loops immediately interact and 
generate a turbulent burst. Two vortex fronts then propagate independently 
up and down from the AB interfaces, arriving to the top and bottom pick-up 
coils practically simultaneously (since the setup is symmetric with respect to 



velocity or of the stabilization field (Finne et a/.|2004b 
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Fig. 14. NMR signal responses of the propagating vortex front and the twisted 
vortex state behind it. Here the turbulent burst is started with the Kelvin-Helmholtz 
instability of the AB interface. (Top panel) It is triggered by increasing $7 by a 
small amount AO across the critical value O c AB at 1.3rad/s. (Main panel) Two 
absorption responses are shown, which have been recorded with the bottom detector 
coil at constant, but different values of magnetic field. The responses are from two 
consecutive identical measurements, to allow direct comparison of signal amplitudes. 
The counterflow peak height (thick line) shows the KH trigger A£l and a rapid 
collapse when the front moves through the coil. The time interval from t = (when 
£l(t) = O c ab) to the start of the collapse measures the flight time of the front from 
the AB interface to the closer end of the detector. The moment when the peak 
height reaches zero corresponds to the point when the front has passed through the 
rear end of the coil. The time required for the collapse, ~ tcf, measures the width of 
the vortex front. The second signal (thin line) is recorded close to the Larmor edge 
and is sensitive to the longitudinal velocity v sz which is generated by the twisted 
vortex state (Fig. 21). Its sudden steep rise at t 30s is caused by the passage of 
the first helical sections of the twisted state through the coil. Its later exponential 
decay reflects the unwinding of the twist, which starts when the front has reached 
the end plate of the cylinder and the vortex ends begin to slip along the flat surface. 



the mid plane of the stabilization field). An example of the NMR readout as 
a function of time is shown in Fig. 14. The KH shear flow instability and the 
associated vortex leak across the AB interface have been extensively described 
in the review by Finne et al. (2006b). 

In Fig. 14 one of the signal traces records the absorption at the counterflow 
peak (cf. Fig. 7). It is at maximum in the initial vortex-free flow at v = 
Vn<t> — v S ct> (where v n ^ = Q r and v S( p = in the laboratory frame) . To trigger 
the KH instability, Q is increased by a small increment across the critical 
rotation velocity fi c AB, which is instantaneously registered as a small increase 
in absorption level, owing to the increased counterflow velocity. Following the 
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instability and the turbulent burst, the vortices subsequently propagate along 
the column, but a response in the counterflow peak height is not observed until 
they reach the closer end of the NMR coil. From thereon the absorption in 
the peak rapidly decreases and drops to zero. Keeping in mind that the peak 
height measures the azimuthally circulating flow in the transverse plane (cf. 
Eq. (1)), its sudden removal requires that the passage of the vortices through 
the coil must occur as an organized sharp front, followed by a highly polarized 
state behind the front. The passage is characterized by the time tqf, which is 
defined in Fig. 14. The propagation velocity Vj of the front is determined from 
its flight time, measured from the AB interface instability (when Q(t) = ^cAb) 
to the arrival of the front at the closer edge of the pick-up coil, i.e. where the 
rapid drop in the counterflow peak height starts. Compared to the flight time, 
the AB interface instability and the turbulent burst can be considered as 
instantaneous. 

The second signal trace in Fig. 14 records the absorption in the Larmor 
peak. It is near zero in the initial vortex-free state, displays a sharp maximum 
after the collapse of the counterflow peak, and then decays to a small, but 
finite value, which is a characteristic of the final equilibrium vortex state. The 
transient maximum is the new feature, which arises from flow in the axial 
direction, created by a helically twisted vortex bundle (Eltsov et al. 2006b). 
With increasing wave vector Q of the helix, the axial flow at a velocity v sz and 
the absorption in the Larmor peak increase monotonically ( Kopu|[2~006 ) . Thus 
the maximum Larmor peak height h tw in Fig. 14 is reached when v sz reaches 
its largest value inside the detector coil. This happens when the most tightly 
spiralled section of the twisted cluster (which is just behind the front, as seen 
in Fig. 21) passes through the middle of the detector coil. When the front 
arrives at the bottom end plate of the cylinder the twist starts to relax, since 
the vortices have to obey the boundary condition on the flat end plate, where 
they slip to reduce their length and winding. The unwinding produces the 
exponentially relaxing absorption with time constant TL ar - The signal decay 
continues down to the absorption level characteristic of the equilibrium vortex 
state, with an array of rectilinear vortex lines. At their other end the vortices, 
after crossing the AB interface, have a continuation as doubly quantized A- 
phase vortices (Hanninen et al. 2003). 

To interpret measurements on the front velocity it is of importance that the 
structure of the twisted state and the front itself are known. Information on 
these characteristics can be obtained by measuring quantitatively the various 
features denoted in Fig. 14 as a function of temperature. 



3. 3. Velocity of vortex front 



A striking consequence from the twisted state is the appearance of superflow 
directed along the helically spiralling vortex cores. This situation is reminiscent 
of a "force-free" vortex configuration, where all the flow is directed along the 
vortex core. Such a structure is expected to be stable up to some instability 
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Fig. 15. (Left) Magnitude of the twist as a function of temperature. The mea- 
surements are performed using the bottom spectrometer. Setups 1 and 2 refer to 
measurements with the detector coils positioned in two different sets of positions 
along the sample cylinder; in Fig. 13 setup 2 is depicted, (o, □): The measured ra- 
tio of the maximum amplitude h tw of the Larmor peak in the twisted state to the 
amplitude h^fl = 0) of the Larmor peak in the nonrotating sample is plotted on 
the left axis. (♦): The maximum value of the twist wave vector Q, obtained from 
simulation calculations, is plotted on the right axis. The solid curve shows the fit 
Q R = 0.7(1 — a') /a. The dashed curve shows the minimum Q at which the vortex 
front still propagates in a thin steady-state configuration. (Right) Apparent thick- 
ness of the vortex front TcpVf as a function of temperature, as determined from 
measurements triggered with the Kelvin-Helmholtz instability (Fig. 13). The solid 
line is the prediction of the model in Eq. (11). 



limit, similar to the Glaberson limit of a rectilinear vortex array in parallel flow 
(Sec. 2.4). In the twisted state, which is uniform in the axial and azimuthal 
directions, the superflow has both an azimuthal component at the velocity 
and an axial component at v sz , which depend on the radial coordinate r and 



are described by the expressions (Eltsov et al. 2006b): 



(Q + Qv )r v - Qfir 2 

V ^ r) = i + W 2 ' Vsz{r)= ugv • (9) 

Since the net flow through the cross section of the cylindrical container should 
vanish, from this condition for v sz one finds that vq = (Q / Q)[Q 2 R 2 / ln(l + 
Q 2 R 2 ) — 1]. The axial flow is directed along the vortex expansion direction 
close to the cylindrical wall and in the opposite direction closer to the center. 
In practice, there has to exist also a radial velocity v sr , since any laboratory 
example of the twisted state is nonuniform. In the case of a propagating vortex 
front, the wave vector Q has its maximum value close to the rear end of the 
front and decreases to zero at the bottom and top end plates of the sample. 
As seen in Fig. 14, the twisted state prominently changes the line shape of 
the NMR spectrum and it is the axial superflow which here has the strongest 
influence. 

As discussed in Sec. 2.2, in vortex-free counterflow the end point of a single 
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vortex moves on the cylindrical side wall roughly with the longitudinal velocity 
v Lz m aQR, while its azimuthal velocity is ~ —(1 —a')QR (in the rotating 
frame). Thus the wave vector of the spiral trajectory is Q = Iul^/C^l,*)! ~ 
(1 — a!)/{Ra) = (/R. If this value is used as an estimate for the wave vector 
of the twisted state, it follows that the helical winding of the vortex bundle 
becomes tighter with decreasing temperature. The tighter twist increases the 
flow velocities in Eq. (9), i.e. flow parallel to the vortex cores is enhanced, 
which will ultimately destabilize the "force-free" twisted-cluster configuration. 
Nevertheless, twisted-cluster propagation appears to persist even below 0.2 T c . 
The twist is removed by the slip of the vortex ends along the flat end plates of 
the cylinder, which generates the exponentially relaxing absorption in Fig. 14 
with the time constant T^ ar . 

In Fig. 15 (left panel) the measured temperature dependence of the mag- 
nitude of the twist is plotted in terms of the maximum height of the Larmor 
peak /i tw , normalized to the height of the Larmor peak at Q = at the same 
temperature (left vertical axis). Two experimental setups with slightly dif- 
ferent specifications were used in these measurements. The line shape of the 
NMR absorption spectrum in the Larmor region depends both on the magni- 
tude and homogeneity of the magnetic polarization field. In the two setups the 
homogeneities varied by a factor of two, which is believed to explain the dif- 
ferences in the absolute values between the two data sets (for details we refer 
to Eltsov et al. (2008)). Nevertheless, it is seen here that the twist increases 



towards low temperatures as expected, but only until a maximum at 0.45T C , 
whereas below 0.45 T c it abruptly starts to decrease. 

The non-monotonic temperature dependence of the twist is confirmed in 
numerical calculations: The value of the twist wave vector behind the front, 



as determined from a fit of the calculated velocity profiles to Eq. (9) (Eltsov 



et al. 2008) and plotted in Fig. 15 (left panel, right vertical axis), also peaks 



at 0.45T C . Two reasons can be suggested for the change in temperature depen- 
dence at 0.45 T c . First, the twist can relax via reconnections between neigh- 
boring vortices in the bundle, which become more frequent with decreasing 
temperature below 0.4 T c (cf. Fig. 25). Secondly, the source of the twist is 
at the vortex front, while the sink is at the end plate of the cylinder, where 
the twist vanishes because of the boundary conditions. From there the re- 
laxation of the twist advances in a diffusive manner along the twisted bun- 
dle. The effective diffusion coefficient increases as the temperature decreases 
(Eltsov et al. 2006b) and thus the faster diffusion limits the maximum twist 
in a finite-size sample at low temperatures. However, overall the stability of 
twisted-cluster propagation appears to be a complicated question at tempera- 
tures below 0.3 T c , where it controls the average number of vortices threading 
through each cross section of the column behind the front. 

The properties of the twist in Fig. 15 roughly agree with the estimate 
that the front velocity can be approximated with the longitudinal velocity 
of a single vortex expanding in vortex-free rotation, v^ z ~ aQR. However, 
this simplification suffers from the following difficulty: Ahead of the front 
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the vortex-free superfluid component is at rest and the effective counterfiow 
velocity might really be approximated with v = Qr, but behind the front 
the density of vortices is close to equilibrium and v s ^ rs v nt j,. In this case a 
vortex, which has fallen behind in the motion, feels a much reduced counter- 
flow and continues to fall more behind. Therefore the thickness of the front 
should increase with time. The explanation to this dilemma is that behind 
the front the superflow induced by the twisted vortex bundle has to be taken 
into account. The longitudinal expansion velocity should now be modified to 
vlz = oi [v n< j,(R) — v S cj>(R)] + (1 — a')v sz (R). Since here v sz (R) is oriented in the 
direction of the front propagation and v s ^(R) < v n <f,(R) in the twisted state, 
the longitudinal expansion velocity V t of the vortices in the tail of the front 
is enhanced. This velocity can be estimated taking v sz (R) and v s ^(R) from 
Eq. (9): 



V t = aQR 



a' 



a QR 



Q 2 R 2 



1 - 



1 + Q 2 R 2 ) log(l + Q 2 R 2 ) 



(10) 



Eq. (10) has a maximum as a function of the Q vector. If (1 — ac')/a < 1.9 



(which corresponds to T > 0.46 T c according to the measurements of Bevan et 



\al\ (1997)), the maximum value of V t is less than the velocity of the foremost 
vortices Vf ~ aQR. In these conditions the thickness of the front increases 
while it propagates. When T < 0.46 T c , a wide range of Q values exists for 
which formally V t ^ Vf. The minimum possible value of Q is shown in Fig. 15 
(left panel) as the dashed curve. In these conditions the front propagates in a 
steady-state "thin" configuration. 

Experimentally, the decay time of the counterfiow peak tcf in Fig. 14 can 
be used to extract the front thickness. The decay starts when the head of the 
front arrives at the closer edge of the detector coil and it is over when that 
part of the front leaves the far edge of the detector coil where the counter- 
flow is not sufficient to generate a non-zero absorption response. The product 
tcf Vi has the dimension of length and can be called the apparent thickness of 
the front. At higher temperatures the actual thickness of the front grows with 
time. Here the apparent thickness depends on the distance of the observation 
point from the site of the turbulent burst and on the rate at which the thick- 
ness increases and vortices fall behind. With decreasing temperature the front 
starts to propagate as a thin steady-state structure and ultimately its appar- 
ent thickness decreases to equal the height of the pick-up coil (h c = 9 mm, 
Fig. 7) and remains thereafter approximately constant. 

Measurements of the apparent thickness of the front are presented in the 
right panel of Fig. 15. At T > 0.45T C , rcF^f > h c and the apparent thickness 
increases with increasing temperature. At 0.45 T c the apparent front thickness 
becomes comparable with the height of the detector coil and thereafter at lower 
temperatures remains at that value. Assuming that initially at the turbulent 
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Time (s) 



Fig. 16. Measurement of vortex front propagation. The NMR absorption in the 
counterfiow peak is monitored as a function of time, after triggering the Kelv- 
in- Helmholtz instability at t = (c/. Fig. 14). The instability is started in constant 
conditions in the vortex-free state at 1.2rad/s, by increasing 0, in one small step 
above the critical value ^ C AB) when a magnetic field stabilized A-phase layer is 
present (topmost panel), or by increasing the magnetic field stepwise above Hab 
at constant O (two lower panels). The two signal traces denote the top (thin noisy 
line) and the bottom (thick line) detector. Since the flight paths for the upper and 
lower sample sections are almost equal, the two traces display almost identical flight 
times. 



burst the front is infinitely thin we can write 

h h + h c h h 

where h h is the distance from the site of the turbulent burst to the nearest edge 
of the pick-up coil and is the expansion velocity at the position in the front 
where the NMR signal from the counterfiow vanishes. Given that the latter 
condition roughly corresponds to v s $ ~ (l/2)t> n we take V t * = (V t + Vf)/2 if 
V t < Vf and simply V* = Vf otherwise. Using V t from Eq. (10) and the simple 
estimates QR = (1 — a') /a and Vf = aQR, we get from Eq. (11) the solid line 
in Fig. 15 (right panel), which is in reasonable agreement with experiment. 

The rapid change in the counterfiow peak height during the passage of 
the vortex front through the detector coil provides a convenient signal for 
measuring the propagation velocity Vf. Examples of these signals are shown in 
Fig. 16. They have been measured at different temperatures to illustrate how 
the temperature dependence of Vf is expressed in the practical measurement. 
All three examples have been measured using the externally triggered Kelvin- 
Helmholtz instability to start the turbulent burst. In Fig. 16 the vortex fronts 
have traveled a distance of ~ 4 cm, before they pass through the detector coil, 
and thus have already acquired their steady-state thin-front configuration. 

Measurements on the front velocity Vf are shown in Fig. 17. As a function of 
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Fig. 17. Normalized velocity V{/£IR of vortex front propagation in a rotating col- 
umn. All externally controlled variables are kept constant during this measurement 
where the initial state is vortex-free rotation and the final state an equilibrium ar- 
ray of rectilinear vortices. (Main panel) The assignments of different hydrodynamic 
regimes refer to the dynamics in the front motion. The open circles denote measure- 
ments in which the turbulent burst is started from the AB interfaces in the middle 
of the column. The squares refer to the case where the front is started at the orifice 
and then moves upward through the entire column. The large filled diamonds mark 
results from numerical calculations. The dashed line represents the mutual friction 



dissipation a(T) measured by Bevan et al. (1997) and extrapolated below 0.35 T c 
with exp(-A/T) ( |Todoshchenko et al.\\2002\ . The dash-dotted curve takes into ac- 
count the twisted vortex state behind the propagating front (Eq. (12)). The solid 
curve displays the theoretical model from Sec. 3.5 which includes corrections from 
the twisted vortex state, turbulent energy transfer, and quantum bottleneck. (Left 
panel) This semi-log plot shows that the analytically calculated model provides a 
reasonable fit to the low temperature data. 



temperature two different regimes of front propagation can be distinguished, 
the laminar and turbulent regimes. The crossover between them is gradual and 
smooth. This is in sharp contrast to the sudden onset of bulk turbulence (as a 
function of temperature around £ ~ 1) in injection measurements in the same 
circular column, when a bundle of closely spaced seed vortex loops escapes 
across the AB interface in a Kelvin- Helmhoits instability event (Fig. 2). A 
sharp transition with a clearly denned critical velocity is the usual case, for 
instance in all measurements with mechanical vibrating objects (in the regime 
( > 1) as a function of drive and flow velocity (Vinen and Skrbek 2008). 
The smooth crossover here in front propagation may be a special property 
of the circular column where vortex polarization along the rotation axis is 
always > 90 %. Nevertheless, smooth crossovers have been observed before, for 
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instance from linear to turbulent wave acoustics in second sound propagation 



in a circular cylinder as a function of driving amplitude (Kolmakov et al. 



2006). 



Above 0.4T C in the laminar regime the results in Fig. 17 are consistent with 
the earlier measurements of Finne et al. (2004c). Here the single- vortex dy- 
namics apply, when inter-vortex interactions can be neglected and Vf ~ aQR 
(Sec. 2.2). At closer inspection it is noticed that the data for the normalized 
front velocity Vf = Vf/QR lie on average below a(T) (dashed curve). The 
reason is that behind the front the vortices are in the twisted state and not 
as rectilinear vortex lines in an equilibrium array. If we integrate the kinetic 
energy stored in the twist-induced flow (Eq. 9) over the sample cross section, 
the result is found to be 



t; filam = 2[l/log(l + C 2 )-CV- (12) 

After including this reduction from the twisted state, the agreement is im- 
proved (dash-dotted curve). In the turbulent regime below 0.4T C , the data 
deviate with decreasing temperature more and more above the extrapolations 
from the laminar regime. Eventually at the lowest temperatures the measure- 
ments become temperature independent, with a peculiar transition from one 
plateau to another at around 0.25 T c . These features are attributed to turbu- 
lent dynamics and are analyzed in the next sections in more detail. 

The measured properties of the propagating vortex front are confirmed 
qualitatively in numerical calculations. They show that in the laminar regime 
the thickness of the front grows with time, but with decreasing temperature 
the twist increases (Fig. 15) and finally at about 0.45 T c the thin steady-state 
front configuration is established, with a time-independent thickness roughly 
equal to the radius of the sample. Below 0.45 T c the twist decreases with de- 
creasing temperature, but remains within the limits where the twist-induced 
superflow is sufficient to maintain a thin time-invariant front configuration. 
The calculated front velocity in Fig. 17 approximately agrees with the mea- 
surements down to about 0.22 T c . It is therefore instructive to analyze the 
calculations, to identify where and by what mechanisms turbulent losses oc- 
cur in the rotating column. This will be discussed in the next section. 



3.4- Numerical calculation of turbulence in vortex front propagation 

The promising agreement of the calculated vortex front velocity Vf with the 
measurements in Fig. 17 suggests that numerical calculations could provide 
useful guidance for the interpretation of the measurements in Sec. 3.3 and for 
constructing the analytic theory in Sec. 3.5. Our main question is the following: 
Why are the numerical results on Vf deviating with decreasing temperature 
more and more above the mutual-friction controlled extrapolation from the 
laminar regime? For simplicity we split the discussion of vortex motions to 
three different length scales: (i) large scales of order ~ R, where turbulent 
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Fig. 18. Calculated vortex propagation at 0.3 T c and 0.4 T c . The motion is started 
from the bottom end of the cylinder, by placing the equilibrium number of vortices 
in the form of quarter loops between the bottom end plate and the cylindrical side 
wall. On the right at 0.4 T c the front has traveled for 80s to a height z ~ 28mm 
above the bottom end plate in a cylinder of 3 mm diameter rotating at 1 rad/s. The 
zoom on the far right shows the vortices in the front and immediately below in 
more detail. On the left at 0.3 T c the same distance is covered in 140 s. Here the 
vortices appear more wrinkled, owing to short-wavelength Kelvin wave excitations. 
The vortices are tightly twisted below the front, but become straighter and smoother 
on approaching the bottom end plate. Many vortices can be seen to connect to the 
cylindrical side wall also below the front. Their average length is shorter than the 
distance from the front to the bottom end plate, although the number of vortices 
N(z) threading through each cross section of the cylinder below the front is roughly 
constant and comparable to that in the equilibrium vortex state: N(z) < N eq . The 
average polarization along the vertical axis is high, ~ 90%. 

fluctuations become visible as variations in the number and distribution of 
vortices, (ii) the inter-vortex scale ~ £, where the presence of Kelvin-waves 
on individual vortex lines can be seen to grow, and (iii) small scales, where 
reconnections between neighboring vortices might occur and excite turbulent 
fluctuations on individual vortex lines. In contrast to the more usual studies 
of turbulent tangles, which generally monitor their free decay as a function 
of time when the external pumping is switched off, we are here dealing with 
highly polarized steady-state motion along the rotating column. This mo- 
tion appears to remain intact down to the lowest temperatures, but becomes 
dressed with turbulent fluctuations in growing amounts towards decreasing 
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Fig. 19. Axial distributions of counterflow energy £(z) (see Eq. (8)) and polarization 
Pz(z) (both in normalized units), calculated for the vortex propagation in Fig. 18. 
On the left £(z) is seen to drop steeply within the narrow vortex front at the tem- 
peratures of 0.4 and 0.3 T c . Here the front is in steady-state time-invariant motion 
and its velocity provides a direct measure of dissipation. Above 0.45 T c some vor- 
tices tend to fall more and more behind during the motion and the shape of the 
front becomes more extended with time. As seen in Fig. 18, small-scale structure 
from Kelvin waves accumulates increasingly on the vortices below 0.45 T c , but as 
shown in the panel on the right, the polarization p z {z) remains always high behind 
the front. 



temperatures. 

Two examples are shown in Fig. 18 which illustrate the development with 
decreasing temperature. These configurations at 0.4 T c and 0.3 T c in an ideal 
column of radius R = 1.5 mm rotating at Q — lrad/s have been calculated 
using the techniques described in Sec. 2.5. Similar results at 0.4 T c and higher 
have been reported in the review by Finne et al. (2006b). The calculations are 
started from an initial configuration with roughly the equilibrium number of 
vortices placed as quarter-loops between the bottom end plate and the side 
wall of the cylinder. During the subsequent evolution the propagating vortex 



front is formed and the twisted cluster starts to acquire its shape (Hanninen 



2006). 



Above 0.45 T c the motion is laminar with relatively smooth vortices, which 
only twist at large length scales. Below 0.45 T c the thickness of the front (in 
the axial direction) settles at A(r) ~ r d, where the parameter d ~ 1. The 
large-scale characteristics as a function of temperature from 0.3 to 0.6 T c are 
displayed in Fig. 19 for the setup of Fig. 18. Here the counterflow energy £{z) is 
obtained by integrating the momentary distribution of counterflow v(r, 0, z, t) 
over each cross section z of the column, while the polarization p z (z) is similarly 
derived by integrating over the vortices threading through this cross section. 
With decreasing temperature the front acquires more and more turbulent fea- 
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T/T c 


0.3 


0.4 


Mutual friction parameter a 


0.040 


0.18 


Mutual friction parameter a' 


0.030 


0.16 


Front velocity Vf 


0.12 OR 


0.22 QR 


Front velocity Vf (mm/s) 


0.18 


0.33 


Total Reconnection rate (event /s) 


300 


130 


Number of vortex lines/cross section 


125 


150 


Reconnection rate per line 1/r [events/(line s)] 


2.4 


0.87 


Front shift 5 during r (mm) 


0.075 


0.38 


Interline separation t (mm) 


0.20 


0.19 


Ratio 5/e 


0.36 


2.0 



Table 1 

Comparison of vortex front propagation at 0.3 T c and 0.4 T c . The bottom line of 
the table shows that below 0.4 T c front propagation is rapidly moving into the 
quantum regime. The calculations are for a cylinder of radius R = 1.5 mm and 
length h = 40 mm, which is filled with 3 He-B at a liquid pressure of P = 29 bar 
and rotates at an angular velocity O = lrad/s. The number of vortex lines is the 
average through each cross section of the cylinder over the length of the twisted 
cluster. The maximum resolution of the calculations is 0.05 mm. 

tures which become visible as increasing small-scale structure: Kelvin waves, 
kinks, and inter-vortex reconnections. These small scale fluctuations exist on 
top of a strongly polarized vortex orientation, which is preserved even at 0.2 T c . 
The high polarization along the rotation axis and the organized configuration 
of the twisted vortex cluster is expected to reduce loop formation and self- 
reconnection on individual vortex lines, but it also suppresses reconnections 
between neighboring vortices from what one would find in turbulent tangles. 
The amount of twist, in turn, is reduced by the unwinding from the slip of the 
vortices along the flat bottom end plate of the column. 

The comparison of the two examples in Fig. 18 is continued further in 
Table 1. As a rule we expect that with decreasing mutual- friction damping 
turbulent disturbances are expected to cascade down to ever smaller length 
scales. In Table 1 this is examined by denning the quantity 5 which measures 
on an average the distance over which the vortex front moves during the 
time interval between two inter-vortex reconnections. Since the axial motion 
of the front slows down and the reconnection rate increases with decreasing 
temperature, 8 decreases rapidly below 0.4 T c . Thus 8 provides a measure 
between mutual friction and reconnections. The bottom line in Table 1 uses 
the dimensionless ratio 8/1 to characterize the relative importance of mutual 
friction dissipation to reconnection losses. When 8/1 > 1, the energy loss from 
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Fig. 20. Axial and radial distributions of the velocity | v s (r, <f>, z,t) | and its fluctu- 
ations in vortex front propagation. Both quantities are normalized with respect to 
£IR and are expressed in the laboratory frame in terms of their averages, {v s )t t( /> (top 
panel) and [((v s — {v s ) tj( j } ) 2 )t,<f>]^ (bottom panel), by integrating over the azimuthal 
coordinate <j> and time t. The different contours are plotted as a function of z at 
fixed radial value, r = 0, R/5, 2 R/5,.... They are calculated for the example at 
0.3 T c in Fig. 18, over the time interval from 60 to 80s, after starting the front 
motion. At t = 80s the front has climbed to a height z(80s) ~ 16mm (here placed 
at z/R = 0). Both plots are generated from vortex configurations, which are saved 
every 0.5 s. The steep change in the top panel in the interval —1 < z/R < signi- 
fies the front with its narrow thickness Az ~ R. Similarly in the bottom panel the 
amplitude of fluctuations (sampled at 2 Hz) reaches a sharp maximum in the front 
region —1 < z/R < at large radii r > 3 R/5. This peak is twice larger than the 
flat values from behind the front in the region of the twisted cluster. At 0.4 T c (the 
case of Fig. 18 right) the fluctuation peak from the front is 14% lower. 

reconnections, which presumably excite Kelvin waves, is not important in the 
total energy balance of the propagation. At 0.4 T c , where 5/£ ~ 2, this is what 
one expects. If S/£ < 1, as is the case below 0.3 T c , then it becomes possible 
that inter-vortex reconnections might play a role in the total energy balance. 
Table 1 thus hints that new phenomena are expected to emerge on length 
scales approaching the inter-vortex distance i. 

The characteristics of the propagating vortex front and the trailing twisted 
cluster are further illustrated in Figs. 20 - 22. In Fig. 20 it is the modulus 
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Fig. 21. Axial and radial distributions of the axial velocity component v sz (r, 0, z, t) 
and its fluctuations. The profiles depict the same calculations at 0.3 T c as Fig. 20. 
This component is generated by the helically wound sections in the twisted cluster. 
The slowly unwinding twist at the bottom end plate (on the left of the figure) 
causes the characteristic linear increase in | v sz | towards the right where new twist 
is continuously formed by the spirally winding motion of the vortex front. Note 
that v sz changes sign, since it is directed antiparallel to the propagation direction 
of the front in the central parts r < 2 R/3 and parallel at larger radii. Its maximum 
magnitude is in the center just behind the front where it is the dominant component 
in | v s | . Its largest fluctuations are at large radii within the front region. 

of the total velocity of the superfluid component, | v s (r, (ft, z,t) |, whose [r,z] 
profiles are displayed. The top panel shows the velocity (f s )t,<?i m the labora- 
tory coordinate frame, averaged over time t and azimuthal angle 0, while the 
bottom panel displays its mean fluctuation amplitude [((v s — (f s )t,^) 2 )t,</>] 5 - In 
Fig. 21 it is the axial component v sz and in Fig. 22 the radial component v sr 
which are analyzed in the same fashion. The largest component is omitted, 
since its profiles are so similar to those of | v s |. In fact, at large radii r > R/5 
the profiles of the two velocities as a function of [r, z] are almost indistinguish- 
able, while at small radii r < R/5 the axial velocity v sz in Fig. 21 is the main 
contributor to \v s \. 

The vortex front becomes clearly denned in Figs. 20 - 22. For instance, 
in Fig. 20 the steep almost linear rise in | r? s | in the interval —1 < z/R < 
signifies the transition from the non-rotating state = at z > to 
almost equilibrium rotation at v S( f, w fir at z < —R. This is the vortex front 
with its narrow thickness ~ R and strong shear flow, created by the vortices 
terminating within the front on the cylindrical side wall perpendicular to the 
axis of rotation. Below the front the number of vortices threading through 
any cross section of the cylinder remains roughly constant. This is seen from 
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Fig. 22. Axial and radial distributions of the radial velocity component v sr (r, 0, z, t) 
and its fluctuations in vortex front propagation at 0.3 T c in Fig. 18. This component 
provides the return currents to the axial component v sz . It is an order of magnitude 
smaller than either the azimuthal v S( f, or the axial v sz component. It is directed 
inward towards the central axis in the narrow front region, while the outward flowing 
radial return currents are distributed more evenly along the length of the twisted 
cluster. The fluctuations in v sr are largest close to the cylinder axis and steadily 
decrease as a function of r. As a function of z the maximum fluctuations are just 
behind the front where v sr changes sign. 

the fact that at 0.3 T c the maximum value of v s , (v sr p) maiX ~ 0.6 QR remains 
stable over the entire length of the twisted cluster, starting from z/R < — 1 
immediately behind the front. At higher temperatures some vortices tend to 
fall behind the front and thus at 0.4 T c one finds that v s slowly increases to 
{v S (/>)ma,x ~ 0.75 QR at z/R < —5. Note also that radially the azimuthal flow 
Vstj) increases monotonously behind the front up to the edge of the cluster at 
~ 0.9 R and then slightly decreases towards the cylindrical wall. 

The fluctuations of | v s | around its mean value, as shown in the bottom 
panel of Fig. 20, are substantial, of order 30 % in the sharp peak created by the 
vortex front. The same applies to other velocity components, their fluctuations 
also tend to be largest in the front region — 1 < z/R < and typically only half 
as large over the length of the twisted cluster. In Figs. 20 - 22 the fluctuations 
are sampled at a relatively slow rate of 2 Hz. As seen in Fig. 18, in the front 
the vortices are not perfectly distributed: It is this disorder in the structure of 
the front and the variations in its number of vortices which gives rise to the 
fluctuation peak. In Fig. 23 the radial distribution of the fluctuations in \ v s \ 
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Fig. 23. Radial distribution of velocity fluctuations within the propagating vortex 
front from Fig. 20. The normalized square of the velocity deviation from the average 
is plotted on the vertical scale, ((v s — (v s )t ! <f>) 2 }t,<f>/(&'R) 2 , versus the square of the 
normalized radial position, (r/R) 2 . The averages in the angled brackets, besides 
having been integrated over time t and azimuthal coordinate (or front motion), 
represent an average over the interval — 1 < z/R < 0. The fluctuations are seen to 
increase rapidly towards large radii, up to where the vortex-free annulus starts at 
R v « 0.87 R. 



in the front region — l<^/i?<0is analyzed. As expected, they grow rapidly 
towards large radii. In Sec. 3.5 we make use of this radial distribution. 

However, not only large-scale disorder contributes to velocity fluctuations, 
but also Kelvin waves start to expand on individual vortices at temperatures 
< 0.3 T c , as seen in Fig. 18 (Hanninen 2006). In Fig. 24 a momentary vortex 
configuration (when the front has reached the height z ~ 28 mm in Fig. 18) is 
broken down in curvature radii 1Z C and the average curvature {Ti^ 1 ) is plotted 
as a function of z. At T > 0.4 T c the sharp peak at the front is caused by the 
vortices curving to the cylindrical side wall with TZ C « R. However, at 0.3 T c 
the Kelvin wave contribution at shorter length scales becomes dominant and 
the average radius of curvature drops to (1Z C ) ~ 0.4 mm. The characteristic 
Kelvin wave frequency ~ nk 2 corresponds at this length scale to 0.4 Hz 
(which is less than the sampling frequency of 2 Hz) . Thus the dominant Kelvin 
waves are included in the velocity fluctuations in Figs. 20 - 22. Nevertheless, in 
these figures Kelvin waves represent only a small part of the total fluctuations 
((v B — ( v s)t,<f>) 2 )t,<i> m the sharp peak created by the front. This is seen from 
the fact that the peak is localized in the z direction within the vortex front 
at — 1 < z/R < 0, while the Kelvin-wave excitations in Fig. 24 extend half 
way down the twisted cluster to z/R ~ — 10, as can also be directly seen from 
Fig. 18. 

Fig. 25 shows that also reconnections between neighboring vortices become 
rapidly more frequent at temperatures < 0.3 T c (see also Hanninen ( 2006[ )). 
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Fig. 24. Distribution of curvature radii along vortices during front propagation, 
calculated for the setup in Fig. 18 . The curvature is defined as the inverse of the local 
vortical curvature radius 1Z C , i.e. = \d 2 s(^,t)/d^ 2 \, where £ is the coordinate 
along the vortex core. The curvature is calculated on an equidistant linear grid 
along each line vortex. On the vertical axis the average curvature (T*^ 1 ) is plotted 
as a dimensionless quantity (7£ c _1 ) R as a function of z at different temperatures. 
The averages are taken over the cross section of the cylinder and over an interval 
Az = 2 mm in length centered at z. The sharp peak at the very front (z ~ —2 mm) 
at temperatures T > 0.4 T c represents the vortices curving to the cylindrical side 
wall with 1Z C ~ R. In contrast, at 0.3 T c small-scale curvature dominates, the average 
radius of curvature has dropped to 1Z C ~ 0.4 mm and extends from the front well 
inside the twisted cluster. 



The reconnections peak in the region around z/R < —1, where the front 
ends and the twisted cluster starts. The increasing reconnection rate with 
decreasing temperature is one factor which helps to reduce the helical twist 
behind the front and might be one reason why the measured twist appears to 
be reduced towards low temperatures below 0.45 T c (Fig. 15). 

A further characteristic of front propagation in the rotating column is that 
fluctuations in vortex length are small, since these are mainly restricted to the 
transverse plane. The total vortex length L(t) in the moving front and in the 
twisted cluster behind it increases linearly in time and thus the fluctuations 
can be expressed with respect to a linearly increasing fitted average L &t (t). 
The average deviation from the mean, (AL) = [((L — Lfit) 2 )] 5 turns out to 
be small: In the conditions of Fig. 18 at 0.3 T c one obtains AL/L ~ 10 -3 . 
This is as expected since fluctuations in length are energetically expensive in 
polarized vortex motion. 

Finally, we point out two features which are seen in the calculations, which 
both grow in prominence with decreasing temperature below 0.3 T c , but which 
have not yet been searched for in measurements. The calculations are started 
with the number of seed vortices equal to or larger than N eq , but nevertheless, 
the number of vortex lines does not remain stable at N eq during steady-state 
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Fig. 25. (Main •panel) Cumulative number of all vortex reconnection events during 
vortex front propagation in the setup of Fig. 18, as a function of time. The re- 
connections occur between two different vortices and mainly in the twisted cluster 
where they reach a maximum close behind the front. This process does not change 
the number of vortices. The vertical arrows indicate the moment when the front 
reaches the top end plate (at height h = 40 mm). The reconnections are here seen 
to increase rapidly at low temperatures, reaching a rate of 300 s -1 or ~ 1 recon- 
nection/ (vortex sec) at 0.3 T c . (Inset) The roughly constant reconnection rate from 
the main panel (when the front is in steady-state propagation before the end plate 
is reached), plotted as a function of temperature. 



propagation^"] 

The first feature concerns the number of vortices which in free steady-state 
propagation thread through any cross section of the column behind the front. 
This number is roughly constant, but less than in equilibrium (for instance in 
Fig. 18 at 0.3 T c , there are about 130 vortices per cross section while N eq ps 160 
vortices). In the simplest model of the twist in Eq. (9) the number of vortices 
in the twisted state is smaller than in the equilibrium vortex state by the factor 
« (QR) 2 / ([1 + (QR) 2 ] In [1 + (QR) 2 ])- This reduction is of correct order of 
magnitude compared to that seen in the calculations at > 0.3 T c . In addition 
to improving the stability of the twisted state, the reduction in the number of 
twisted vortices decreases the axial flow velocity v sz (Fig. 21), which has the 
effect of reducing the longitudinal front velocity Vf and counteracts its increase 
above the laminar extrapolation in Fig. 17. While the number of vortices 
per cross section of the twisted cluster decreases with decreasing temperature 
below 0.3 T c , the density of vortices, nevertheless, remains constant over the 
cross section (but less than the equilibrium value: n v < n Vieq = When 
the front finally reaches the upper end plate in Fig. 18 and the twist starts to 
relax, simultaneously N gradually recovers and approaches N cq from below. 



In the equilibrium vortex state iV e q ~ n(R — d eq ) 2 n v , where the width of the 



— 1/2 

vortex-free annulus around the central cluster is d eq ~ n v (jRuutu et al. 



1998b). 
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The second feature of the calculations is that in steady state propagation 
the average length of the vortices is less than that of the twisted cluster (which 
means that the total number of individual vortices in the column may be well 
above N eq ). For instance at 0.3 T c in Fig. 18, the average length is ~ 15 mm 
and only ~ 6 mm for those vortices with both ends on the cylindrical side wall. 
Nevertheless, the polarization of all vortices along z is high, ~ 90 %. After the 
front has reached the upper end plate in Fig. 18 and the twist starts relaxing, 
the short vortices provide an ample storage from which to add more vortices 
so that iVeq is reached. 

To summarize the calculations, we return to our starting point, namely 
the increasing deviation of the calculated front velocity above the laminar 
extrapolation below 0.4 T c (Fig. 17). The analysis of the calculated vortex 
configurations shows increasing turbulent disorder at large length scale ~ R, 
growing Kelvin-wave amplitudes on individual vortices on scales ~ i and the 
presence of inter-vortex reconnections. Combined these changes from the in- 
creasing turbulent influence with decreasing temperature make up for the 
difference, in the presence of a still finite mutual friction dissipation. An ad- 
ditional sink of energy is the cascading of Kelvin-wave excitations below the 
resolution limit of the numerical calculations (usually ~ 0.05 mm), which are 
lost from the energy balance. The calculations are in reasonable agreement 
with measurements down to 0.3 T c , but at lower temperatures changes in the 
twisted cluster propagation appear to occur, to maintain stability and polar- 
ization. These have not been adequately studied and may lead to revisions of 
the current model at the very lowest temperatures. 



3. 5. Analytical model of turbulent front 



The measured front velocity in Fig. 17 displays below 0.4 T c two plateaus 
which are separated by a transition centered at 0.25 T c . The numerical anal- 
ysis in Sec. 3.4 places this transition in the temperature region where the 
Kelvin-wave cascade acquires growing importance, i.e. in the regime where 
sub-intervortex scales start to influence energy transfer. The peculiar shape 



of the measured ff(T) curve prompted L'vov et al. (2007a) to examine en- 
ergy transfer from the quasi-classical length scales (where bundles of vortices 
form eddies of varying size with a Kolmogorov spectrum) to the quantum 
regime (where Kelvin waves expand on individual vortex lines). The result is 



a bottleneck model (L'vov et al. 2008) which matches energy transfer across 



the crossover region from super-intervortex length scales to sub-intervortex 
scales. It can be fitted with reasonable parameters to the measured Vf(T) 



curve (Eltsov et al. 2007) 



3.5.1. Two-fluid coarse-grained equations 

To describe the superfluid motions on length scales exceeding the inter- 
vortex distance I we make use of the coarse-grained hydrodynamic equation for 
the superfluid component. This equation is obtained from the Euler equation 
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for the superfluid velocity U = U s , after averaging over the vortex lines which 
are assumed to be locally approximately aligned, forming vortex bundles which 



mimic the eddies in viscous flow (see the review by Sonin (1987)), 



dU 



+ (U- V)E/ + V//=/, 



mf • 



(13) 



Here fi is the chemical potential and f m f the mutual friction force (3), which 
can be rewritten as follows: 



mf 



-a'(U - U n ) xw + awx[wx(U - U n )] 



(14) 



We use w = V x U for the superfluid vorticity; uj = oj/u is the unit vec- 
tor in the direction of the mean vorticity; U n is the velocity of the normal 
component (which is fixed). In flow with locally roughly aligned vortices the 
mutual friction parameters define the reactive (oc a') and dissipative (oc a) 
forces acting on a bunch of vortex lines, as it moves with respect to the normal 
component in a field of slowly varying vortex orientation. 

We shall work in the rotating reference frame where U n = 0. In this frame, 
the reactive first term in Eq. (14) renormalizes the inertial term U x uj on 
the left hand side (LHS) of Eq. (13), introducing the factor 1 — a' (Finne et 



Z pOOlj ). The relative magnitude of the two non-linear terms, the ratio of the 
inertial term and of the friction term [the latter is the second term on the 
RHS of Eq. (14)], is the Reynolds number in this hydrodynamics. It proves 
to be the flow velocity independent ratio of the dimensionless mutual friction 
parameters: ( = (1 — a') /a (which was introduced by Finne et al. (2003) in 
the form q = 1/C)- To arrive at a qualitative description of the front we follow 



L'vov et al. (2004) and simplify the vectorial structure of the dissipation term 



and average Eq. (14) over the directions of the vorticity uj (at fixed direction of 
the applied counterflow velocity v). With the same level of accuracy, omitting 
a factor 2/3 in the result, we get 



fmf (/mf) w /| w | = ~~aUJ eS V 
LU cS = \J{\U}\ 2 ) . 



au cS v 



(15) 



Here and henceforth (...) denotes "ensemble averaging" in the proper theo- 
retical meaning. From the experimental point of view (...) can be considered 
as time averaging in a frame which moves with the front (over a time interval, 
during which the front propagates a distance equal to its width). The result- 
ing mean values, to be discussed in Eqs. (16) and (19), can be considered as 



functions of (slow) time and space. We follow L'vov et al. (2004) and neglect 
the fluctuating turbulent part of u} c g in Eq. (15). We thus replace \u\ by its 
mean value, the vorticity from rotation. In this approximation Eq. (15) takes 
the simple form: 

fmf = -r U , r = a uj e g . (16) 
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3.5.2. Flow geometry, boundary conditions and Reynolds decomposition 

a. Axial flow geometry and flat idealization. As a simplified model of 
the rotating cylinder we consider first the flat geometry with x as the stream- 
wise, y as the cross-stream, and z as the front-normal directions. A more 
realistic axial symmetry will be discussed later. We denote the position of the 
front with Z(t) and look for a stationary state of front propagation at constant 
velocity —Vf. 

Z{t) = -V f t . (17) 

The normal velocity is zero everywhere, U n (r,t) = 0, which corresponds 
to co-rotation with the rotating cylinder. The vortex front propagates in the 
region z < Zit\ where the superfluid velocity is — V^, while far away behind 
the front at z > Z(t) the superfluid velocity U tends to zero (i.e. towards the 
same value as the normal fluid velocity). 



b. Reynolds decomposition and definitions of the model. Following 



the customary tradition [see, e.g. Pope ( |2000[ )], we decompose the total veloc- 
ity field U into its mean part V and the turbulent fluctuations v with zero 
mean: 

U = V + v, V=(U), (v) = Q. (18a) 

In our model 

V = HV, V^Qr. (18b) 

The following mean values are needed: the mean velocity shear S(r,t), the 
turbulent kinetic energy density (per unite mass) K(r,t) and the Reynolds 
stress W(r,t). Their definitions are as follows: 



S(r,t) 



dV 



K(r,t) = U\v\ 2 ), W(r,t) = -{v x v z ) 



dz 7 '"' ~ 2 

Also we will be using the kinetic energy density of the mean flow, 



K v {r,t) 



l\V(r,tf 



and the total density of the kinetic energy 

)C{r,t)=K v {r,t) + K{r,t) . 
Clearly, the total kinetic energy in the system, £(t), is: 



S{t) = J JC(r,t)dr . 



(19a) 



(19b) 



;i9c) 



(20) 



3.5.3. Balance of mechanical momentum and kinetic energy 



a. Balance of mechanical momentum. Averaging Eq. (13), with the 
dissipation term in the form of Eq. (16), one gets for the planar geometry: 
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dV dW 

w -^ + a^V = 0. (21) 

b. Mean-flow kinetic energy balance. Multiplying Eq. (21) with V one 
gets the balance for the kinetic energy of the mean flow, K v , defined in 
Eq. (19b): 

8K d r n 

—^-—VW+SW + auj cS K v = 0. (22) 
at oz l J 

The second term on the LHS of this equation describes the spatial energy 
flux and does not contribute to the global energy balance of the entire rotating 
cylinder. The next term, S W, is responsible for the energy transfer from the 
mean flow to the turbulent subsystem; this energy finally dissipates into heat. 
The last term, au> e f[K v , represents direct dissipation of the mean flow kinetic 
energy into heat via mutual friction. 

c. Turbulent kinetic energy balance. Let us take Eq. (13) for the ve- 
locity fluctuations, with the dissipation term in the form Eq. (16), and we 
get an equation for the energy balance which involves triple correlation func- 
tions. It describes the energy flux in (physical) space. In the theory of wall 
bounded turbulence ( Pope|[2000 ) these triple correlations are traditionally ap- 



proximated with second order correlation functions. For more details about 



the closure procedures we refer to L'vov et al. (2006a, 2007b) and references 
cited therein. 

To be brief, notice that the total rate of kinetic energy dissipation in the 
vortex front has two contributions in well developed turbulence, 

£diss = £diss,l + ^diss,2 • (23a) 

The first term, £diss,l> arises from mutual friction which acts on the global 
scale. It can be estimated from Eq. (13) as 

^diss,l 

~ au eS K(z,r) , (23b) 

where K(z,r) is the turbulent kinetic energy from Eq. (19a). The second con- 
tribution, £diss,2, originates from the energy cascade toward small scales, where 
the actual dissipation occurs. In well developed turbulence of viscous normal 
fluids this dissipation is caused by viscosity. It dominates at the smallest length 
scales, known as the Kolmogorov microscope. In superfluid turbulence the vis- 
cous contribution is absent. Instead, at moderate temperatures it is replaced 
by mutual friction. When mutual friction becomes negligibly small at the 
lowest temperatures, the turbulent energy, while cascading down to smaller 
length scales, is accumulated into Kelvin waves at some crossover scale. Ulti- 
mately the energy then cascades further down where it can dissipate, e.g. by 
the emission of excitations (phonos in 4 He II and quasiparticles in 3 He-B). In 
any case, the energy has to be delivered to small length scale motions owing 
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to the nonlinearity of Eq. (13). Clearly in steady state conditions, the rate 
of energy dissipation at small scales, £diss,2, is equal to the energy flux, Eq ux , 
from the largest scales, where the energy is pumped into the system [from 
the mean flow via the shear S, as follows from Eq. (22)]. The nonlinearity of 
Eq. (13) is of standard hydrodynamic form and the associated energy flux can 
be estimated by dimensional reasoning, as suggested by Kolmogorov in 1941 



see e.g. Pope (2000) 



e disSt2 = e fiux ~bK 3 / 2 (r)/£(r). 



(23c) 



Here b is a dimensionless parameter, which in the case of classical (viscous) 



wall bounded turbulence can be estimated as 6 c i ~ 0.27 [see e.g. L'vov et al. 



(2006a)]. In Eq. (23c) C(r) is the outer scale of turbulence (which defines the 
length scale of the largest eddies containing the main part of the turbulent 
kinetic energy). Clearly near the centerline of the cylinder, C(r) is determined 
by the thickness A(r) of the turbulent front at given radius r: C(r) ~ A(r). 
Near the wall of the cylinder the size of the largest eddies is limited by the 
distance to the wall, R — r. Therefore in this region, £ ~ R — r. In the whole 
cylinder, C should be the smaller of these two scales, so that one can use an 
interpolation formula 



c 



-1/ 



Air 



(R-r) 



(23d) 



The resulting energy balance equation, which accounts for the energy dis- 
sipation in Eqs. (23), can be written (in cylindrical coordinates r, z, and ip) 
as follows: 



8K bK 3 ' 2 
— + a u ca K H — 

Ot L 



(24) 
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Here K, W, and S are defined in Eqs. (19a). Generally speaking, they are 
functions of r, z, and time t. The outer scale of turbulence, £ depends on r 
and t in our approximation (23d). 

The right-hand-side (RHS) of Eq. (24) represents the energy flux (per unite 
mass) from the mean flow to the turbulent subsystem. This expression rigor- 
ously follows from Eq. (13) and is exact. The two terms on the LHS of Eq. (24) 
(proportional to a and b) describe the energy dissipation (23). The last two 
terms on the RHS of Eq. (24) are proportional to the phenomenological di- 
mensionless parameter g w 0.25, which was estimated by L'vov et al. (2007b). 



These terms model turbulent diffusion processes (in the differential approxima- 
tion) with the effective turbulent diffusion parameter L g \/K, now expressed 
in cylindrical coordinates. The first term in parenthesis describes turbulent 
diffusion in the direction normal to the front, which is the main reason for the 
propagation of the turbulent front. The second term, which vanishes in the flat 
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geometry, describes the turbulent energy flux in the radial direction toward 
the centerline of the sample cylinder. This term plays an important role in 
the propagation of the front as a whole. The reason is that the central part of 
the front, where the direct conversion of the mean flow energy to turbulence 
is small, does not contain enough energy to exist without the radial flux of 
energy. 



d. Total kinetic energy balance. Adding dK v /dt from Eq. (22) and 
dKj dt from Eq. (24) one can see that the term S W, responsible for the 
transfer of energy from the mean flow to the turbulent flow, cancels and we 
arrive at a balance equation for K, = K v + K: 



OK „ bK 3 l' 2 a „„„, 



(25a) 
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Integrating this relation over the cylinder, we see that the energy diffusion 
terms do not contribute to the total energy balance: 



S(t) = J K(r,t)dr 



dS 



- / dr 



auj eB /C 



bK 3 / 2 
C 



(25b) 
(25c) 



To summarize this section, we note that Eqs. (21) and (24) allow us, at 
least in principle, (i) to describe the propagating turbulent front in a rotating 
superfluid, (ii) to find the front velocity Vf, and (iii) to describe the structure 
of the front: its effective width A(r), the profiles for the mean shear S(r,z), 
the Reynolds stress W(r,z), and the kinetic energy K(r,z). Here we present 
only some preliminary steps in this direction, based on a qualitative analysis 
in Sec. 3.5.4 of the global energy balance (25a), and discuss the role of the 
radial turbulent diffusion of energy in Sec. 3.5.5. In Sec. 3.5.6 we finally specify 
the model, accounting explicitly for the bottleneck at the classical-quantum 



crossover, as described by L'vov et al. (2007a), and explain how it is influenced 



by the temperature dependent mutual friction. 



3.5. 4- Qualitative analysis of global energy balance 



The total kinetic energy £ is dissipated by the propagation of the turbu- 
lent front at constant velocity Vf, as expressed in Eq. (8). This means that 
during time t the length of the cylinder with unperturbed mean flow V = Q r 
decreases by Vft. Using Eq. (8) we can write the overall energy budget (25c) 
as: 



V f 



4 

7rQ 2 R* 



dr 



au eS )C(z,r) 



z,r) 



C(r) 



(26a) 
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This equation requires some corrections, because in its derivation we were a 
bit sloppy. First of all, approximation (16) for the mutual friction force / m f 
cannot be applied for the analysis of laminar flow, where the actual orienta- 
tions of vorticity and velocity are important. To fix this we present the RHS 
of Eq. (26a) as a sum of two contributions, 



and for Vf ; i am we use Eq. (12). To get the corrected equation for Vf )tU rb we 
replace in Eq. (26a) JC with its turbulent part K and replace b with b = 6(1— a') 
in order to account for the correction to the nonlinear term in Eq. (13) from 
the reactive part of the mutual friction (14) (proportional to a'). Integrating 
the result over the azimuthal angle ip and making use of the axial symmetry 
of the problem, we get: 



The hydrodynamic equation (13) is not applicable near the wall, where R — r 
is less than the mean inter-vortex distance t. Therefore this region is excluded 
from the integration (26c). 

In the limit of a fully developed turbulent boundary layer (TBL), the tur- 
bulent kinetic energy and Reynolds stress are independent of the distance in 
the z direction to the boundary of turbulence. Here we consider the turbulence 
in the front to be bounded in the z direction, i.e. we assume that the front 
thickness A(r) < R and ignore the influence of the cylindrical container wall. 
Therefore it is reasonable in the qualitative analysis of Eq. (26c) to ignore the 
z dependence of K(z,r) and W(z,r), replacing these objects with their mean 
values across the TBL: 



K(r, z) K(r) , W(r, z) W(r) , oo cS (r, z) cJ cff (r) . (26d) 



As seen from Eq. (26d), we use the same approximation also for uj c q. Now we 
can trivially integrate Eq. (26c) with respect to z: 



Here A(r) is the characteristic width of the TBL at the distance r from the 
central axial. To perform the next integration over r, one needs to know the 
r-dependence of A, uj cS , and K. Notice first that A(r)uJ e s(r) has the dimen- 
sion of velocity. Therefore one expects that in the self-similar regime of a fully 
developed TBL, A(r)uJ e s(r) has to be proportional to the characteristic veloc- 
ity, which is Qr. In other words, we expect that A(r)u c s(r) oc r. As we will 
see below in Sec. 3.5.5, A(r) oc r, while u c s is r-independent. With the same 
kind of reasoning we can conclude that K and W have to be proportional to 



V{ — Vf.lam + turb > 



(26b) 




(26c) 




(26e) 
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r 2 , since they have the dimensionality of velocity squared: 

A(r)uJ cS = attr, K{r) = far) 2 , W(r) = c(fir) 2 . (27) 
Now integrating (26e) one gets 

t* = ^ ~ (2c) 3 / 2 6(l - o')A(i2/€) + ^ , (28a) 
A(i?/£) = 0.2 + d[ \n(R/£) - 137/60 + U/R + ...] , (28b) 

where we used d = A(r)/r. At Q = lrad/s, the ratio R/l ~ 17. This gives 
A(R/£) ~ ^4(17) ~ 1.8. We take 6 = 6 c i and choose the parameters a = 0.2, 
c = 0.25, and d — 2 to fit the measurement in the region (0.3-0.4) T c . With 
these parameters Eq. (28) gives Vf « 0.16 in the limit T — > (when a = a' = 0) 
and a very weak temperature dependence up to T ~ 0.45 T c . 

3. 5. 5. Role of the radial turbulent diffusion of energy 

According to experimental observations, in steady state the front propa- 
gates as a whole with the velocity Vf, independent of the radial position r. To 
make this possible, the turbulent energy has to flow from the near-wall region, 
where the azimuthal mean velocity and consequently the energy influx into 
the turbulent subsystem are large, toward the center, where the influx goes to 
zero. 

To clarify the role of the radial energy flux, consider Eq. (24), averaged in 
the z direction: 

a Z„K(r) + _£» r£(r) ^» W (r) = ^ W(r) . (29) 

On the RHS of this equation we replaced S(r) with its natural estimate 
Qr/A(r). According to Eq. (23d), close to the cylinder axis C(r) 
In this region Eq. (29) has a self-similar scaling solution (27), in which u cS 
is indeed r-independent, A(r) oc r, and for the Reynolds-stress constant c in 
Eq. (27) one finds the relationship: 



c = u^Ai^ b fc ^ 
c n R 2V2 



A(R) 



R 



(30a) 



In the vicinity of the wall, where R — r R, the outer scale of turbulence 
is C ~ R — r and goes to zero when r — > i?. To account for the influence of 
the cylindrical wall in this region we assume that, similar to a classical TBL, 
the two pair velocity correlations W(r) and K(r) have the same dependence 
on distance from the wall. Therefore their ratio is approximately constant: 
W(r)/K(r) ks const. If so, the energy balance Eq. (24) dictates 



W(r) ~ K(r) ~ (R-r) 2 , 



(30b) 
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Fig. 26. Scaled front velocity vs temperature, similar to Fig. 17, but emphasizing 
the comparison of measurements in Sec. 3.3 to analytical results in Sec. 3.5 in the 
turbulent regime at low temperatures. The thin and thick solid lines show consec- 
utive model approximations which sequentially account for dissipation in turbulent 
energy transfer Eq. (28b) and bottleneck effect in Sec. 3.5.6. (Insert) Value of pa- 
rameter b(T) in Eq. (28), which was used for the thick solid line in the main panel. 

with some relation between constants, similar to Eq. (30a). Actually, when R— 
r becomes smaller than the inter-vortex distance £, the whole hydrodynamic 
approach (13) together with Eq. (29) fails. Therefore we cannot expect that 
W(r) and K(r) really go to zero when r — > R. Rather they should approach 
some constant values dominated by the vortex dynamics at distances R—r ~ £. 

To summarize, we should say that Eq. (29) and its solution (30) cannot 
be taken literally as a rigorous result. They give a qualitative description of 
the kinetic energy profiles. Indeed, the predicted profiles are in qualitative 
agreement with numerical simulation calculations, as seen eg. in Fig. 23: Here 
K(r) goes to zero when r —>■ 0, increases at small r approximately as r 2 , as 
concluded in Eqs. (27) and (30a), reaches a maximum at r/R < 1 and then 
decreases in qualitative agreement with Eq. (30b). Similar comparison with 
other numerical results in Sec. 3.4 makes it believable that the analysis of 
the global energy balance in Sec. 3.5.4 and the predicted plateau for t> f in the 
T — > limit are reasonable. 

Nevertheless, both in the region of lower and higher temperatures the ex- 
periment shows deviation from this "plateau" (Fig. 26). The reason for this 
deviation at T > 0.35 T c is that turbulence is not well developed near the 
cylinder axis where the shear of the mean velocity, responsible for the turbu- 
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lent excitation, decreases. Therefore in the intermediate temperature region 
only part of the front volume is turbulent. When the temperature decreases, 
the turbulence expands toward the axis. We should therefore account also for 
the laminar contribution to the front velocity (12) near the axis. Notice fur- 
thermore that the simple sum of laminar and turbulent contributions to Vf 
in Eq. (26b) oversimplifies the situation, by not accounting for the turbulent 
motions in Eq. (26c). We will not discuss this issue, but just suggest an in- 
terpolating formula between the laminar and turbulent regimes, which has a 
shorter intermediate region than Eq. (26b): 



v f = V^lam + ^turb* ( 31 ) 

where Vf t i am and Vf t tur b are given by Eqs. (12) and (28). This interpolation is 
shown in Fig. 26 as a thin line for T < 0.3 T c and as a thick line for T > 0.3 T c . 
The agreement with measurements above 0.3 T c is good, but there is a clear 
deviation below 0.25 T c , where a < 10~ 2 . Below an explanation is outlined 
which takes into account the quantum character of turbulence, since as shown 



by Volovik (2003), individual vortex lines become important below 0.3 T c in 
the region where the transition to the lower plateau occurs in Fig. 26. 

To appreciate the last aspect note that the mean free path of 3 He quasipar- 
ticles at the conditions of the measurements at 29 bar pressure and T ~ 0.3 T c 
is close to £, while at 0.2 T c it exceeds R. This change from the hydrodynamic 
to the ballistic regime in the normal component may influence the mutual fric- 
tion force acting on individual vortices. However, what is important for the 
interpretation of the measurements is that the effect of the normal component 
on the superfluid component becomes negligibly small. Therefore it actually 
does not matter what is the physical mechanism of this interaction: ballis- 
tic propagation of thermal excitation with scattering on the wall or mutual 
friction that can only be described in the continuous-media approximation, 
as given by Eq. (14). For simplicity we use the hydrodynamic approxima- 
tion (14) in Sec. 3.5.6 to describe the turbulent kinetic energy dissipation at 
low temperatures. 

3.5.6. Mutual friction and bottleneck crossover from classical to quantum 
cascade 

a. Bottleneck at zero temperature. At low temperatures, when mu- 
tual friction becomes sufficiently small, the energy flux toward small length 
scales (or large wave vectors k) propagates down to the quantum scale 
I and vortex discreteness and quantization effects become most important. 
Even though some part of the energy is lost in intermittent vortex recon- 
nections, the dominant part proceeds to cascade below the scale t by means 
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by both slow vortex filament motions and fast vortex reconnection events. 
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As shown by L'vov et al. (2007a), the important point for the rate of energy 



dissipation (and consequently for the turbulent front velocity) is that Kelvin 
waves are much less efficient in the down-scale energy transfer than classical 
hydrodynamic turbulence: in order to provide the same energy flux as in the 
hydrodynamic regime, the energy density of Kelvin waves at the crossover 
scale i has to be A 10 / 3 times larger than that of hydrodynamic motions. For 
3 He-B, A = ln(£/a ) ~ 10, where a is the vortex-core radius. Assuming that 
the energy spectrum is continuous at the crossover scale £ and that no other 
mechanisms intervene between the classical and Kelvin-wave cascades, then 
to maintain the same value of energy flux there must be a bottleneck pile-up 
of the classical spectrum near this scale by the factor A 10 / 3 . To account for 
this phenomenon of energy pileup, we construct the "warm cascade" solution 
which will be described in what follows. 

For a qualitative description, we first define the hydrodynamic kinetic en- 
ergy density (in the one-dimensional /c-space) £ k , related to the total kinetic 
energy S as follows: 

£ = f S k dk. (32a) 



The Kolmogorov-1941 cascade of hydrodynamic turbulence (with fc-independent 
energy flux, e k =>- e=const) is described with the spectrum: 



S. 



£ 



2/3^-5/3 



(32b) 



The energy flux carried by the classical hydrodynamic turbulence with the 
K41 spectrum (32b) cannot adequately propagate across the crossover region 
at i. Therefore, hydrodynamic motions on larger length scales (smaller wave- 
vectors) will have increased energy content up to the level Sye, as required 
when the same energy flux has to be maintained by means of Kelvin waves. 
As a result, for k < l/£ the spectrum of hydrodynamic turbulence S k will not 
have the K41 scale-invariant form S k given by Eq. (32b). To get a qualitative 
understanding of the resulting bottleneck effect, we use the so called "warm 
cascade" solutions found by Connaughton and Nazarenko (2004). These so- 



lutions follow from the Leith (1967) differential model for the energy flux of 



hydrodynamical turbulence, 



k\ 13 F k 



dk 



TP — k 



(33) 



where F k is the 3-dimensional spectrum of turbulence. The generic spectrum 
with a constant energy flux can be found as the solution to the equation 

e k = e: 

24e ( T \ 3 / 2 l 



V 

7rp/ 



2/3 



(34) 



_ 1 1 |/C| H/2 

Here the range of large k values belongs to the thermalized part of the spec- 
trum, with equipartition of energy characterized by an effective tempera- 
ture T, namely T/2 of energy per degree of freedom, thus, F k = T up and 




10-2 



io-- 



io- 



kf 



10 



I0 2 



Fig. 27. Energy spectra in the classical, k < 1/1, and quantum, k > 1/1, ranges 
of scales. The two straight lines in the classical range indicate the pure K41 scaling 
£^ 41 oc AT 5 / 3 of Eq. (32b) and the pure thermodynamic scaling £f. oc k 2 . In the 
quantum range, the solid line indicates the Kelvin wave cascade spectrum oc fc~ 7 / 5 , 
whereas the dash-dotted line marks the spectrum corresponding to the non-cascad- 
ing part of the vortex tangle energy oc k^ 1 . 



£k = Tk 2 J rep. At low k, Eq. (34) coincides with the K41 spectrum (32b). 

This "warm cascade" solution describes reflection of the K41 cascade and 
stagnation of the spectrum near the bottleneck scale which, in our case, cor- 
responds to the classical-quantum crossover scale. To obtain the spectrum in 
the classical range of scales, it remains to find T, by matching Eq. (34) with 
the value of the Kelvin wave spectrum at the crossover scale Sk ~ n 2 /£. This 



gives T/p 



~ K 



.ii 



/A 5 e) 



1/4 



Obviously, the transition between the classical and quantum regimes is not 
sharp and in reality we should expect a gradual increase of the role of the 
self-induced wave-like motions of individual vortex lines with respect to the 
collective classical-eddy type of motions of vortex bundles. Thus, the high- 
wavenumber part of the thermalized range is likely to be wave rather than 
eddy dominated. However, the energy spectrum for this part should still be of 
the same k 2 form which corresponds to equipartition of thermal energy. This 
picture, as explained below, relies on the assumption that the self-induced wave 
motions have small amplitudes and, therefore, do not lead to reconnections. 

The resulting spectrum, including its classical, quantum and crossover parts, 
is shown in Fig. 27 as a log-log plot. It is important to note that in the 
quantum range k > l/£, in addition to the cascading energy associated with 
Kelvin waves, there is also energy associated with the tangle of vortex fila- 
ments (shown in Fig. 27 with the dash-dotted line). The energy spectrum of 
this part is ~ |&|~\ which is simply the spectrum associated with a singular 
distribution of vorticity along ID curves in 3D space ( |Araki et a/.|2002[ ). It does 
not support a down-scale cascade of energy. The cascading and non-cascading 
parts have similar energies at the crossover scale, that is the wave period and 
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Fig. 28. On the left the solutions of the differential Eq. (35) for different values 
of the dimensionless mutual friction parameter 7 = V / y k^,£ k + have been plot- 
ted. The smallest wave vector k + corresponds to the outer scale of turbulence: 
C ~ 2ir/k + . The wave vectors are normalized with respect to k + : x = k/k+ and 
F(x) = F x k + /Fk + . The dashed lines denote slopes of £k with K41 scaling oc /c -5 / 3 
and thermodynamic scaling £ k oc k 2 . On the right the slope F' = dF(x)/dx is shown 
at x = 1 for different values of mutual friction 7. The horizontal dashed line marks 
the K41 value of slope, 11/3. 



the amplitude are of the order of the characteristic time and size of evolv- 
ing background filaments. In other words, the scales of the waves and of the 
vortex "skeleton" are not separated enough to treat them as independent com- 
ponents. This justifies matching the classical spectrum at the crossover scale 
with the Kelvin wave part alone, ignoring the vortex "skeleton" . This is valid 
up to an order-one factor and justifies the way of connecting the "skeleton" 
spacing £ to the cascade rate e. 



b. Effect of mutual friction on the bottleneck crossover. Returning 
back to the propagating turbulent vortex front, recall that in the measure- 
ments of Sec. 3.3 the inertial interval R/£ is about one decade. Therefore the 
distortion of the energy spectrum owing to the bottleneck reaches the outer 
scale of turbulence. This leads to an essential suppression of the energy flux 
at any given turbulent energy or, in other words, to a decrease in the effec- 
tive parameter b, which relates e and K in the estimate (23c). The effect is 
more pronounced at low temperatures when mutual friction is small; thus b(T) 
should decrease with temperature. We analyze this effect with the help of the 
stationary energy balance equation for the energy spectrum Sk in fc-space 

£ = -iw&.*--(i-<0^^, (35) 

where T(T) is the temperature dependent damping in Eq. (16), and the energy 



flux over scales e(k) is taken in the Leith (1967) differential approximation. 



In addition to Eq. (33), we included in Eq. (35) the mutual friction correction 
factor (1 — a') and substituted F k = £k/k 2 . 
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Figure 28 displays the set of solutions for Eq. (35). We use C/£ = 12 as the 
ratio of the outer and crossover scales and characterize the bottleneck with 
the boundary condition £ k /[k 3 d(£ k /k 2 )/dk] = —4 • 10 5 at the crossover scale. 
One goal of these calculations is to find the slope of the function F k at the 
beginning of the inertial interval k = k + = 2n/C which, according to Eq. (35), 
characterizes the rate of energy input into the system at fixed value of the 
total energy £ , i.e. the phenomenological parameter b in the estimate (23c). 
The dashed line in the right panel shows the Kolmogorov value of this slope, 
— 11/3, which is associated with the classical value of b = b c \ xs 0.27. Thus 
the ratio 36 c i|F'|/ll can be interpreted as the effective value of b for a given 
value of damping T. To relate T with the temperature dependent value az(T), 
we substitute the r independent u e f[ from Eq. (27) to Eq. (16). This gives 
T = a(a/d)Q. The resulting function b(T), shown in the inset of Fig. 26 
(for VI = lrad/s), decreases from its classical value b c \ 0.27 to ~ 0.1 at 
T < 0.2 T c . Now, after accounting for the temperature dependence of b(T) in 
Eq. (28), we get the temperature dependence for the propagation velocity V{ 
of the quantum-turbulent front shown in Fig. 26 by the bold-solid line below 
0.3 T c . This fit is in good agreement with the measured data. We thus have to 
conclude that in this particular measurement the rapid drop in the dissipation 
rate on entering the quantum regime can be explained as a consequence of the 
relatively close proximity of the outer and quantum crossover scales in this 
measuring setup. 

3. 6. Summary: turbulent vortex front propagation 

In 3 He-B strong turbulence is restricted to the lowest temperatures below 
0.4 T c . Depending on the type of flow, turbulence varies in form and losses. The 
usual reference point is an isotropic and homogeneous turbulent vortex tangle. 
Here we have discussed an opposite extreme, the conversion of metastable ro- 
tating vortex-free counterflow to the stable equilibrium vortex state in a long 
circular column. The task is to explain how increasing turbulence with de- 
creasing temperature influences this type of polarized vortex motion in steady 
state propagation. The deterministic part of the motion takes place in the form 
of a spiraling vortex front followed by a helically twisted vortex cluster. When 
mutual friction decreases, turbulent losses start to contribute to dissipation, 
concentrating in the propagating front and immediately behind it. 

Experimentally the conversion of a metastable vortex-free state to an equi- 
librium array of rectilinear vortex lines can be arranged to occur at constant 
externally adjusted conditions. This is done with different types of injection 
techniques, which trigger the formation of a propagating vortex front and the 
trailing twisted cluster behind it. In a long rotating column the propagation 
can be studied in steady state conditions. At temperatures above 0.4 T c the 
propagation is laminar, but below 0.4 T c a growing influence from turbulence 
appears. This is concluded from the measured propagation velocity which 
provides a measure of the dissipation in vortex motion. In addition to the in- 
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creasing turbulent dissipation below 0.4 T c , the measurement shows a peculiar 
transition at ~ 0.25 T c between two plateaus, and a temperature independent 
finite value of dissipation on approaching the T — > limit. 

So far measurements on vortex propagation exist only for a column with 
circular cross section. This is an exceptional case of high stability (if the cylin- 
drical symmetry axis and the rotation axis are sufficiently well aligned). Here 
the crossover from the laminar to the turbulent regime is smooth as a func- 
tion of mutual friction dissipation. The same smooth behavior is also seen in 
the response to spin down, after a step-like stop of rotation. This is in stark 
contrast to the spin down of a column with square cross section, as we will see 
in Sec. 4. Similar differences seem to apply in viscous hydrodynamics where 
flow in a circular pipe is believed to be asymptotically linearly stable for all 



Reynolds numbers, in contrast to a pipe with square cross section (Peixinho 



and Mullin|2006 ) . In these measurements the formation and decay of turbulent 



plugs is monitored. Vortex plugs and fronts have also been observed in pipe 
flow of superfluid 4 He through long circular capillaries (Marees et al. 1987). 

Numerical vortex dynamics calculations have been used to examine the tur- 
bulent contributions to dissipation in vortex front propagation. This analysis 
demonstrates the increasing role of turbulent excitations on sub-intervortex 
scales with decreasing temperature. Analytical arguments have been devel- 
oped which explain that the transition to the lower plateau in the measured 
overall dissipation is caused by the difficulty to bridge the energy cascade from 
the quasi-classical to the quantum regime as a function of decreasing mutual 
friction dissipation. This bottleneck-scenario of L'vov et al. (2007a) in energy 
transfer was thus directly inspired by the experimental result. It is expected 
to apply foremost to polarized vortex motion in the rotating column where 
vortex reconnections are suppressed. This situation differs from the case of 
an ideal homogeneous vortex tangle, where Kozik and Svistunov ( 2008a| b) 
suggest that a bottleneck is avoided owing to the high reconnection rate. 

As for the second leveling off at the lowest temperatures, at present time 



three different types of measurements Bradley et al. (2006), Eltsov et al. 



(2007), Walmsley et al. (2007a)] conclude that dissipation in vortex motion 
remains finite at the lowest temperatures. This applies for both the fermion 
superfluid 3 He-B and the boson case of 4 He. Although the mechanisms for 
dissipation in the T — ► limit in these two superfluids are different and still 
under discussion, the results demonstrate that coherent quantum systems can 
be inherently lossy even in the T = state. 
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4. Decay of homogeneous turbulence in superfiuid 4 He 

4-1. Introduction and experimental details 



4-1.1. Quasiclassical and ultraquantum types of superfiuid turbulence 

In this section we switch to nearly homogeneous and isotropic superfiuid 
turbulence with an emphasis on the recent experiments in Manchester. A truly 
homogeneous turbulence is of course impossible to achieve in real experiments 
because the tangle is always confined by the container walls and can have 
other inhomogeneities specific to the particular process of its generation. 

The modern paradigm of homogeneous isotropic turbulence at high Reynolds 
numbers Re in classical liquids is that there is a broad range of wavenumbers 
k within which the energy, pumped at large scales, gets continuously redis- 
tributed without loss towards smaller length scales thus arranging a steady- 
state distribution in /c-space of the Kolmogorov type (Kolmogorov 1941a[b 
Batchelor|[l953| |Frisch||1995| ) (Eq. (32b)) 



Ce 2 / 3 \k\ 



-5/3 



(36) 



where the Kolmogorov constant was found to be C ~ 1.5 (Sreenivasan 1995[ ). 
In the steady state, the energy flux e is equal to the dissipation of the kinetic 
energy through the shear strain in the flow at short length scales; its rate 
integrates to 

£ = -v c \U 2 ch (37) 

where v c \ is the kinematic viscosity and u c \ is the r.m.s. vorticity. In realistic 
systems, the distribution Eq. 36 is truncated at small k = k\ by either the 
length scale of forcing or the container size h, and at large k = k 2 by the 
Kolmogorov dissipation scale A c i(i?e) ~ k 2 l which decreases with increasing 
Re = (k 2 /ki)^ 3 . In superfluids, the inertial cascade is expected to operate if 



the mutual friction parameter a is sufficiently small (Vinen 2000). At large 
length scales > £ the cascade is classical. However, at short length scales < £ 
the cascade becomes quantum as the discreteness of the vorticity in superfluids 
adds new behaviour. At sufficiently low temperatures, when the energy reaches 
very small scales <C £ without dissipation, it is expected that the quantum 



cascade takes the form of a non-linear cascade of Kelvin waves (Svistunov 



1995, Kozik and Svistunov 2004), eventually being truncated at some quantum 
dissipation scale A q (a). While the theory of this regime at very high k is 
now established, no direct observations exist so far. The most complicated 
is, of course, the transitional region between these two clear-cut limits, the 
Kolmogorov and Kelvin-wave cascades. The question being currently debated 
is whether the energy stagnates at k < £~ 1 due to the poor matching in 



the kinetic times of the two cascades (L'vov et al. 2007a), or gets efficiently 



converted from 3-d classical eddies to 1-d waves along quantized vortex lines 



with the help of various reconnection processes (Kozik and Svistunov 2008a). 



When comparing the two superfiuid isotopes from the experimental point 
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of view, the T = limit in superfluid 4 He is achievable at some T < 0.5 K 
( Vinen||2000 Walmsley et a/.|2007a ) while for 3 He-B much lower temperatures 
T < 0.5 mK are required. Owing to the two or three orders of magnitude 
smaller core size in 4 He, a$ ~ 0.1 nm, the Kelvin- wave-cascade is expected to 
extend over a broader range of length scales, down to ~ 3nm (Vinen 2001 



Kozik and Svistunov 2008b). It has been suggested (Vinen 2002, Vinen and 
Niemela 2002) that in 3 He-B Kelvin waves should become overdamped at 
frequencies ~ 10 kHz (corresponding to wavelengths ~ 2 /im) due to the 
resonant scattering on core-bound quasiparticles ( Kopnin and Salomaa|1991 )). 

Flow on a scale greater than the inter-vortex distance £, which is initially 
typically between ~ 10 and 100 /im in experiments, can be obtained by me- 
chanical stirring of the liquid. So far the following methods have been used: 



counter-rotating agitators with blades (Maurer and Tabeling 1998), pipe flow 



(Smith et al. 


1999 


Roche 


2007 


), flow through an orifice (Guenin and Hess 


1978 


), towed 


grids ( 


Smith et al. 


1993 


Stalp et al. 


1999, 


Niemela et al. 


2005). 



Hanninen et a/.|2007b Vinen and Skrbek 2008 ), as well as wires ( Bradley|2000 
Nichol et al. 2004), and microspheres (Schoepe 2004), plus most recently im- 



pulsive spin-down to rest of a rotating container ( Walmsley et a/.|[2007a ). It is 
also possible to initiate large-scale flow in superfluid 4 He without any moving 
parts or rotation of the cryostat: by either running thermal counterflow in 
wide channels at T > 1 K (Barenghi and Skrbek 2007, Chagovets et al. 2007) 



or a jet of injected current (Walmsley and Golov 



2008a). 



Actually, the turbulence can take two very different forms depending on 
whether the forcing is at scales above or below I. So far we were discussing 
the flow on classical scales > £, where the energy cascades towards shorter 
length scales like in the Richardson cascade in classical turbulence; the large 
quasiclassical eddies being the result of correlations in polarization of vortex 
lines. On the other hand, when forced on quantum scales < £, the resulting 
uncorrelated tangle has no classical analogs and should have completely differ- 
ent dynamics first described by Vinen ( 1957a|b~flc 1958) and later investigated 
numerically by Schwarz (1988). In both cases, the dissipation of flow energy 



is through the motion of vortex lines; its rate per unit mass can be assumed 



Stalp et al. 


1999 


Vinen 


2000, 


Stalp et al. 


2002) 



£ 



-v(kL)' 



(3* 



where k 2 L 2 is an effective total mean square vorticity and the parameter v is 
believed according to various models to be approximately constant for a given 
temperature and type of flow. This formula is the quantum analog of the clas- 
sical expression Eq. (37). As we will show below, the efficiency of the vortices 
in dissipating energy in the T = limit, expressed through the "effective 
kinematic viscosity" is, is different for these two regimes of turbulence. 

To measure the values of is, one can monitor the free decay of homogeneous 
tangles of both types. In any tangle, the quantum energy associated with the 
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quantized flow on length scales r < £ is £ q = jL/ p s (per unit mass), where 
the energy of vortex line per unit length is 7 = Bp s n 2 and B m ln(£/ao)/47r is 
approximately constant. This is the same energy that is shown in the quantum 
range k > l/£ (including both non-cascading vortex "skeleton" and Kelvin- 
waves) in Fig. 27. If the total energy is mainly determined by S q , from Eq. (38) 
we arrive at the late-time free decay 

L = Bv~ l r l . (39) 

This type of turbulence, without any motion on classical scales, will be called 



ultraquantum (in the literature one can find e.g. "Vinen turbulence" (Volovik 



2003, 2004) and "random" or "unstructured" vortex tangle). For 4 He, K4 ■ 
2irh/m^ a>o ~ 0.1 nm and B ~ 1.2, while for 3 He-B, K3 = 2 r nhj{2m^) 1 ao • 
13 — 65 nm and B rs 0.7. 

Now suppose that the tangle is not random but structured due to the pres- 
ence of flow on classical length scales r > £, and the additional energy of this 
classical flow S c is much greater than S q . This type of turbulence will be called 
quasiclassical (e.g. "Kolmogorov turbulence" and "structured" or "polarized" 
vortex tangle). For the Kolmogorov spectrum between wavenumbers k\ and 
k 2 (ki -C k 2 ), while the size of the energy-containing eddy stays equal to the 



size of container h (i.e. k\ ~ 2ir/h), the late-time free decay follows (Stalp et 



al. 1999) 



L = {ZCf^K-^u- 1 ' 2 ^' 2 . (40) 

Eventually, when the energy flux from the decay of classical energy S c will 
become smaller than that from the quantum energy S q (which should typically 
happen when only a couple of vortex lines are left in the container, i.e. L ~ 
h~ 2 ), the quasiclassical regime will cross over to the ultraquantum one, so 
L oc t~ 3 / 2 decay will be replaced by L oc f -1 . Deviations from the described 
scenario are also possible if some fraction of the classical energy is stored in 



non-cascading "thermal spectrum" (L'vov et al. 2007a). 



In this Section we first review the relevant experimental techniques of de- 
tecting and generating turbulence. Next we describe the more recent new 
developments at the very lowest temperatures. The results of these measure- 
ments will be discussed in detail in the last part of this section, where we 
conclude with a general discussion of the significance of the observed temper- 
ature dependence of v. 



4-1.2. Techniques of measuring the vortex line density L 



The attenuation of second sound in superfluid 4 He, which is proportional to 
the total length of vortex lines, has been the prefered technique of measuring 



L in the temperature range 1 - 2 K f 


Vinen 


1957a 


Stalp et al. 


1999 


Niemela 


et al. 


2005 


Barenghi and Skrbek 


2007 


') since 


; Hall and Vinen 
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1956) 


. however 



it cannot be extended below 1 K. Very recently, scattering of negative ions 
off vortex lines has been succesfully utilized at temperatures 80 mK 1.6 K 



, showing good agreement with previous measurements 



(Walmsley et al.|2007a 
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of v at overlapping temperatures (Stalp et al. 2002). Actually, the first use of 
ions for the observations of a vortex tangle was reported almost fifty years ago 
( Careri et a/.|1960 ), and numerous further experiments revealed their potential 
for investigating turbulence in superfluid 4 He (see the reviews by Tough ( 1982 ) 
and Donnelly (1991) for more references on early studies). 

It is worth mentioning possible alternative techniques of detecting turbu- 
lence in 4 He in the T = limit that are currently being developed, such as 
calorimetric (Ihas et al. 2008) and others (Vinen 2006). The design of towed 
grid for turbulence generation, suitable for experiments at T < 0.5 K in 4 He is 
presently being attempted by McClintock at Lancaster and Ihas at Gainsville. 



Optical vizualization of individual vortex lines in the turbulent tangle (Bewley 
e/ al.|2006l|2008l|Zhang et a/.|[2005| [Rellergert et al |2008[ ISergeev et a/.||2006| 



(see the contribution by Van Sciver and Barenghij ([2008 ) in this volume) is 
another promising development, although the prospects of extending this tech- 
nique down to at least T = 0.5 K seem distant. As we mentioned in Sec. 2 and 
3, computer simulations have contributed a great deal to our understanding 
of certain processes in the dynamics of tangles and of the Kelvin-wave cascade 
at T = in particular (iSchwarz 11988 , Kivotides et al. 2001, Tsubota et al. 



2000 , Vinen et al. 2003 , |Hanninen et al. 2007b ) (see also the contribution by 



Tsubota and Kobayashi in this volume); however the proper modelling of the 
inertial cascade pumped at large (quasiclassical) length scales and dissipated 
at short (quantum) scales requires a considerable range of length scales and 
hence remains a challenge. 



4-1.3. General properties of injected ions in liquid helium 



Schwarz 1975 


Fetter 1976 


Donnelly 1991 


Borgh- 



esam 



2007[ )) are convenient tools to study elementary excitations and quan- 



tized vortices in superfluid helium. They were used to detect vortex arrays 



in rotating helium ( Careri et al. 1962 , Yarmchuk et al. 



1979), quantized vor- 



tex rings (Rayfield and Reif 1964), and vortex tangles (Careri et al. 1960) in 
4 He. To create a negative ion, an excess electron is injected into liquid he- 
lium, where it self-localizes in a spherical cavity ( "electron bubble" ) of radius 
~ 2 nm at zero pressure. To create a positive ion, an electron is removed from 
one atom, which results in a positively-charged cluster ion ( "snowball" ) of ra- 
dius ~ 0.7nm. These objects are attracted to the core of a quantized vortex, 
resulting in trapping with negligible escape rate provided the temperature is 
low enough (T < 1.7 K for negatives and T < 0.4 K for positives in 4 He 
at zero pressure). Because of this interaction, both species have been used a 
great deal to study quantized vortices although negative ions are often more 
convenient, especially if one wants to cover a wider range of temperatures. The 



binding energy of a negative ion trapped by a vortex is about 50 K (Donnelly 
1991). In superfluid 3 He, thanks to much larger diameters of vortex cores this 



binding is much weaker, hence there were no observations of ion trapping on 
vortex cores yet. Still, the anisotropy of ions motion in 3 He-A made it possible 
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to use negative ions to detect textural changes around vortex cores ( Simola et 



al. 1986). 



In this section, we focus on negative ions in superfluid 4 He. Their mobility 
is limited by scattering thermal excitations and is hence rapidly increasing 
with cooling. Below some 0.8 K it is so high that they quickly reach the 
critical velocity for creation of (depending on pressure and concentration of 
3 He impurities) either rotons (and then move with the terminal Landau ve- 
locity ~ 60 m/s, shedding off rotons continuously) or quantized vortex rings 
(and then get trapped by such a ring and move with it as a new entity - a 
charged vortex ring) ( |Hendry et al. 1990 1988 Winiecki et al. 2000, Berloff 
and Roberts 2000a[b ). Above 1 K, the trapping diameter a of an ion by a 
vortex line is inversely proportional to the electric field, a oc E~ x . In a field 
E = 33 V/cm it is decreasing with temperature from a = 100 nm at T = 1.6 K 
to a = 4nm at T = 0.8 K (Ostermeier and Glaberson 1974). 



Rayfield and Reif (1964) produced singly-charged singly-quantized vortex 



rings in 4 He at T = 0.4 K and, by measuring the dependence of their self- 
induced velocity v oc k 3 K~ x , Eqs. (41) - (42), on their energy K, confirmed 
the value of the circulation quantum k = 1.00 x 10~ 3 cm 2 s _1 . Schwarz and 



Donnelly ( 1966 ) investigated the trapping of charged vortex rings by rectilin- 



ear vortices in a rotating cryostat and found the trapping diameter to be of 
order of the ring radius, meaning that the interaction is basically geometrical. 
They wrote: "quantized vortex rings are very sensitive "vortex-line detectors," 
making them suitable probes for a number of problems in quantum hydrody- 



namics 



n 



Guenin and Hess (1978) used charged vortex rings to detect turbu- 
lent vortex tangles in 4 He at T = 0.4 K, created by forcing a jet of superfluid 
through an orifice at supercritical velocities. 

While charged vortex rings are more convenient than free ions because of 
their much greater trapping diameter, their dynamics is more peculiar (Fig. 29; 



all figures in Sections 4.1-4.3 represent the Manchester experiments (Walmsley 



et al. 2007a, Walmsley and Golov 2008a)). The energy K, velocity v and 
impulse P of a quantized vortex ring of radius 1Z with a hollow core of radius 



ao and no potential energy in the core are (Glaberson and Donnelly 1986 



Donnelly 1991) 



K 



-p s K 2 TZ 



K 



hi 



, 8TZ 

m — 

a 

8K 



1 



(41) 



(42) 



AirTZ v a 
P = np s KTZ 2 . (43) 
By integrating these equations, one can calculate the trajectories of the charged 



vortex rings subject to a particular electric field (Walmsley et al. 2007b) 



4- 1-4- Measuring the vortex line density L by ion scattering 

In the absence of vortex lines the ions would propagate either freely (these 
dominate at T > 0.8 K) or riding on small vortex rings (dominate at T < 
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Fig. 29. Time of flight (corresponding to the leading edge of the pulse of collector 
current shown in inset) of charged vortex rings as a function of electric field at 
T = 0.15 K in both the horizontal and vertical direction. The solid line show the 
times of flight for rings with initial (as injected into the drift volume) radii and 
energies of 0.53 /jm and 21 eV, calculated using Eqs. (41), (42), and (43). 

0.6 K). Without other vortex lines in their way, the time of flight of both 
species over a particular distance is a well-defined function of temperature, 
electric field, and initial radius of the attached vortex ring (if any). Hence, 
after injecting a short pulse of such ions, a sharp pulse of current arrives at 
the collector as shown in Fig. 30. The interaction with other existing vortices, 
that happened to be in the way of the propagating ions, is characterized by a 



trapping diameter" a ( Schwarz and Donnelly 1966 , Ostermeier and Glaberson 



1974) and leads to the depletion of the pulses of the collector current. This 
is used to measure the average density of vortex lines L between the injector 
and collector through 



L{t) = (a/i)- 1 ln(/(oo)//(t)). 



(44) 



To calibrate the value of a for either free negative ions or ions trapped on 
a vortex ring at different temperature and electric field, one can measure 
the attenuation of the pulses of the collector current when the cryostat is at 
continuous rotation at angular velocity Q (thus having an equilibrium density 
of rectilinear vortex lines L = 2Q/k) as in figure 31. This is best done in the 
direction perpendicular to the rotation axis. 

At low temperatures in quantum cascade Kelvin waves of a broad range of 
wavelengths are excited, however the main contribution to the total length L 
converges quickly at scales just below t ( Kozik and Svistunov|2008b| ). Hence, a 
probe ion with the trapping diameter a <C £, moving at a speed (v ~ 10 cm/s) 
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Fig. 30. Examples of current transients produced by short pulses for three different 
temperatures (the temperatures, pulse durations and mean driving fields are indi- 
cated in the legend). The solid lines show the transients without a vortex tangle 
in the ions' path, while the dashed ones represent the transients suppressed by the 
vortex tangle which has been decaying a specified time (indicated near curves) af- 
ter stopping generation. The charge carriers are either free ions (T = 1.60 K and 
T = 0.90 K) or charged vortex rings (T = 0.15 K). The electronics time constant is 
0.15 s, hence the time of arrival of the fastest peak (T = 0.90 K) cannot be resolved. 



much greater than the characteristic velocities of vortex segments (~ k/£ < 
10 _1 cm/s), should sample the full length L. Experimental measurements with 
different o = 0.4-1.7 fim (using charged vortex rings in a range of driving 
electric fields) indeed produce consistent values of L. 

There are evidences that small concentrations of trapped space charge 
(when the ratio of trapped space charge density n to the vortex line length L, 



does not exceed n/L ~ 10 cm ) do not affect the tangle dynamics (Walmsley 



and Golo"v||2008b ) . For example, experiments generating quasiclassical tangles 



at T > 0.7 K by a jet of ions (and hence resulting in a substantial trapped 
space charge) revealed that the late-time decay is identical to that for the tan- 
gles generated by an impulsive spin-down to rest without injecting any ions. 
Still, one might worry that the very presence of trapped charge, as it interacts 
with the beam of probe ions, might affect the result of measuring L. Hence, 
in order not to introduce extra turbulence and not to contaminate the tangle 
with any ions, the measurement was performed by probing each realization of 
turbulence only once - after a particular waiting time t during its free decay. 
Then the contaminated tangle was discarded and a fresh tangle generated. 
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Fig. 31. Dependence of pulse amplitude on the angular velocity of rotation, Q, in 
examples of calibration measurements. The temperatures, driving fields and trap- 
ping diameters a are indicated. The charge carriers are either free ions (T = 1.60 K) 
or charged vortex rings (T = 0.78 K and T = 0.50 K). 

4-1.5. Design of ion experiment 

The cell used in the Manchester experiments had cubic geometry with sides 
of length h = 4.5 cm. A schematic drawing of the cell is shown in Fig. 32. The 
relatively large size of the cell was important to enhance the efficiency of ion 
trapping and the time resolution of vortex dynamics, and was really instru- 
mental to ensure that the continuum limit £ <C h holds even when the vortex 
line density drops to as low as just L ~ 10 cm -2 . This also helped ensure 
that the presence of the walls, which might accelerate the decay of turbulence 
within some distance ~ £, does not affect the dynamics of the turbulent tangle 
in the bulk of the cell. In order to probe the vortex densities along the axial 
and transverse directions, there were two independent pairs of injectors and 
collectors of electrons. Both injectors and collectors were protected by electro- 
static grids, enabling injecting and detecting pulses of electrons. The injectors 



were field emission tips made of 0.1 mm diameter tungsten wire (Golov and 
Ishimoto]|1998 ). The threshold for ion emission was initially ~ -100 V and -210 



V for the bottom and left injectors respectively; however, after over two years 
of almost daily operation they changed to some -270 V and -520 V. The fact 
that the two injectors had very different threshold voltages helped investigate 
the dependence (or rather lack of it) of the radius of initial charged vortex 
rings on the injector voltage. The six side plates (electrodes) that make up 
the cube can be labeled "top" , "bottom" , "left" , "right" , and two "side" elec- 
trodes. The top, bottom, left and right electrodes had circular grids in their 
centers. All grids were made of square tungsten mesh with period 0.5 mm and 
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Fig. 32. A side cross-section of the Manchester cell. The distance between opposite 
electrode plates is 4.5 cm. An example of the driving electric field is shown by dashed 
lines. Such a configuration was used in all the Manchester experiments described 
in this section. It was calculated for the following potentials relative to the right 
electrode: left electrode at -90 V, side, top and bottom electrodes at -45 V, and 
the right electrode and collector at producing a 20 V/cm average driving field 
in the horizontal direction. To inject ions, the left tip was usually kept at between 
-500 to -350 V relative to the left grid. When ions traveling vertically across the 
cell (injected from bottom tip and detected at the top collector) were required, the 
potentials on the electrodes were rearranged as appropriate. 

wire diameter 0.020 mm, giving a geometrical transparency of 92%. The grids 
in the bottom, left and right plates had diameter 10 mm and were electrically 
connected to those plates. The grid in the top plate had diameter of 13 mm 
and was isolated from the top plate. The injector tips were positioned about 
1.5-2 mm behind their grids. The two collector plates were placed 2.5 mm 
behind their grids and were typically biased at +10 — 1-25 V relative to the 



grids. Further details can be found in Walmsley et al. (2007b) 



4-1.6. Tangles generated by an impulsive spin-down 

This novel technique of generating quasiclassical turbulence, suitable for 
any temperatures down to at least 80 mK, relied on rapidly bringing a rotat- 



ing cubic-shaped container of superfmid He to rest (Walmsley et al. 2007a). 



The range of angular velocities of initial rotation Q was 0.05-1.5 rad/s. In clas- 
sical liquids at high Re, spin-down to rest is always unstable, especially at high 
deceleration and in axially asymmetric geometries. It is known that, within a 
few radians of initial rotation upon an impulsive spin-down to rest, a nearly 
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Fig. 33. Cartoon of the vortex configurations produced by spin-down in the experi- 
mental cell (side view) at different stages. (1) Regular array of vortex lines during 
rotation at constant Q before deceleration. (2) Immediately after stopping rotation 
(0 < < 10), turbulence appears at the outer edges but not on the axis of rotation. 
(3) After about 30 rad of initial rotation (fit ~ 30), 3D homogeneous turbulence is 
everywhere. (4) Then (Ot ~ 10 3 ) the 3D turbulence decays with time. Shaded areas 
indicate the paths of probe ions when sampling the vortex density in the transverse 
(L t , 3) and axial (L a , 4) directions. 



homogeneous turbulence develops with the energy-containing eddies of the 
size of the container. In a cylindrical container, centrifugal and Taylor- Gortler 
instabilities usually break up at the perimeter, thus facilitating slowing down 
of the outer region of the initially rotating liquid. Simultaneously, because of 
the Ekman pumping of non-rotating liquid into the central axis via the top 
and bottom walls, the central core of the initial giant vortex slows down too 
(Donnelly 1991). In a cubic container, the turbulence becomes homogeneous 
much faster. Some residual rotation of the central region in the original di- 
rection might still survive for a while but, as we show below, the generated 
turbulence is in general pretty homogeneous. 

This behaviour is well-documented for classical liquids (van Heijst 1989), 
and one expects similar processes to occur in a superfluid liquid providing the 
process of initial multiplication of vortices does not affect the dynamics. As 
spin-down experiments always begin from already existing dense rectilinear 
vortex arrays of equilibrium density L = and rapid randomization and 

multiplication of these vortices is expected due to the lack of axial symmetry of 
the container, as well as surface pinning (friction) in the boundary region, these 
seem to be sufficient for the superfluid to mimic the large classical turbulent 
eddies. In Fig. 33 we show the four different stages of the evolution from a 
vortex array to a decaying homogeneous tangle upon an impulsive spin-down 
to rest, which are in agreement with the observations outlined below. 

Before making each measurement, the cryostat was kept at steady rotation 
at the required Q for at least 300 s before decelerating to full stop, then 
waiting a time interval t and taking the data point. Then the probed tangle 
was discarded and a new one generated. Hence, different data points represent 
different realizations of the turbulence. The deceleration was linear in time 
taking 2.5 s for Q = 1.5 rad/s and 0.1 s for Q = 0.05 rad/s. The origin t = 
was chosen at the start of deceleration. 
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Fig. 34. Cartoon of the vortex configurations produced by a pulse of injected ions at 
T < 0.5 K in the experimental cell (side view) at different stages. (1, < 1 s) A pulse 
of charged vortex rings is injected from the left injector. While most make it to 
the collector as a sharp pulse, some got entangled near the injector. (2, ~ 5 s) The 
tangle spreads into the middle of the cell. (3, ~ 20 s) The tangle has occupied all 
volume; from now on it is nearly homogeneous (as probed in two directions). (4, up 
to 1000 s) The homogeneous tangle is decaying further. The shaded areas indicate 
the trajectories of ions used to probe the tangle along two orthogonal directions. 



4-1-7. Current-generated tangles 

An alternative technique of generating turbulence, by a jet of injected ions 
that does not require any moving parts in the cryostat, has been developed 
(Walmsley and Golov 2008a) (Fig. 34). In early experiments with injected 
ions at low temperatures (Bruschi et al. 1966, Bowley et al. 1982), it was 



observed that a pulse of negative ions through superfluid 4 He leaves behind a 
tangle of vortices. Walmsley et al. (2007a) found that the properties of these 
tangles can be quite different. As we explain below, the tangles generated 
after long injection at high temperatures possess the properties of developed 
quasiclassical turbulence while those produced after short injection at low 
temperatures have signatures of random tangles with little large-scale flow. 

At high temperatures T > 0.7 K, while moving relative to the normal 
component, the ions experience viscous drag through scattering of excitations 
- this means that the ions entrain the normal component along. They can also 
get trapped by an existing vortex line (provided T < 1.8 K), after which the 
vortex line will be pulled along with the current of such ions. The existing 
vortices and hence the superfluid component, through the action of mutual 
friction, will be pulled by the already entrained normal component too, and 
vice versa. Then one can expect the injected current to produce a large-scale 
jet-like flow of liquid helium which is an efficient way of driving quasiclassical 
turbulence at large scales. This is similar to jet flow through an orifice - one of 



the traditional means of generating turbulence in classical liquids (Chanal et 



al. 2000). On the other hand at low temperatures the ions always bring small 



vortex rings along, hence, upon formating a tangle, can directly pump up the 
vortex length L without introducing substantial large-scale flow - at least for 
not too-long injections. 



;4 



75 



10 4 



o 



1(T- 



10 



TTTTTT| 1 1 < — I I I I I | 1 1 1 ^^pl I II | ~~I 1 1 — I I I I I | 1 1 1 — I I I I 

, * 0.5 rad/s ^pjkg - '■, Q = 1 ,5 rad/s 



□ □□ ° 

□ ° D □ 

°n 0,15rad/s 



+ V + + + + + 

- 

0,05 rad/s 



-3/2 



+ + + ° o 





ton. 



T=0.15K 



10' 



10' 



10 

t, s 



io J 



10* 



Fig. 35. Lt(t) at T = 0.15 K for four values of Q. The average driving fields used 
for n = 1.5 rad/s: 5 V/cm (o), 10 V/cm (y), 20 V/cm (o), 25 V/cm (A). At other 
values of f2 the electric fields used were either 20 V/cm (0.05 rad/s) or 10 V/cm (0.5 
and 0.15 rad/s). The dashed line shows the dependence t~ 3 / 2 . The horizontal bars 
indicate the initial vortex densities at steady rotation, L = 2£l/n, at f2 = 1.5 rad/s, 
0.5 rad/s, 0.15 rad/s and 0.05 rad/s (from top to bottom). 



4-2. Experimental results 

4-2.1. Quasiclassical turbulence generated by spin-down to rest 

Above 1 K, the L oc t~ 3 ^ 2 free decay was monitored by second sound for 



towed grid turbulence in Oregon (Stalp et al. 1999). The recent spin-down 



experiments at Manchester are in good agreement with them in the overlap- 
ping temperature regime. Below 1 K, only scattering of ions off vortex lines 
has been used to measure L so far. 

In Fig. 35, the measured densities of vortex lines along the horizontal axis 
(transverse, L t ) are shown for four different initial angular velocities Q. During 
the transient, which lasts some ~ lOOQ^ 1 , L t (t) goes through the maximum 
after which it decays eventually reaching the universal late-time form of L oc 
£-3/2 _ p or ^ > q 5 rac i/ S; the values of L at maximum were too high to be 
detected. The initial vortex densities at steady rotation, L — 2Q/k, are shown 
by horizontal lines. Similarly, the densities of vortex lines measured along the 
vertical axis (axial, L a ) are shown in Fig. 36. 

To illustrate what is happening near the vertical axis of the container at 
different stages of the transient after a spin-down, in Fig. 37 we show five 
records of the current to the top collector arriving after a short (0.1 s) pulse 
of probed ions was fired from the bottom injector. Each is characteristic of 
a particular configuration of vortex lines near the rotational vertical axis of 
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Fig. 36. L a (t) at T = 0.15 K for four values of f2. The average driving field used was 
20 V/cm in all cases except at 0.5 rad/s where both 10 V/cm (v) and 20 V/cm 
(A) were used. The dashed line shows the dependence i -3 / 2 . The horizontal bars 
indicate the initial vortex densities at steady rotation, L = 20,/n, at = 1.5 rad/s, 
0.5 rad/s, 0.15 rad/s and 0.05 rad/s (from top to bottom). 

the container during the transformation from an array of parallel lines to a 
homogeneous decaying vortex tangle. One can see three different characteristic 
times (vertical dashed lines in Fig. 37) of arrival of ions via different means. 
The first peak at w 0.4 s (determined by the time constant of the current 
preamplifier, 0.15 s) corresponds to the ions trapped on the rectilinear vortex 
lines which can slide along those lines very quickly, provided those lines are 
continuous from the injector to collector as during steady rotation. The second 
peak at ~ 1 s corresponds to the coherent arrival of the ballistic charged vortex 
rings from the bottom to the top. The third broad peak at times ~ 3 s (but 
with a long tail detectable until ~ 40 s) corresponds to the charge trapped on 
the vortex tangle and drifting with the tangle. Hence, we have the following 
regimes (and curves in Fig. 37) labeled by the time t after the spin-down: 

(1) Qt < 0, steady rotation. The nearly equal first and second peaks tell that 
there exist a vortex array without disorder (otherwise the second peak 
would get suppressed and the third would appear). 

(2) Qt = 1.5. The first peak gets enhanced three-fold while the second one 
is still there (and no third peak) - meaning more trapped ions can now 
reach the collector along the rectilinear vortex lines while there is not 
much turbulence in this region yet. 

(3) Qt = 6. The first peak has disappeared in favour of the second one which 
has got broadened - at this stage the rectilinear vortices should have be- 
come scrambled in the Ekman layers near the top and bottom walls/grids 
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Fig. 37. Records of the current to the top collector, injected from the bottom injector 
as a 0.1 s-long pulse at the time t after an impulsive spin-down from = 1.5 rad/s 
to rest. T = 0.15 K, E = 20 V/cm. 

while the ballistic charged vortex rings are still the dominant transport 
of charge. 

(4) Qt = 12. Now both fast peaks have disappeared completely while the 
third broad peak has emerged - this means that a turbulent tangle has 
finally reached the axial region. 

(5) Qt — > oo. The sharp second peak has recovered but all others vanished 
- after the turbulence has decayed only ballistic charged vortex rings 
carry the charge, neither the rectilinear array nor the turbulent tangle 
contribute to the transport any more. 

In Fig. 38 the measured densities of vortex lines along the horizontal, L t , 
and vertical, L a , axes are shown, by solid and open symbols respectively. To 
stress the scaling of the characteristic times with the initial turn-over time Q^ 1 
and the universal late-time decay oc t~ 3 / 2 , the data for different Q are rescaled 
accordingly. We can see that at all fl the transients are basically universal. 
Immediately after deceleration, L t jumps to pa 10 4 f2 3 / 2 cm _2 s 3 / 2 , indicating 
the appearance of the turbulent boundary layer at the perimeter, while L a is 
stable at L; pa 2Q/k. Only at Qt pa 3, the latter starts to grow, signalling the 
destruction of the rotating core with vertical rectilinear vortices. After passing 
through a maximum at Qt = 8 and Qt = 15 respectively, L t and L a merge 
at Qt ~ 30 and then become indistinguishable. This implies that from now 
on the tangle density is distributed nearly homogeneously. The scaling of the 
transient times with the turn-over time Q" 1 indicates that transient flows are 
similar at different initial velocities Q, which is expected for flow instabilities 
in classical inviscid liquids. Eventually, after ilt ~ 100, the decay takes its late- 
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Fig. 38. J7~ 3 / 2 L t (t) (filled symbols) and J7~ 3 / 2 L a (t) (open symbols) vs. ttt for four 
values of 17 at T = 0.15 K. The straight line oc i~ 3 / 2 guides the eye through the 
late-time decay. 

time form L oc t~ 3//2 expected for quasiclassical isotropic turbulence, whose 
energy is mainly concentrated in the largest eddies bound by the container 
size h, but homogeneous on smaller length scales. We hence assume that the 
turbulence in 4 He at this stage is nearly homogeneous and isotropic, and can 
apply Eq. (40) to extract the effective kinematic viscosity v. The cross-over 
to the L oc t~ x regime at late time has never been observed, probably because 
the measured values of L never dropped below 10 cm -2 ^> h~ 2 ~ 0.05 cm -2 . 

At 0.08 < T < 0.5 K the measured L(t) were independent of temperature. 
In Fig. 39, we compare the transients L t (t) at low (T = 0.15 K) and high (T = 
1.6 K) temperatures. The prefactor in the late-time dependence L oc t -3 / 2 at 
T = 0.15 K is almost an order of magnitude larger that that for T = 1.6 K. 
This implies that at low temperatures, the steady-state inertial cascade with 
a saturated energy-containing length and constant energy flux down the range 
of length scales requires a much greater total vortex line density. This means 
that in the T = limit the effective kinematic viscosity v is approximately 
70 times smaller than at T = 1.6 K (provided C ~ 1.5 and k\ « 2n/h are 
independent of temperature). 

At all temperatures the transients L t (t) after a spin-down are, in first ap- 
proximation, universal, i.e. the timing of the maximum is the same fit ~ 7 and 
the amplitudes of the maximum are comparable. This supports our approach 
to turbulent superfTuid helium at large length scales as to an inviscid classi- 
cal liquid, agitated at large scale and carrying an inertial cascade down the 
lengthscales, independent of temperature. On the other hand, as the temper- 
ature increases and interactions with the viscous normal component become 
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Fig. 39. n- 3 / 2 L t (t) vs. nt for T = 0.15 K (filled symbols) and T = 1.6 K (open 
symbols). Dashed and solid lines oc i~ 3 / 2 guide the eye through the late-time decay 
at T = 1.6 K and 0.15 K, respectively. 



stronger, changes in the shape of transients might be expected. Indeed, one can 
see that the slope of L(t) after the maximum but before reaching the ultimate 
late-time decay L oc t~ 3 / 2 {i.e. for 10 < Qt < 100) is changing gradually with 
increasing temperature from being less steep than £~ 3 / 2 at T = 0.15 K to more 
steep than t~~ 3 / 2 at T = 1.6 K (Fig. 39). At temperatures around T = 0.85 K, 
it nearly matches t~ 3//2 , thus making an erroneous determination of v possible 
by taking this part of the transient for the late-time decay L oc t~ 3 / 2 . Indeed, 
in the first publication (Walmsley et al. 2007a), parts of some transient at 
T = 0.8 — 1.0 K for as early as Qt > 15 were occasionally used to be fitted by 
L oc i~ 3 / 2 that often resulted in overestimation of the value for v. To rectify 
this, we have re-fitted the datasets used in the original publication as well as 
subsequent measurements with L oc t~ 3 / 2 for spin-downs using the following 
rules: for T < 0.5 K only points for Qt > 300 were used, between 0.5 K and 
1.0 K only points for Qt > 150 used and at T > 1.0 K only Qt > 75 were used. 
This resulted in slight systematic reduction of the extracted values of v{T) 
at temperatures 0.8-1.2 K from those published in Walmsley et al. ( 2007a[ ) ; 
what looked as a rather steep drop in v{T) at 0.7-0.8 K, now occurs at 0.85 
- 0.90 K and is somewhat reduced in magnitude. 



4-2.2. Quasiclassical turbulence generated by an ion jet 

As mention ed above, ion jets can produ ce quasiclassical tangles decaying 
as L oc t~ 3//2 (Walmsley and Golov 2008a). Examples of the decay of such 
tangles alongside with those obtained by a spin-down at T = 1.6 K and 
T = 0.8 K are shown in Figs. 40 and 41, demonstrating good quantitative 
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Fig. 40. Free decay of a tangle produced by a jet of free ions from the bottom injector 



(•) (Walmsley and Golov 2008a), as well as by an impulsive spin-down to rest (□) 
( |Walmsley et al.||2007a| from 1.5 rad/s and 0.5 rad/s, at T = 1.60 K. All tangles 
were probed by pulses of free ions in the horizontal direction. The line L oc t~ 3 / 2 
corresponds to Eq. (40) with v = 0.2k. 



agreement between the late-time decays. At temperatures 0.7 - 1.6 K the 
late-time decays of quasiclassical turbulences generated by these two different 
techniques were identical within the experimental errors. 



4-2.3. Ultraquantum turbulence generated by a jet of charged vortex rings 

In the temperature range < T < 0.5 K, all tangles produced by a pulse 
of injected current of duration 0.1-1 s, intensity 10~ 12 -10~ 10 A and in the 
driving field — 20 V/cm revealed the late-time decay L oc i -1 , all with the 
same prefactor (Fig. 42) (Walmsley and Golov 2008a). This universality of 
the prefactor for all initial vortex densities is a strong argument in favour 
of the random character of the tangles whose decay is described by Eq. (39) 
(ultraquantum turbulence). During the injection, the tangle originates near 
the injector, presumably as the result of colliding many charged vortex rings 
all of the same radius R ~ 0.5 /mi, and then spreads into the centre of the 
cell in 3-5 s, eventually filling all container and becoming nearly homogeneous 
after ~ 20 s. This can be seen at the transient in Fig. 42 (*) where the tangle 
was initiated at the bottom injector but then sampled along the horizontal 
axis of the container. The dynamics of the tangle spreading was found to be 
independent of the driving field 0-20 V/cm. 

The tangles produced by the bottom and left injectors had very nearly 
identical late-time decays L oc (Fig. 42). The corresponding values of u/k, 



;4 



81 




Fig. 41. Free decay of a tangle produced by a jet of free ions from the bottom 
injector (•) (Walmsley and Golov 2008a| ), as well as by an impulsive spin-down to 



rest dWalmsley et aZ.||2007a[ ) from 1.5 rad/s, at T = 0.85 K. All tangles were probed 



by pulses of free ions in the horizontal direction. The ion jet data are the average of 
nine measurements at each particular time but the spin-down data show individual 
measurements. 

inferred using Eq. (39), are 0.120 ± 0.013 and 0.083 ± 0.004. As the scattering 
diameters a of charged vortex rings produced only by the left injectors could 
be calibrated directly in situ on the arrays of rectilinear vortex lines at steady 
rotation, the absolute value of v = 0.08k for these tangles probed along the 
horizontal axis is treated as more reliable. 



4-3. Discussion: dissipation in different types of turbulence 



4-3.1. Dissipation in 4 He at different temperatures 

Let us now summarize what is known about the dissipation of various 
types of turbulence in superfluid 4 He at different temperatures. As the relevant 
parameter, we plot in Fig. 43 the effective kinematic viscosity v (as denned in 
Eq. (38)) which is a function of the mutual friction parameter a(T) but can 
be different for different types of flow. One can see that at high temperatures, 
T > 1 K, all experimental points group around ~ 0.1 k, give or take a factor 
of two. The fact that different types of experiments seem to produce slightly 
different absolute values and temperature dependences might have various 
reasons. Firstly, experimental techniques rely on the means of calibrating the 
sensitivity to the absolute value of L and on the knowledge of other parameters 
in the model. The calibration is normally performed on an array of rectilinear 
vortex lines in the direction perpendicular to the lines but not on tangles of 
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Fig. 42. Free decay of a tangle produced by beams of charged vortex rings of different 
durations and densities, T = 0.15 K. The injection direction and duration, and 
driving field are indicated. Probing with pulses of charged vortex rings of duration 
0.1-0.3 s were done in the same direction as the initial injection, except in one case 
(*). The line L oc t^ 1 corresponds to Eq. (39) with v = 0.1k. 



vortex lines. For example, recent refined calculations (Chagovets et al. 2007) 
suggested to correct all previous second sound measurements of L on isotropic 
tangles by a factor of 37r/8 ~ 1.2, which we apply here. Secondly, the exact 
values and dependence on temperature of the effective parameters relating 
L(t) to v (such as B m 1.2 in Eq. (39), and C « 1.5 and k\/h m 2tt in 
Eq. (40)) are not known, thus complicating precise comparison of v for different 
temperatures and types of flow. Thirdly, at T > 1 K the normal component can 



also become turbulent 1 


Schwarz and Rozen 


1991. 


Barenghi et al. 


2002 


2000 


L'vov et a/.|2006b 


), hence its velocity field might vary for different 



of generating turbulence. 

Theoretically, one can consider two limiting cases: either completely locked 
(v n = v s ) turbulent flow of both components or an absolutely laminar normal 
component (v n = in a local reference frame) that slows down the segments 
of quantized vortices moving past it. The former is favourable at high mutual 
friction a(T) and low viscosity of the normal component T) n (T) (i. e. at the high- 
temperature end) and is expected to be described by the kinematic viscosity 
r] n / p (dotted line in Fig. 43), where p = p n + p s , while the latter is favourable 
at low a and high r/ n (i.e. at the low-temperature end) and is expected to 
follow the predictions of mutual-friction-controlled models which assume a 
laminar normal component (e.g. line-connected stars and asterisks for random 
tangles in Fig. 43). The experimental situation is obviously somewhere in 
between (and the measured v should probably tend to the smaller value of 
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Fig. 43. The effective kinematic viscosity v{T) for various types of turbulence in 
4 He. (i) Quasiclassical turbulence: v{T) inferred from the L oc i -3//2 dependent free 
decay of tangles using Eq. (40) and produced by impulsive spin-down (open triangles 
Walmsley et~aT\ ( [2007a| )); an ion jet (open circles |Walmsley and Golov| (12008a)); 
towed grids of two different designs (open squares Stalp et al. ( 2002 1 , Niemela et al. 



(2005|). The solid theoretical curve represents Eq. (45) (Vinen and Niemela 2002) 



and the dashed curve is a fit to the Walmsley et al. (2007a) data (open triangles 



and circles) by Kozik and Svistunov ( 2008b |.(ii) Ultraquantum turbulence: v(T) 



inferred using Eq. (39) from L oc t dependent free decay of tangles produced 



by colliding charged vortex rings (filled up- and down- triangles Walmsley and 



squares 



Golov ( 2008a )); counterflow (filled right-pointing triangles Vinen (1957b) and filled 



Schwarz and Rozen (1991 )); ultrasound (filled diamonds 



Milliken et al. 



(1982|); computer simulations (line-connected stars Tsubota et al. (2000)); as well 
as from the analysis of the measured tangle density in applied counterflow (filled 
left-pointing triangles Vinen (1957b])); and computer simulations (line-connected 



asterisks Schwarz (1988)). Note that in simulations (* and *) the normal component 



was artificially clamped to laminar flow, while in most experiments at T > 1 K it 
is involved in turbulent motions to different extent. The values of the kinematic 
viscosities, r? n / p n of the normal component and rj n / (p n + p s ) of the "coupled normal 



and superfluid components" are shown by dotted lines (Donnelly and Barenghi 



1998). 



those predicted by the two limits), although the particular means of driving 
might tip the balance towards either of the limits. This might be one reason 
why different techniques show systematic disagreement in the temperature 
range above 1 K. 

One can speculate that quasiclassical turbulent flow generated by mechan- 
ical stirring is less prone to such uncertainties, since both the superfluid and 
normal large-scale flows are generated simultaneously and are probably locked. 
On the other hand, the generation of quasi-random tangles was attempted by 
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a variety of techniques: those continuously pumped by counterflow in narrow 
channels might have the normal component more laminar, while those with 
freely decaying turbulence in wide channels were shown to evolve towards cou- 



pled superfluid and normal turbulence at high temperatures (Barenghi and 



Skrbek 2007), especially so if created with ultrasound which equally pumps 
turbulence in both components. In any case, on approaching T = 1 K from 
above, when the normal component progressively becomes more viscous and 
less dense, all experimental and numerical values of v(T) for quasi-random 
tangles seem to converge, which supports these reasonings. 

At temperatures below 1 K, the values of v(T) appear to split rapidly: 
those for ultraquantum turbulence apparently stay at nearly the same level, 
~ 0.1 k, as at high temperatures, while those for quasiclassical turbulence keep 
decreasing until the zero-temperature limit of ~ 0.003 k. The mean free path of 
excitations in 4 He rapidly increases with decreasing temperature and becomes 
comparable with i ~ 100 pm at 0.7 K and with the container size h = 4.5 cm 
at 0.5 K. At these temperatures, the normal component is too viscous to 
follow the turbulent superfluid component and hence the convenient model 
of a laminar normal component becomes adequate here (in this respect this 
regime of 4 He is similar to that in superfluid 3 He-B below 0.4 T c discussed 
in previous sections). It is indeed comforting that numerical simulations of 



tangles initiated at short scales (without large-scale turbulence (Tsubota et 



q/.|2000 )) are in good agreement with experiments for random tangles at these 



temperatures (Walmsley and Golov 2008a[ ). As yet, there are no satisfactory 
computer simulations of homogeneous quasiclassical turbulence in the limit of 
zero temperature. 

Between 1.1 - 1.6 K, the values of ^(T) for quasiclassical turbulence gen- 



erated with towed grids of different designs (Stalp et al. 1999, Niemela et al. 



2005), impulsive spin-down to rest, and an ion jet are reasonably consistent 
between each other and show an increase from v = 0.05 k to 0.2 k. There 
was an early attempt (Stalp et al. 1999) to relate these values with the kine- 
matic viscosity rj n / p = 10~ 4 cm 2 /s ~ 0.1 k, assuming locked turbulence of the 



superfluid and normal components. However, it was soon realized (Stalp et 



al. 2002) that the contribution of quantized vortex lines to dissipation con- 



trolled by mutual friction and reconnections is also important. This certainly 
dominates below 1.2 K, where rj n /p begins to increase with cooling (hence, 
making the normal component effectively decoupled from vortices on short 
length scales) while the measured v{T) decreases with cooling. 

In the mutual-friction-dominated regime and provided the normal compo- 
nent is only locked to the superfluid velocity at scales > £, the actual dissipa- 
tion is through the self-induced motion of vortex lines relative to the normal 
component - whatever the kinetics of the quantum cascade. In this regime, a 
segment of a vortex line bent at radius 1Z and hence moving at self-induced 
velocity v ~ BkTZ^ 1 dissipates energy at the rate ~ p s otB 2 n i 7l~ 2 per unit 
length, where B « 1.2 is introduced in Eq. (39). This approach leads to the 
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formula (Vinen 2000 Vinen and Niemela 2002) 



v = sB 2 c 2 ,an . 



(45) 



The multiplier s is meant to account for the degree of correlation between the 
motion of individual segments of vortex lines and the surrounding superfluid 
in a structured quasiclassical tangle (s ~ 0.6), compared to a random tangle 
of the same density L (s = 1). In the simulations of counterflow-maintained 



tangles in the local induction approximation by Schwarz ( 1988 ), the measure of 
the relative abundance of small scale kinks on vortex lines, c 2 . = L~ l (TZ~ 2 ), was 
found to increase from 2 to 12 as a decreased from 0.3 to 10 -2 with decreasing 
temperature from 2 to IK. This formula nicely agrees with the simulations 
by Tsubota et al. (2000). It also captures the trend and magnitude of v{T) 



for all experimental quasiclassical tangles as well as for the random tangle 



maintained by counterflow (Vinen 1957a) (where the normal component is 



expected to be nearly laminar); although to fit it to the data of quasiclassical 
tangles in Fig. 43 (solid line) we used s = 0.2, apparently owing to partial 
locking between the flows in the normal and superfluid components. 

As vortex segments with smaller radii of curvature TZ lose their energy 
faster, for a given energy flux down the length scales there is exists a dissi- 



pative scale A q below which the cascade essentially cuts off ( 


Svistunov 


1995 


Vinen 


2000 


L'vov et al. 


2006b 


Kozik and Svistunov 


20081 


i). For a 


devel- 



oped non-li near Kelvin-wave cascade 
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6k oc k 



(A q « 



with the amplitude spectrum 
( |Kozik and Svistunov||2004[ ), this implies that (7^ 2 ) ~ £~ 2/5 A q 8/5 . 
In this sense the Schwarz's parameter c 2 . ~ (£/A q ) 8//5 quantifies the range of 
wavelengths involved in the quantum cascade and becomes another quantum 
analog of the classical Reynolds number: Re = (/i/A c i) 4 / 3 . With decreasing 
a, the value of A q progressively decreases until, at T < 0.5 K (Vinen 2001 



L'vov et al. 2006b, Kozik and Svistunov 2008b), it reaches the wavelength 



A p h ~ 3 nm at which phonons can be effectively emitted. Below this tempera- 
ture, as the dissipation of Kelvin waves due to mutual friction becomes negligi- 
ble compared to the emission of phonons, v{T) is independent of temperature. 
However, in the presence of the quantum cascade at larger wavelengths, it is 
the cascade's kinetics that controls the energy flux, and it might happen that 
the vortex length L already saturates at higher values of a, thus signalling the 
onset of the regime with v{T) = const as T — > 0, Eq. (38). 

Let us now discuss the specific models for the low-temperature behaviour 
of ^(T). For the dissipative mechanisms that involve individual quantized vor- 
tices, such as those maintained by vortex reconnections, simple considerations 



(Vinen and Niemela 2002) usually predict v to be of order k. However, more 



involved models, such as that of a "bottleneck" for quasiclassical turbulence, 



might suggest smaller values (L'vov et al. 2007a). If certain configurations of 
vortices are inefficient in transferring energy down to short wavelengths {e.g. 
either there exists a non-cascading part of the spectrum or reconnections are 
less efficient in polarized tangles), the resulting accumulation of an extra con- 
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tribution to L leads to a reduction of the parameter v. This is what is observed 
for quasiclassical tangles in 4 He on approaching the T = limit. 

An important question remains: how is the energy of classical eddies passed 
over to the shorter quantum length scales? Two main mechanisms are cur- 
rently discussed: (i) the excitation of Kelvin waves through purely non-linear 
interactions in classical eddies or (ii) vortex reconnections (each leaving a 
sharp kink on both vortex lines, hence effectively redistributing the energy to 
smaller length scales). Both predict an increase in L(T) at constant energy 
flux down the cascade (i.e. the decrease of f(T)) with decreasing temperature. 
The former considers the accumulation of extra "non-cascading" quasiclassical 
vorticity at length scales above I due to the difficulty in transferring energy 
through wavenumbers around l~ x if reconnections do not ease the process (re- 
ferred to as the "bottleneck" model in Sec. 3) (L'vov et al. 2007a). The latter 
mainly associates the excess vortex length L with the contributions from the 
new self-similar structures produced by vortex reconnections on length scales 
shorter than £, when vortex motion becomes progressively less damped at low 
temperatures (Kozik and Svistunov 2008b). This model relies on four types 
of processes: reconnections of vortex bundles, reconnections of neighbouring 
vortices, self-reconnections of a vortex, and non-linear interactions of Kelvin 
waves, which bridge the energy cascade from the Kolmogorov to the Kelvin- 
wave regimes. It can be successfully fitted to the experimental v(T), as seen 
in Fig. 43. In short, both models predict an enhancement of vortex densities 
but on different sides of the cross-over scale I. 

Finally we may ask why is u(T = 0) larger for ultraquantum (random) tan- 
gles than for quasiclassical tangles in the T = limit? Whatever the model, 
this means that, for the same total density of vortex lines L, the rate of energy 



dissipation in random tangles is larger. Within the bottleneck model (L'vov et 



aI1 |2007a[ ), this seems straightforward as these tangles have no energy on clas- 



sical scales at all. Within the reconnections/fractalization scenario (Kozik and 



Svistunov 2008b), this is explained by the fact that reconnections in partially 



polarized tangles are less frequent and less efficient. In the framework of this 
latter model, to explain why there is apparently no temperature dependence of 
v{T) for ultraquantum tangles below T = 1 K (although no experimental data 
exist for the range 0.5 - 1.3 K), one should assume that for this type of tangle 
all significant increase in the total line length L caused by the fractalization 
with decreasing temperature has already occurred above 1 K. 



4-3.2. Comparing turbulent dynamics in 3 He-B and A He 

It would be instructive to provide a comparison of turbulent dissipation 
in superfluid 4 He and 3 He-B. Unfortunately at present, with only few mea- 
surements on varying types of flow, this remains a task for the future. The 
Lancaster measurements on the decay of an in homogeneous tangle generated 
with a vibrating grid by Bradley et al. (2006) revealed a L oc t~ 3 ^ 2 late-time 
decay (Fisher and Pickett 2008). Assuming that the origin of the r 3/2 de- 
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Fig. 44. Normalized measures of dissipation rate, the effective kinematic viscosity 
v and the normalized velocity of the propagating vortex front Vf = Vf/CIR, as a 
function of mutual friction dissipation a in 4 He and 3 He-B: (i) v(T) for quasiclassical 
turbulence in 4 He: spin-down, open triangles ( Walmsley et al.|2007a ); ion jet, open 
circles (Walmsley and Golov||2008a ); towed grid, open diamonds ( Stalp et aZ.||2002 



Niemela et al. 20051; (ii) v(T) for ultraquantum turbulence in He: colliding charged 



vortex rings, filled triangles (Walmsley and Golov 2008a I; numerical simulations, 
line-connected stars ( Tsubota et aZ.||2000 ) , and asterisks ( Schwarz||1988). (iii) v(T ) 
estimate for ultraquantum turbulence in 3 He-B: dashed line (Bradley et al. 20061. 
(iv) V{ in 3 He-B: experiment, open squares (Eltsov et al. 20071; and numerical 
simulations, line-connected squares (Sec. 3). 



pendence is the same as in quasiclassical turbulence with a saturated energy- 
containing length, and that this length is equal to the spread of the turbulent 
tangle over a distance of 1.5 mm, using Eq. (40) they extract the temperature- 
independent value of v « 0.2k at 0.156 < T/T c < 0.2. This is different from 
the 4 He value v ~ 0.003 k. However, it is questionable whether for this tur- 
bulence the condition for the applicability of Eq. (40) is met, namely that 
the classical energy S c dominates over the quantum energy S q . From their re- 
cent measurements (Bradley et al. 2008) of the initial velocity of mean large- 
scale flow, u ~ 0.5mm/s, and vortex density, L ~ 10 4 cm -2 , we estimate 
£ c ~ u 2 /2 ~ 10- 3 cm 2 /s while S q ~ 0.7k 2 L ~ 10" 2 cm 2 /s; i.e. £ c < S q even 
in the early stages of the decay. Hence, it is possible that the dynamics of 
the decay of their localized tangles is that of an ultraquantum turbulence but 
accelerated due to diffusion or emission of vortex rings into space; thus the 
decay L(t) is steeper than oc t~ x . More experiments with grid turbulence in 
3 He-B are thus desirable, which could shed light on the temperature depen- 
dence of turbulent dynamics, e.g. if there is a similar decrease in z/(T) below 
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a ~ 10~ 3 as measured for 4 He, or to provide an independent characterization 
of the energy-containing length scale. It would also be important to reproduce 
the results for 4 He in the T = limit with another method of generating 
turbulence (e.g. by means of towed or vibrating grids). 

The Helsinki experiments on the propagating vortex front (Sec. 3) provide 
a measure of dissipation via the front velocity Vf. The observed relative en- 
hancement of the rate of decay below 0.4 T c (a < 1) proves the efficiency of the 
turbulent cascade in dissipating the energy of large-scale flow. The leveling off 
of this rate below 0.28 T c (a < 10~ 2 ) is qualitatively similar to that observed in 
4 He at T < 0.5 K (a < 10~ 5 ), and might hint at the T = limit in the inertial 
cascade. However, the respective values of a, at which the leveling-off occurs, 
differ substantially. Various reasons can be given to explain this difference, eg. 
that different types of flow have their own specific u(a) dependence (evidence 
for this is the difference in the ^(T) values of quasiclassical and random tan- 
gles). For instance, to test the interpretation of the temperature dependence 
of Vf(T) on approaching the T — ► limit, given in Sec. 3.5.6 in terms of the 
proximity of the outer and quantum crossover scales, it would be instructive 
to conduct similar measurements in containers of different radii. 

In front propagation the energy dissipation rate is proportional to Vf(T), 
while in experiments described in this section the rate is proportional to v(T). 
Thus it is interesting to compare the temperature dependences of Vf(T) and 
v(T) in the two superfluids. We plot in Fig. 44 the effective kinematic viscos- 
ity v/k and the normalized velocity of the propagating front Vf = Vf/QR as 
a function of mutual friction dissipation a. Note that the role of turbulence 
comes into play in a very different manner in these two types of flow: as nearly 
homogeneous and isotropic turbulence in 4 He in the Manchester experiments 
and, roughly speaking, via the deviation from laminar vortex front propaga- 
tion, owing to turbulent excitations, in 3 He-B in the Helsinki experiments. 

It is striking that for a < 0.1 (i.e. the regime where developed turbulence 
becomes possible) and before leveling off at T — > 0, the data for all experi- 
ments and numerical simulations discussed here approach a similar slope of 
. For v(a), in the spirit of Eq. (45), this implies that the parameter c| is 



a 



0.5 



roughly proportional to a - even well below the values of a of the calcula- 
tions by Schwarz (1988). On a finer scale, particular models of matching the 
classical and quantum cascades, as discussed in the previous subsection, might 
produce specific dependences of 02(a), which perhaps can be tested in future 
experiments in more detail, provided that the quality of data is improved. 

We do not attempt to compare the absolute values of vj k and V{ in Fig. 44. 
Instead, we concentrate on the a— dependences, especially on the values of a 
at which the dissipation rates level off as T — ► 0. The fact that the measured 
3 He-B data level off at higher values of a might suggest that the quantum 
cascade in 3 He-B is not as developed as in 4 He and that the ultimate dissipative 
mechanism, independent of a in the a — > limit, is stronger and takes over 
at larger length scales in 3 He-B. Indeed, the numerically calculated velocity 
of front propagation ff(a), which does not explicitly incorporate any excess 
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dissipation beyond mutual friction dissipation and hence is equally suitable 
for 3 He-B and 4 He, is steadily decreasing with decreasing a, following the 
same trend as the experimental and numerical v(T) for 4 He. In contrast, the 
measured Vf(ct) levels off below a ~ 1 • 10~ 2 , which seems to imply that the 
T = regime prevails at a < 1 • 10~ 2 in 3 He-B, while this happens only below 
a < 1 • 10~ 5 in 4 He. The fermionic nature, the large vortex core diameter, and 
the much lower absolute temperatures lead to different inherent dissipation 
mechanisms in 3 He-B. For example, the energy loss during each reconnection 
event can be substantial in 3 He-B but not in 4 He. This is supported by the 
fact that the values of v measured in Lancaster for 3 He-B with a ~ 30 nm 
(at P = 12 bar) level off (i.e. branch away from the common trend of v(a) 
marked by the numerically calculated dependences) at larger values of a ~ 
3 • 1CT 2 , than the measured front velocity V{ for a ~ 16 nm (at P = 29 bar) 
at a ~ 1 • 1CT 2 . We conclude by noting that further work on turbulence in 
different types of flow in both 4 He and 3 He-B is highly desired. 

4-4- Summary: decay of quasiclassical and ultraquantum turbulence 

The newly developed technique of measuring the density L of a vortex tan- 
gle by the scattering of charged vortex rings of convenient radius ~ 1 /im 
has made it possible to monitor the evolution of tangles in superfluid 4 He 
down to below 0.5 K, i.e. deep into the zero-temperature limit. The dynamics 
of two very different types of tangles can be studied. Quasiclassical tangles, 
mimicking the flow of classical liquids on large length scales were generated 
by an impulsive spin-down from angular velocity fl to rest of a rotating cubi- 
cal container with helium. After a transient of duration ~ 100 /Q the turbu- 
lence becomes nearly isotropic and homogeneous and decays as L oc t~ 3 / 2 , as 
expected for turbulence posessing a Kolmogorov cascade of energy from the 
energy-containing eddies of constant size, set by the container size h. This was 
studied in isotopically pure 4 He in a broad range of temperatures of 80 mK - 
1.6 K, corresponding to the range of mutual friction a from ~ 10~ 10 to 10 _1 . 
Identical results were obtained for the free decay of quasiclassical tangles gen- 
erated by a central jet of ions, although only at temperatures above 0.7 K so 
far. 

Alternatively, non-structured (ultraquantum) tangles of quantized vortices, 
that have little flow at scales above the inter- vortex distance and hence no clas- 
sical analog, can be obtained by colliding many small quantized vortex rings of 
radius R <C £ at temperatures below 0.5 K, i.e. in the zero-temperature limit. 
These tangles took about 10 s to spread from an injector into all experimental 
volume of size h — 4.5 cm, followed by free decay with the universal dynam- 
ics L oc t -1 , independent of the initial conditions. The relatively fast rate 
of spreading is indeed surprising and might be due to the small polarization 
(mean velocity) of the tangle. 

Quantitative measurements of the free decay for both types of tangles al- 
low the extraction of the "effective kinematic viscosities" v - the flow-specific 
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parameter linking the rate of energy dissipation with the total density of vor- 
tices L. It turned out that at high temperatures T > 1 K, where the quantum 
cascade is not well-developed, the values of v for both types of turbulence 
are comparable. However, in the zero-temperature limit, where the dissipa- 
tion can only take place at very short length scales <C I, to which the energy 
can only be delivered by a cascade of non-linear Kelvin waves on individual 
vortex lines, the saturated value of v for quasiclassical turbulence seems to 
be smaller than that for ultraquantum turbulence. On the microscopic level 
this is most probably related to the fact that in the presence of quasiclassical 
eddies the partial mutual alignment of vortex lines in "bundles" slows down 
the process of exciting Kelvin waves for a given vortex density L. For instance, 
this mutual alignment is expected to reduce the frequency and efficiency of 
vortex reconnections, which are believed to be the defining process of the dy- 
namics of non-structured (ultraquantum) tangles as they fuel the Kelvin-wave 
cascade. At present, two microscopic models are developed to describe the 
energy transfer from classical to quantum scales (cascades) in quasiclassical 
tangles: one relying on reconnections and the other assuming that additional 
"non-cascading " vortex density helps maintain the continuity of the energy 
flux down to the dissipative length scale. 

We have also compared the rates of dissipation as a function of mutual 
friction in superfluid 4 He and 3 He-B. Both similarities and differences in the 
behaviour are spotted and discussed. However, a final quantitative verdict on 
the role of specific mechanisms cannot be delivered at this stage, pending 
further development of theoretical models and of experiments on different 
types of turbulent flow in both superfluids. 

5. Summary 

In the last five years we have witnessed important advances in the under- 
standing of the appearance, growth, and decay of different types of superfluid 
turbulence, especially in the fundamentally important limit of zero tempera- 
ture where the intrinsic processes within the superfluid component set the dy- 
namics of the tangle of quantized vortices. As these are absolutely undamped 
on a broad range of length scales down to wavelengths of order 10 nm - 1 /im, 
depending on the type of superfluid, a principally new microscopic dynamics 
emerges - set by the instabilities of individual vortices, their reconnections 
and one-dimensional cascades of energy from non-linear Kelvin waves (wave 
turbulence) along individual vortex lines. Amazingly, even in this limit, var- 
ious turbulent flows behave classically on large length scales, and often the 
observed rate of decay is no different from that at high temperatures - owing 
to the fact that the energy-containing and dissipative lengths and times are 
well-separated. 

However, also the discrete nature of quantized superfluid vorticity becomes 
important in such processes as the growth of an initial seed vortex tangle out of 
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a single- vortex instability, the efficiency of dissipating the energy of large-scale 
eddies through short-wavelength Kelvin waves on individual vortex lines, or 
the dynamics of non-structured tangles possessing no large-scale flow. In the 
helium superfluids these underlying processes are separately observable and 
are ultimately expected to become the corner stones of a detailed theoretical 
framework. This should make it possible to develop a consistent picture of 
turbulence in helium superfluids, which describes the nonlinear turbulent dy- 
namics of discrete line vortices in a macroscopically coherent quantum liquid 
of zero viscosity. 
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